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The modern formulation of exclusive reactions within Quantum Chromody- 
^ , namics is reviewed, the emphasis being placed on the pivotal ideas and meth- 

\ ods pertaining to perturbative and non-perturbative topics. Specific prob- 

' lems, related to scale locality, infrared safety, gluonic radiative corrections 

(Sudakov effects), and the role of hadronic size effects (intrinsic transverse 



momentum), are studied. These issues are more precisely analyzed in terms 
of the essential mechanisms of momentum transfer to a hadron while remain- 
ing intact. Different factorization schemes are considered and the conceptual 
lacunas are pointed out. The quite technical subject of renormalization-group 
Q^' evolution is given a detailed account. By combining analytical and numerical 

1^ ■ algorithms, the one-gluon exchange nucleon evolution equation is diagonalized 

and next-to-leading eigenfunctions are calculated in terms of Appell polyno- 
mials. The corresponding anomalous dimensions of trilinear quark operators 
^ ' are found to form a degenerate system whose envelope shows logarithmic 

' large-order behavior. Selected applications of this framework are presented, 

focusing on the helicity-conserving elastic form factors of the pion and the nu- 
cleon. The theoretical constraints imposed by QCD sum rules on the moments 
of nucleon distribution amplitudes are used to determine a whole spectrum of 
optional solutions. They organize themselves along an "orbit" characterized 
by a striking scaling relation between the form- factor ratio R = \Gy^\/G^ 
and the projection coefficient i?4 on to the corresponding eigensolution. The 
main reasons for the failure of the present theoretical predictions to match 
the experimental data are discussed and workable explanations are sketched. 

ll.lO.Hi, 12.38.Aw, 12.38.Bx, 12.38.Cy, 12.38Lg, 13.40.Gp, 13.40Hq, 14.40.Aq 

Typeset using REVT^ 



*E-mail: stefanis@tp2.ruhr-uni-bochum.de 



1 



Contents 



I PREFACE 

ig INTRODUCTION 

[mi STANDARD CONVOLUTION SCHEME 

FACTORIZATION 

g SHORT DISTANCE PART 

LARGE DISTANCE PART: PERTURBATIVE ASPECTS . . . 

[1| Renormalization 

Meson evolution equation 

1^ Nucleon evolution equation 

LARGE DISTANCE PART: NON-PERTURBATIVE ASPECTS 

NUCLEON DISTRIBUTION AMPLITUDES 

DELTA DISTRIBUTION AMPLITUDES 



ig ELECTROWEAK FORM FACTORS 

NUCLEON FORM FACTORS 

g NUCLEON-DELTA TRANSITION FORM FACTOR 

CHARMONIUM DECAYS 



V| MODIFIED CONVOLUTION SCHEME 

MODIFIED FACTORIZATION 

g MODIFIED NUCLEON FORM FACTOR 

Screening of as singularities 

1^ Numerical form-factor analysis 

1^ Higher-order nucleon distribution amplitudes . . 

11 Global pattern of nucleon distribution amplitudes 



VII CONCLUSIONS AND OUTLOOK 



APPENDIXES 



3 EIGENFUNCTIONS OF THE NUCLEON EVOLUTION EQUA 
TION 



[B] ELASTIC FORM FACTORS OF THE NUCLEON 



I. PREFACE 



Quantum Chromo dynamics (QCD) is that part of the Standard Model which is sup- 
posed to describe the strong interactions at the microscopic level. However, despite the 
numerous phenomenological successes of this theory in regimes where perturbation theory 
applies, several conceptual and technical challenges concerning the large-distance domain 
still remain. 

Understanding the questions and fixing the problems is possibly not enough to yield 
convincing quantitative answers towards an analytic description of confinement and the 
formation of hadronic bound states. Especially the theoretical analysis of exclusive processes, 
in which mtac^ hadrons appear in the initial and final states, involves the detailed calculation 
of hadronic wave functions - incalculable within a pcrturbativc framework. Nevertheless, 
it is tempting to analyze the current status of knowledge by exposing physical scenarios 
pertinent to understanding some basic features in the transition from the perturbative to 
the non-perturbative phase. Probably the greatest technical barrier here is the mutation of 
light current quarks - the degrees of freedom of high-energy QCD - to massive constituent 
quarks, operative after gluon bosonization in low-energy effective theories. 

By now the pure perturbative QCD option to large-momentum transfer exclusive pro- 
cesses has been proved elusive. The troubles may be a refiexion of the poorly known complex 
substructure of hadrons (intrinsically linked to the non-perturbative QCD vacuum) ; or the 
faults may lie elsewhere, e.g., in the limited knowledge of higher-order perturbative correc- 
tions and the possible lack of convergence of the perturbative expansion. A common key 
element of these difficulties is infrared (IR) sensitivity due to incomplete factorization of 
short from large-distance effects. As a consequence, either there are intrusions in forbid- 
den kinematical regimes - invalidating the initial factorization assumptions - or in order to 
preserve IR finiteness, severe IR-cutoff prescriptions have to be employed - resulting in a 
depletion of the perturbative contributions. 

This report focuses on these problems from a user's perspective, identifying the principal 
drawbacks, and trying to set benchmarks for future developments. We concentrate on two 
major subjects: (i) hadron distribution amplitudes and elastic form factors, and (ii) imple- 
mentation of the renormalization group within appropriate factorization schemes. To set 
the stage, we discuss and review problems relating to the momentum (mass) scales involved 
in form factor calculations: scale locality, IR safety, gluonic radiative corrections (Siidakov- 
type form factors), and the role of hadronic size effects (intrinsic transverse momenta). In 
addition, we give a systematic non-perturbative analysis for determining hadron distribu- 
tion amplitudes in conjunction with QCD sum rules. These issues are related to distinct 
mechanisms by which a large space-like momentum is transferred to an intact hadron and we 
use detailed calculations to investigate how these mechanisms influence the predictions for 
-^7r(Q^), G^m(Q^)) ^m{Q^)j G^((5^) relative to existing data. Other phenomenological 
applications in this context are also discussed. 

II. INTRODUCTION 

Lack of precise knowledge about the large-scale inter-quark forces and the mathematical 
intractability of the nonlinear structure of the QCD dynamics preclude an ab initio analytical 
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Figure 1. Perception of strong-interaction dynamics with momentum flow. 



calculation of quantities like the hadron wave functions. Of the various approaches now being 
studied, those employing quark/gluon condensates, the order parameters of non-perturbative 
QCD, are among the most successful - albeit not without their loose ends. In particular, 
nonlocal condensates which ascribe a nonzero average virtuality to vacuum quarks and 
gluons (a technical review with earlier references is given in f^) may enable one to do more 
realistic calculations of hadron wave functions, but such techniques are still in an embryonic 
phase of development. Thus far the non-locality of quark/gluon condensates appears only 
in model form, though it was found by lattice calculations that it decreases exponentially 
with the coordinates. 

Frameworks which are capable of simulating the non-perturbative regime of QCD, either 
in the continuum (e.g., via QCD sum rules or numerically on the lattice 0], have been 
widely used during the last decade. Although modern lattice calculations in the strong- 
coupling regime are very promising, there are conceptual and technical limitations: e.g., 
the properties of a single proton cannot be determined properly at present. Especially the 
computation of higher-order moments (corresponding to higher derivatives acting on local 
operators) of light-cone distribution amplitudes appears very difficult at present On the 
other hand, since hadron wave functions enter only integrated quantities, like form factors, 
their extraction directly from the data is still in its infancy, and much more experimental 
input is needed to make significant progress here. 

In the past few years, there have been several theoretical attempts to emulate low- 
energy QCD by using effective chiral (x) actions, e.g., the Nambu-Jona-Lasinio model (for a 
recent review, see, e.g., and references cited therein). Such approaches implement chiral 
symmetry breaking and are, in principle, deducible from the instanton model of the QCD 
vacuum [Q, albeit a rigorous derivation is still rudimentary, and the range of their validity 
is limited to momentum values below an ultraviolet (UV) cutoff scale of about 1 GeV. 

The perception of strong-interaction dynamics at different resolution (energy) scales is 
illustrated in Fig. 0. 

Exclusive processes in QCD are therefore particularly interesting because they provide 
the possibility of analyzing hadron wave functions in terms of their quark and gluon de- 
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grees of freedom. In order to take advantage of perturbative QCD in exclusive processes,! 
a short- distance amplitude has to be isolated and proved to be IR finite. This provides a 
useful handle on strong-interaction amplitudes and one is in the position to make quanti- 
tative computations in a systematic fashion. Formally, this corresponds to the procedure 
of multiplicative separation (factorization) of regimes, ubiquitous in many areas of theoret- 
ical physics. It is the property of cross sections (in inclusive processes) or amplitudes (in 
exclusive processes) that high-momenta ( "hard" ) and low-momenta ( "soft" ) regions can be 
disentangled in such a way that the factorized parts depend only on the dynamics specific 
for the corresponding scale. 

Proofs of factorization theorems both for inclusive cross sections and exclusive 

amplitudes, have reached a certain level of sophistication, but the subject is still open (for 
a comprehensive review, we refer to lll|]). For instance, the generalization of factorization 
theorems beyond leading twist faces severe problems due to incomplete Bloch-Nordsieck 
(BN) cancellation of IR divergences |T^,n]. Along similar lines of argument, uncanceled 



IR divergences near the boundary of phase space (the so-called endpoint region) have to 
be re-summed by improving renormalization group techniques to render the perturbative 
calculation sound. Completing this procedure, the second step is to interphase the factorized 
parts to renormalization. 

According to these ideas, the asymptotic, i.e., large Q^-behavior of electroweak form 
factors can be calculated with perturbative QCD by simply assuming free valence quarks 
entering and exiting the hard-scattering region, without any recourse to confinement. It 
is clear that as the external momentum becomes smaller, the resolution scale decreases, 
quantum modes with corresponding wavelengths emerge and the pure perturbative picture 
of quasi- free quarks bound together by single-gluon exchange breaks down. At this point, 
fluctuations of the background fields of the non-trivial QCD vacuum interfere and a self- 
consistent computation must take them into account. Where this transition of the pure 
perturbative phase to the non-perturbative one takes place is a point of ongoing controversy 

A useful framework for describing exclusive processes along these premises is the convo- 
lution scheme of Brodsky and Lepage |T9| , |20| , |2l| . [A somewhat different approach to hard 
processes was developed independently by Efremov and Radyushkin |ll6| , ^ , p3| .] Within 
this scheme, the reaction amplitude becomes the product (the convolution) of two or more 
factors, each depending only on the dynamics specific for that particular momentum (or dis- 
tance scale), and the evolution of the factorized parts is renormalization group controlled. 
More important, the hard part of the process becomes amenable to the methodology of QCD 
perturbation theory. To enforce IR finiteness, subtractions of IR singularities are imposed. 
The guiding idea is that the time scale involved in the hard part of the amplitude and that 
one for the formation of the intact hadron in the final state(s) are disparate, so that the 
corresponding dynamics are uncorrelated (impulse approximation). As a consequence of the 
sudden character of the hard interactions, the perturbed quantum mechanical state remains 
unchanged and time evolution is governed by the perturbed Hamiltonian. 



Exclusive" means that the momenta of all individual hadrons participating in the process are 
measured. 
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Figure 2. Inclusive versus exclusive reactions in QCD. 



In leading order of the coupling constant, large- momentum transfer quark-gluon subpro- 
cesses can be adequately described by one-gluon exchange kernels, the justification being 
provided by asymptotic freedom |^^. We will call this approach in the following "the stan- 
dard convolution scheme" of exclusive processes. Detailed applications of this type of QCD 
analysis to hadronic form factors and decay amplitudes are given in later sections of this 
survey. 

Having extracted a (process-dependent) hard-scattering amplitude, the remaining soft 
contributions responsible for the bound-state dynamics are encapsulated in universal (but 
factorization- scheme dependent hadron distribution amplitudes. In axial gauges (e.g., 
the light-cone gauge, = 0), the distribution amplitude is the probability amplitude for 
the hadron to consist of valence quarks with (longitudinal) momentum fractions < Xj < 1, 
Y^iXi = 1 (in an infinite momentum frame) moving coUinearly up to the factorization scale. 
Restricting the flow of soft transverse momenta into the distribution amplitudes amounts 
to subtracting their high-momentum tales, and is tantamount to avoid double counting and 
ensure that all vertices and propagators entering the microscopic processes of a Feynman 
graph are entirely governed by perturbative QCD.i The practical utility of this picture, be- 
sides being physically appealing, resides in the fact that presciently dissecting the reaction 
amplitude in the form of a convolution of soft and hard parts, its evolution behavior, i.e., 
the variation of the parts with momentum, is completely controlled by the large scale of the 
process. In fact, if the momentum transfer is very large, the UV regularization scale of 
the renormalization group can be safely traded for to give rise to renormalization group 
evolution. As the momentum transfer imparted to the hadron by the external electromag- 
netic probe increases, inclusive channels are gradually triggered, and the probability that 
the hadron rebounds without going apart diminishes. This behavior is illustrated in Fig. |^. 

A theoretical tool to calculate hadron wave functions is provided by the operator prod- 
uct expansion [^. In the operator product expansion language, the hard-scattering am- 
plitude (small-distance interaction) corresponds to the Wilson coefficient function, whereas 



more technical discussion of these points is postponed to a following section. 
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the initial- and final-state hadron wave functions (large-distance interactions) derive from 
the matrix elements of local operators with appropriate quantum numbers for each channel. 
A complicating factor in the operator product expansion analysis is that the identification 
of a short-distance part in the exclusive amplitude is not a priori obvious and has to be 
determined process by process. 

At leading-twist level, the soft parts of the factorized exclusive amplitude are represented 
by the valence-state wave functions on the light cone and describe the distribution in lon- 
gitudinal momentum fractions, being averaged over transverse momenta up to values of the 
factorization scale. Contributions from higher Fock states with additional gg-pairs and glu- 
ons (higher twists) are suppressed by powers of the momentum transfer. Such contributions 
correspond in the operator product expansion to composite quark-gluon operators with pro- 
gressively higher dimensions, the inclusion of which introduces new order parameters that 
have to be estimated. 

Applications of QCD sum rules to exclusive processes started already in the early eighties 
by computing the meson and nucleon form factors |2^j2^, as well as the first moment (i.e., 
the decay constant) of the pion distribution amplitude. Following a different strategy, 
V. L. Chernyak, A. R. Zhitnitsky, and I.R. Zhitnitsky (CZ) attempted to reconstruct model 
distribution amplitudes from their few first moments for the pion |2^, the nucleon ||30| | , 



and other hadrons |^T|. The values of the moments were restricted in their approach by 
constraints extracted from QCD sum rules, and model distribution amplitudes as polyno- 
mials in the longitudinal momenta of the valence quarks were derived by means of moment 
inversion. Fixing the second and fourth moment of the pion distribution amplitude, CZ 
proposed a "double-humped" model which brings the value of the pion form factor close to 
the data for a value of the strong coupling constant as as low as 0.3. This was considered 
a remarkable success because the asymptotic prediction yields the value Q^F^^^Q"^) ^ 0.13 
which falls short a factor of three compared to the experimental data |3^ (though the accu- 
racy of the present data is rather poor to be conclusive). A similar analysis for the nucleon 
on the basis of the first- and second-order moments, led the same authors to propose a 
nucleon distribution amplitude which shows considerable asymmetry in the distribution of 
the longitudinal momentum of the valence quarks. 

Soon an alternative nucleon distribution amplitude was suggested by Gari and Stefanis 
(GS) [^, constructed with the aim to yield helicity-conserving nucleon form factors which 
account for the possibility that the electron-neutron differential cross section is dominated 
by Gg, while is asymptotically small or equivalently that there is a sizeable neutron 
Pauli form factor overwhelming the Dirac one at all values. The GS model gives very 
good agreement with the latest high-Q^ SLAC data |]3^ on (or Ff ) and makes realistic 
predictions for the corresponding neutron form factor in the high-momentum region |^ . 



In contrast, as effected in the CZ model overestimates both form factors almost by 
a factor of two in the region of 10-20 GeV^, if realistic values of Aqcd around 200 MeV 
are used.i But, on the theoretical side, a heavy price is paid: some moments of the GS 
model distribution amplitude cannot match the requirements set by the CZ moment sum 



^In the original CZ analysis, the value Aqcd = 100 MeV was used. Such a low value is now de 
facto excluded by experiment. 
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rules in the allowed saturation range ||35|. Moreover, as it was shown later by Chernyak and 



coworkers [36|, this model leads to a prediction for the ^Si pp decay width of charmonium 
which is several orders of magnitude smaller than the experimental value.0 

In 1987 the moment sum rules were re-evaluated by King and Sachrajda (KS) 
spotting the gaps in the CZ analysis and shifting the range of the moment sum rules, albeit 
the gross features of the method were confirmed as was the basic shape of the nucleon 
distribution amplitude. 

A couple of years later, Chernyak, Ogloblin, and I. R. Zhitnitsky (COZ) refined 
their previous moment sum rules for the second-order moments and extended their method 
to third-order moments. Their new moment sum rules comprise 18 terms with restricted 
margins of uncertainty relative to the previous CZ analysis and, in general, comply with the 
results of the KS computation, but contradict those obtained on the lattice for the lowest two 



moments ||3^ . The same authors have also proposed a new model distribution amplitude for 
the nucleon - still restricted to polynomials of second degree in the longitudinal momentum 
- which satisfies all, but 6 of the new moment sum rules, whereas the CZ amplitude and 
the GS one violate, respectively, 13 and 14 of them. The KS amplitude provides almost 
the same quality as that of COZ, with only 7 broken moment sum rules. In the same year, 
Schafer presented a variety of model distribution amplitudes for the nucleon, which 
incorporate polynomials of degree three in the longitudinal momentum, and found that such 
contributions play a token role, if properly incorporated. 

The essence of these investigations is that distribution amplitudes extracted from moment 
sum rules are much broader than the asymptotic solution (to the evolution equation) and 
have a rich structure that is reflected in a quite asymmetric balance in the distribution of lon- 
gitudinal momentum fractions among the valence quarks, accentuated by nodes. The latter 
may be understood in the following way: A valence quark embedded in the non-perturbative 
QCD vacuum rebounds differently from one in "free space" . Hence, quark/gluon condensates 
act somewhat like a dispersive medium in which different components of a wave propagate 



with different speeds and tend to change phase with respect to each other [^. Regions of 



different phases in the nucleon distribution amplitude may be interpreted as evidence for 
binding effects inside the nucleon, amounting perhaps to diquark formation. A full-fledged 
discussion of these issues will be given in the course of this work. 

The next major step in the determination of nucleon distribution amplitudes was done in 
the early nineties by Stefanis and Bergmann with the invention of the heterotic conception 
PT| , |i^ (Fig. ^.i Previously, the CZ model (or its descendants: KS and COZ models) on 
one hand and the GS model on the other hand were treated in the literature as competing 
alternatives - mutually excluding each other. In addition, either way, it was not possible to 
reconcile the theoretical moment sum rules constraints with the experimental data because 
none of these models is able to give, simultaneously, a quantitatively satisfactory agreement 
with the form-factor and charmonium-decay data. 

Perhaps understandably, in view of such distinctive models and predictions, the con- 



^Provided one uses again the favored value of the strong coupling constant = 0.3. 
^ "Heterosis" in Greek means increased vigor due to cross-breeding. 
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Figure 3. Heterotic conception of tlie nucleon distribution amplitude. The heterotic model amalgamates 
typical characteristics of COZ-like and GS-like amplitudes. 



ventional view has been one of fragmentation. But the new idea, underlying the nucleon 
heterotic model, makes it possible to amalgamate the best features of COZ-type and GS- 
type distribution amplitudes into a single mold, thus lifting the disparity between theory 
and experiment. As it will become more transparent in the following sections, the heterotic 
distribution amplitude is a "hybrid" - sort of - and seems to have a foot in each of the 
previous models. This duality is also reflected in its profile which, though distinctive in 
overall shape from both the COZ and the GS distribution amplitudes, bears geometrical 
characteristics typical for both models. 

Similar ideas of heterosis were then applied by the same authors |^ to the distribution 
amplitude of the A"*" (1232) isobar, treating the moment sum rules of Farrar et al. (FZOZ) 
Q and those by Carlson and Poor (CP) |^ in combination. Again a heterotic distribution 
amplitude was determined in between the CP and the FZOZ model distribution amplitudes. 
Using the heterotic distribution amplitudes for the nucleon and the A, the transition form 



factor was calculated within the standard convolution scheme and remarkable 



agreement with the available data was found. Even more, a recent reanalysis by Stuart 
et al. PB|,E^, of the inclusive e — p data in the A(1232) region by the SLAG experiment 
NEll, combined with low data, high data from the SLAG experiment E133, and 
missing mass squared data, finds results for the transition form factor systematically higher 
then the previous data analysis by Stoler |58] and confirms within the errors the "heterotic" 
predictions |^ . 

These developments gave rise to a more general picture (Stefanis and Bergmann) 
49| , |50|j5T| , |52[| , which deals with global features of nucleon distribution amplitudes. The cru- 



cial question is: If it is possible to synthesize one hybrid-type distribution amplitude, namely 
the heterotic one, is it then possible to find also other distribution amplitudes with elements 
belonging to different classes of solutions to the moment sum rules? That the moment sum 
rules are not stringent enough to fix the shape of the nucleon distribution amplitude uniquely 
was already pointed out by Gari and Stefanis , and in more mathematical detail by Ste- 
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fanis A similar analysis for the pion distribution amplitude was given by Mikhailov 

and Radyuslikin [Q. 

To determine the possible variation of distribution amplitudes allowed by moment sum 
rules, one has first to identify appropriate order parameters in order to classify the obtained 
solutions. It turns out that the crucial parameters are: (1) the expansion coefficient B4 in 
the eigenfunctions decomposition of the nucleon distribution amplitude, (2) the form-factor 
ratio R = |GmI/^M5 (3) a criterion in the moment sum rules analysis that accounts 
for the higher stability of the lower moment sum rules relative to the higher ones [0. In 
addition, one can define a "hybridity" angle which quantifies the information about the 
mixing of geometrical characteristics of distribution amplitudes associated with different 
sorts of (asymmetric) longitudinal momentum balance. The upshot of this treatment is a 
pattern of unity in diversity. One is able to track the "metamorphosis" of the heterotic 
distribution amplitude across an orbit in the plane spanned by B4 and R as it transforms 
into the COZ one. Furthermore, the orbit turns out to be finite, starting at small R- 
values, associated with the heterotic distribution amplitude, and terminating at a COZ-like 
distribution amplitude which is a fixed point of the transformation. Interpolating solutions 
with comparable x^- values are also determined, showing various degrees of "hybridity" . In 
effect, the emerging orbit provides at a glance a global, and perhaps more perspicacious, 
perspective on the information encapsulated in available moment sum rules [35,57] for the 
structure of nucleon distribution amplitudes. 

Regardless of the degree of conviction of particular model distribution amplitudes, prob- 
lems are lurking in the endpoint region of phase space, where the typical gluon virtualities 
are much smaller than the external high momentum. Indeed, Isgur and Llewellyn-Smith 



15|, and also Radyushkin [16,17|, have pointed out that the contributions from this region 



in fact dominate the pion and nucleon form factors, rendering a perturbative treatment at 
accessible momentum transfers questionable. Their main conclusion is that the observed 
power-law behavior of the hadron form factors is not due to the perturbative contribution 
via hard-gluon exchange, but rather a reflection of the flnite size of the hadrons. The per- 
turbative contribution takes over, they argue, at momentum-transfer values too high to be 
measured in the foreseen future. Similar thoughts were more recently expressed by Bolz 
and Kroll on a phenomenological basis. These authors claim that the broad flat 

maximum of the (soft) overlap contributions to the pion and nucleon form factors in the 
intermediate momentum-transfer region "mimics" the Q^-dependence of the short-distance 
term. Again the conclusion is that the bulk of the existing data can be explained by the soft 
contribution alone without resort to the small perturbative contribution. Such statements 
could be construed as indicating that perturbative QCD is irrelevant for this sort of exclusive 
processes at laboratory Q'^. 

A serious criticism in this context concerns whether asymmetric distribution amplitudes, 
as those discussed above, are reliable. It was noted by Stefanis |^ that a finite set of low- 
order moments cannot provide a unique solution, and that the ensuing variation is inherent 
in the technique of moments inversion that is unable to render the diversity of possible 
solutions small. It is exactly because of this reason that a "hierarchical" criterion, which 
weighs the moment sum rules according to their order, was used in the "heterotic" approach, 
proposed in |]42| , ^ . Imposing such a global criterion, the sensitivity to disregarded higher- 
order terms is de facto averted [O]. 
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Another objection, intimately tied to the QCD sum-rule method, was raised by 



Radyushkin and expounded in several articles (see, e.g., [T^])- The point is that local con- 
densates assign zero virtuality to vacuum quarks corresponding, equivalently, to infinite 
correlation lengths of vacuum field fiuctuations and, since the higher-order moments are 
sensitive to the size of the non- locality (or inverse quark virtuality), this effect should be 
taken into account. This is perhaps a turning point in the conceptual evolution of QCD sum 
rules, but crucial technical challenges remain, e.g., the operator product expansion has to 
be generalized to nonlocal operators. Nevertheless, following this rationale, Mikhailov and 
Radyushkin and Bakulev and Radyushkin |^7[ (see also |^) were able to obtain pion 
wave functions, modeling the vacuum quark virtuality distribution by a Gaussian around 
an average virtuality ranging from 0.4 GeV^ to 1.2 GeV^ |^0|- ^ more recent analysis 



by Dorokhov, Esaibegyan, and Mikhailov |61|] attempts to calculate quark and gluon virtu- 
alities within the model of the liquid instanton vacuum []8|,|60[|. The wave functions derived 
this way are broader than the asymptotic solution of perturbative QCD, but, unlike the CZ 
amplitude, have no dip in the central region. It goes without saying that the corresponding 
second and fourth order moments are close to their asymptotic values. These investigations 
seem to confirm again the dominance of the soft contribution to the pion form factor at 
intermediate values of momentum transfer. 

Unfortunately, the extension of this method to the nucleon case is not straightforward. 
New nonlocal quark-gluon condensates of higher dimension emerge, which have to be es- 
timated, but cannot be reduced to condensates of lower dimension because they explicitly 
violate the "factorization hypothesis" [Q. 

Putting our comments so far together, we draw two important conclusions: (1) the end- 
point region of exclusive processes is IR sensitive and has to be considered very carefully; 
(2) the moment sum rules have defined a much more profitable goal to work towards than 
an elaborate formalism for calculating reliable distribution amplitudes. But however incom- 
pletely we may comprehend the details of how the complex structure of the non-perturbative 
QCD vacuum works, we all understand that until more advanced techniques are really avail- 
able, the moment sum rules and the herewith derived model distribution amplitudes are 
revealing. In succeeding sections, we focus on the ramifications and implications of issue (1) 
considering assessment (2) as a sound working hypothesis. 

While the non-perturbative non-locality of condensates may eliminate a source of errors 
in the determination of wave-function moments, another type of non-locality, namely in 
transverse configuration space, may help out in restoring the validity of the perturbative 
treatment in the endpoint region by permitting the divergences due to soft gluons to be 
screened by means of (re-summed) Sudakov effects. We do not propose to review these 
extensive investigations in great detail here, but we do sketch the key ideas in connection with 
form-factor calculations and complete this review in Sec. ^ with some concrete examples. 

In Abelian theories, like Quantum ElectroDynamics (QED), the emission of soft real 
(bremsstrahlung) quanta prevents the vanishing of the cross section by giving rise to BN 
cancellation between (soft) real and virtual photons [0. This cancellation does not occur 
within perturbation theory because these two different sorts of soft quanta correspond to 
different powers of the coupling constant, rent by the iteration procedure. Hence, one needs 
closed-form all-order expressions to accomplish cancellation. In the IR regime of non- Abelian 
theories (like QCD), where the ratios of all physical scales relative to the IR-cutoff scale are 
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large, cross sections of non-forward processes involving isolated colored particles, whether 
or not an infinite number of soft gauge bosons (gluons) are included, vanish in the IR limit. 
In contrast, processes involving only neutral (composite) particles have non- vanishing cross 
sections in this limit. It transpires that colored amplitudes at large momentum transfers are 
suppressed by damping exponential (Sudakov) factors - one for each non-color-singlet, near 
mass-shell particle []62| . 

The Sudakov form factor is the probability for no emission of soft photons (gluons) 
with increasing momentum transfer in hard photon-electron scattering within QED (QCD). 
Sudakov found |0 that the one-loop result exponentiates to give a double logarithm of the 



form exp — (Q;/c7r)ln (Q^/yU^) , fi being the small scale of the system; e.g., the invariant 
squared mass of the electron [^] (c = 2), or an auxiliary photon mass in the case of on-shell 
external fermions (c = 4) [Q,^,^ . The dominance of the leading double- logarithmic term 



in the high- momentum limit was confirmed by Mueller |63] , and Collins |^ by re-summing 
non-leading single logarithms to all orders of the fine-coupling constant a. The upshot is 
that, asymptotically, the Sudakov form factor drops to zero faster than any power of Q^. 

The generalization of this type of calculations to QCD was initiated by Cornwall and 
Tiktopoulos [^. They computed the non-Abelian vertex function up to the three-loop 
order and conjectured exponentiation of the leading double-logarithmic result. The formal 
proof of exponentiation was provided by Belokurov and Ussyukina [ |69[] , 
and Steiner 



70 



and by Dahmen 

The inclusion of single- logarithmic terms was supplied by Sen [Q, along 
the lines of Collins' previous QED analysis . An elaborate and systematic re-summation 
approach for leading and non-leading logarithms was conducted by Collins and Soper |72 



(see also |]75| , [71[ ), using eikonalization and renormalization group type equations in the axial 
gauge, with recourse to the Crammer- Yennie method ||75|. They found that the leading 



generic form exp < — cln In 



contributions stem from integration regions where either all four-momentum components 
of gluons are small ("soft region") or from gluons collinear to external (massless) quark 
lines.! When such loop momenta regions overlap they give rise to double logarithms of the 

/ o , r, in which the constant c depends on the color 

ln(MVA?jcD)J J 

algebra via the renormalization group effective coupling /3o = (H — 2nf/3) /4, where nf is 
the number of quark flavors. In our considerations to follow nf is taken to be equal to 3. 

One crucial clue of Sudakov effects in exclusive processes came into view, largely as 
the result of work by Botts and Sterman |7^. They realized that the IR cutoff in the 
Sudakov function can be traded for the inter-quark separation in transverse configuration 
space (after Fourier transformation), thus providing an in situ regularization of IR diver- 
gences without introducing external IR regulators. Conventionally, transverse momenta in 
the hard-scattering amplitude are dispensed with on the assumption that they are negligible 
compared to the large scale Q^. This is true - no doubt - as long as all momentum flows 
within a Feynman graph are large. However, in the endpoint region, as already mentioned, 
some momentum can formally come close to zero, becoming smaller than the neglected 
transverse momentum. Retaining the flow of transverse momenta into the hard-scattering 
region, is actually tantamount to modifying factorization because the momentum region at 



^Whether a gluon belongs to the soft or the collinear set depends on the gauge. 
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Figure 4. Flow chart illustrating two possible ways of realizing factorization of exclusive reactions in 
convolution form within QCD. 

the interface between the true confinement regime - parameterized in the wave functions - 
and the true hard regime - expressed through the short- distance part - is exphcitly taken 
into account in the convolution of the reaction amphtude. In the axial gauge, all Sudakov 
effects can be encapsulated in wave-function-like factors, one for each incoming and out- 
coming quark line, that link the soft wave functions with the hard scattering region. This 
amounts to a finite wave function renormalization (Stefanis The modified convolution 

(factorization) scheme is illustrated in Fig. |^. 

In this vein, Li and Sterman []7^ were able to show that in the pion case there is strong 
Sudakov suppression for large inter-quark separation as increases. Whence, contributions 
from the endpoint region, where perturbation theory fails, become less and less important. 
Indeed, the bulk of the pion form factor accumulates in regions where the coupling constant 
is less than about 0.7 and the transverse distances are moderate. This is a remarkable result, 
for no external regulators are needed to saturate the strong coupling at small momenta, the 
IR protection being intrinsically provided by the transverse inter-quark separation alone. 

As regards three-quark states, the situation is more complicated because several trans- 
verse scales are involved and therefore IR protection is not automatically accomplished. It 
was shown by KroU and Stefanis, and their respective collaborators , that the simple ex- 



tension of the pion analysis to the nucleon case by Li [7£] fails. The main reason is that there 



exist kinematical regimes where none of the Sudakov factors provides suppression for arbi- 
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trary choices of the inter-quark distances. One has to correlate the different transverse scales 
and adopt a common IR cutoff scale that is taken to be the maximum transverse separation 
("MAX" prescription). The underlying physical idea is that due to the color neutrality of a 
hadron, its quark distribution cannot be resolved by gluons with a wavelength much larger 
than an average inter-quark separation scale. Thus, gluons with wavelengths large compared 
to the (transverse) hadron size probe the hadron as a whole, i.e., in a color-singlet state and 
decouple. As a result, quarks in such configurations act coherently and therefore (soft) gluon 
radiation is dynamically inhibited. 

The particular advantage of this approach |[7S|P5| is that (1) it suffices to protect the 
amplitudes from becoming singular for all possible kinematical configurations, and (2) it 
yields a perturbative contribution to the nucleon form factors which saturates, i.e., which 
is rather insensitive to distances of order 1/Aqcd- Other IR regularization choices, for 
instance that adopted by Li, may lead to un-compensated singularities (see for more details, 
in Sect. |Vip. Results for both the proton as well as the neutron form factors were 
obtained, taking also into account the intrinsic transverse momenta in the nucleon wave 



function in the realm of the pion analysis of Jakob and Kroll The dark side of these 
theoretical improvements is that the calculated nucleon form factors fall short by at least a 
factor of two relative to the experimental data. The possible reasons for this depletion will 
be discussed in subsequent sections. 

The remaining part of this report is divided into two main parts. The first part deals 
with hadron distribution amplitudes (Sect. [ITD and form factors (Sect. |rv| ) within the 
standard convolution scheme. In the second part (Sect. [V| ) we treat similar topics within 
the modified convolution scheme. We close this review with a further in-depth discussion 
of the whole picture, attempting at providing a stimulating outlook for the future. Some 
important technical issues and useful formulae are collected in the appendixes. 

III. STANDARD CONVOLUTION SCHEME 

The theoretical tools for the description of exclusive processes are the hard-scattering 
amplitude - which describes the process-dependent quark-gluon interactions within pertur- 
bative QCD - and the probabihty amplitude for finding the lowest-twist Fock state {alias 
the distribution amplitude for the valence state) in a hadron. The total exclusive amplitude 



is then represented by the convolution of these three factors |21,23|, assuming factorization 
of large-momentum fiow regions from those of soft transverse momenta, necessary to form 
bound states. In terms of the operator product expansion, this corresponds to truncation 
at leading twist (i:=dimension-spin) level, namely t = 2 (meson case) and t = 3 (nucleon 
case). Higher-twist components, corresponding to a higher number of partons (quark-pairs 
and gluons), are suppressed by powers of the large scale, i.e., the momentum transfer Q^. 
To clarify the physical meaning of these statements, we now turn to factorization. 

A. FACTORIZATION 

Factorization theorems are of central importance in quantum field theory. The basic idea 
is that one can separate high-momentum from low-momentum dependence in a multiplica- 
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tive way. For example, proving that UV divergences occurring in Feynman graphs can be 
absorbed into multiphcative renormahzation factors (infinite constants) is instrumental in 
establishing renormalizability. The technical challenge is to prove factorization of a partic- 
ular QCD process to all orders in the coupling constant going beyond leading logarithms 
pj| . These difficulties derive from the fact that in QCD a new type of IR-divergence is 
encountered, the collinear divergence, and that in higher orders the self-coupling of gluons 
becomes important in the exponentiation of IR divergences. 

The realization of factorization when applying to elastic form factors can be written in 
the form of a convolution of the hard-scattering amplitude (dubbed Th) and two soft wave 
functions corresponding to the incoming and out-coming hadron (termed $). Generically 
(i.e., absorbing all integrations over internal variables into ®), one has 

F(Q2) = ^ TMQ) ® <f '"(m//i) , (1) 

where m sets the typical (small) virtuality in the soft parts and Q is the (external) large 
scale, characteristic of the hard (parton) subprocesses. 

The matching scale fi at which factorization has been performed is arbitrary and, as- 
suming that ^ ^ m, it can be safely identified with the renormahzation scale - unavoidable 
in any perturbative calculation - by virtue of the renormahzation group equations. In this 
way, F can be rewritten as a function only of the coupling constant operative at that same 
scale, the latter being identified with the large external scale. 

As long as scale locality is preserved, i.e., as long as the variation of the effective coupling 
constant with yu is governed by the same momentum scale, and the limit m — > is finite, 
Eq. (0) is valid because intrusions from the hard into the soft regime are prohibited. Potential 
IR divergences in Th are removed by subtractions on account of the properties of the wave 
function parts. This means that Th in leading order is by definition insensitive to long- 
distance interactions, i.e., it is IR safe. All IR-sensitivity resides in the hadron distribution 
amplitudes <|)^°(°"*) which are peaked around small transverse momenta of the order with 
their large-momentum tails removed by cutting off from above the integration over transverse 
momenta. Their dependence on the cutoff scale is mild and is governed by renormahzation 
group-evolution. Both the subtraction procedure of the large momentum tails in the soft 
parts and the cancellation of IR divergences in the hard part are not unique. They actually 
define the factorization procedure adopted, thus imposing an implicit factorization-scheme 
dependence on the hadron distribution amplitudes [Q. Of course the asymptotic behavior 
of F should not be affected and still be determined by the asymptotic distribution amplitude 
evolved with the leading (i.e., lowest-order) anomalous dimension, associated with vertex 
and quark self-energy corrections. 

B. SHORT DISTANCE PART 

Invoking factorization, the leading order expression for the helicity-conserving hadron 
form factor can be cast in the form 

T(g2) = f\dx] f\dy]<^>*{y„Qy)Tn{x„y„Q)^x„Q,) , (2) 
JO Jo 

where 
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Figure 5. Example of a Feynman graph contributing to the nucleon form factor at tree level. 



[dx] = 6 1 — X! n '^^i (^) 

V i=i / j=i 

{n = 2 for mesons and n = 3 for baryons), and Qx = niin(xjQ) or Qx = niin ((1 — Xi)Q) 
with analogous definitions for the y variable. The form factor is the probability for the 
hadron to absorb large transverse momentum while remaining intact. The bound-state 
(i.e., confining) dynamics is encoded in $, while in the hard-scattering amplitude Th the 
hadron is simulated by on-mass-shell (but off-light- cone energy) valence quarks with neg- 
ligible mass and transverse momentum. To leading order ag, Th is the sum of all Born 
diagrams contributing to the particular process. The transition from the initial to the final 
state is supposed to go via hard gluon re-scattering which involves a factor (as(<5^)/Q^), 
[Q"^ = q]_ = —q^) after one quark was derailed from the initial to the final direction. To 
leading order, the contributions of gg-irreducible diagrams to the pion form factor amount 
to the hard-scattering amplitude 

T(.,,.Q) = 15^11^. (4) 

where the abridged notation x = 1 — x,y = 1 — y has been used, C-p = {N^ — \)/2Nc = 
4/3 is the Casimir operator of the fundamental representation of SU{3)c, and a^iQ"^) = 
47r//3o ln((5^/AQ(^j-,)j is the running coupling constant in the one-loop approximation. 
In the nucleon case, there is a total of 14 Born diagrams with different topologies out of 
which only 8 give non-zero contributions. A complete list of these diagrams is compiled in 



31 1; the Feynman rules for light-cone perturbation theory are given in [21|. A typical dia- 
gram contributing to the nucleon form factor is depicted in Fig. |^. The quark-propagator de- 
nominators along the first and the second quark lines are, respectively: —Q'^xi and —Q'^x^yi^ 
those of the gluon propagators which channel the external momentum flow from the struck 
quark to the spectator quarks are: —Q^Xiyi and —Q^x^y^, where all quark masses and 
/cj^-momenta have been set equal to zero (coUinear approximation). Recall that in the light- 
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cone frame, parton "i" in a hadron has four- momentum pi = XiP + ki = {p'^,p ,p±) , where 
Xi = ^ and ki = (0, k^). Momentum conservation imphes 



n 



^x, = l; ^k^, = 0. (5) 



i=l i=l 
n2 _ ^2 



Quarks are on- mass shell, i.e., p^ = mi but off the light-cone energy: 

■ ''^ 

In a frame where P = (P+ = 1, 0, 0^) and for massless quarks, this relation simplifies to 

Note that the sum over all p~ is not equal P . The difference is a boost-invariant measure 
of how far off energy shell a Fock state is [^,^]. This off-shellness is large in the kinematic 
endpoint region, i.e., when k^^j or is small and, as a consequence, the hadron wave function 
should vanish in these limits. Hence, formally, all wave functions should satisfy the boundary 
conditions 

k^j ipn{xi, k_Li, Ai) ^0 as k^^ ^ oo 

(8) 

ip„{xi, k_Li, Xi) as Xi ^0 

if the free-particle Hamiltonian is to have a finite expectation value (Aj denotes here parton's 
"i" helicity). None of these conditions is generally satisfied in the absence of UV (k_L —>■ oo) 
and IR {x — > 0) regulators (for a further discussion of these subtle issues, we refer to |]83[|). 



C. LARGE DISTANCE PART: PERTURBATIVE ASPECTS 

1. Renormalization 

The function ^(xj, Q) is the probability distribution amplitude for finding the hadron to 
consist of valence quarks each carrying a fraction < Xj < 1 of the hadron's longitudinal 
momentum P+ (P^ = (P° ± P^) / a/2) at transverse separations (relative to the hadron's 
P+ momentum) not smaller than Q~^. We work in the infinite momentum frame, where 
P = (0, 0, P^) and P^ = |P| —>■ oo. In this frame the p^ component is large and conserved 
(the same for p), and the p~ component is small and not conserved (the same for E). 
Furthermore, boosts along the -|-z-direction ("-(-3" -direction) leave the transverse momenta 
unchanged. In leading order, Qx{y) can be replaced by Q, which is then the only large 
scale involved. The small scale of the system, called m in Eq. (|l]), does not enter the game 
explicitly (provided it is much smaller than Q); it nevertheless plays a crucial role in the 
determination of $ and is the non-perturbative counterpart of the perturbative scale Aqcd- 
Physically, the existence of such a scale indicates that gluons with wavelengths larger than 
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m ^ do not "see" individual quarks and decouple. One should eventually be able to relate 



m with typical non-perturbative scales, like the average quark virtuality in the vacuum pQ 
mentioned in the Introduction, or the size of instantons P,|60[1. 

Dispensing with the fc^-dependence in Th, allows to integrate out this variable in the 
wave function parts up to the factorization scale /^^ and replace the (renormalized) wave 
function by the distribution amplitude 

/ 2 \ -C7f//3o n 

V ^^QCD/ i=l 

where 

Ek..]n^, (10) 

and 7f is the anomalous dimension associated with quark self-energy in the light-cone gauge: 

,, = C,(l + 4/'d.^). (11) 

The logarithm in front of the integral in Eq. stems from UV divergences due to gluon 
radiative corrections in the hard-scattering amplitude. Hence, all vertices and propagators 
have to be renormalized. However, Zf, the gluon renormalization constant, and Zi^, the 
quark-gluon renormalization constant, can be absorbed (removed) by replacing as by its 
renormalized value al'^'^ = {Z2 / Zi^)"^ Z§ as so that the only renormalization constant to be 
determined is Z2 which renormalizes the quark propagator (or, by taking its square root, 
each incoming and outgoing quark line). But Z2 can be computed from the gluon radiative 
corrections to the quark-photon vertex by virtue of the QED Ward identity Z^"' = . It 
is evident from our factorization prescription that only momenta Q < Q are included in 
the wave function, meaning in turn that the virtual photon's momentum probing it is also 
limited to Q. Hence, in computing the UV-divergent part of the photon-quark vertex, the 
photon momentum can be ignored, keeping only gluon loop momenta greater than Q. Then 
Z2 — > Z2{Q) and as Q increases, Z2{Q) (i.e., the probability to find a bare quark) decreases. 
The asymptotic result is 

Z^ = hm In , (12) 

where we have used the fact that in leading order k± = minjxjQ}, (or k± <^ Q^^ = 

mm{xiQ}) so that each individual "tilded" wave-function scales can be traded for a common 
factorization scale /i^. This renormalization procedure ||8^ yields c = 1 for mesons (two 
quark lines) and c = 3/2 for baryons (three quark lines). 

The gauge-invariant distribution amplitude ^^^\xi, fj?) is intrinsically non-perturbative 
and - provided the same factorization scheme is used to cast exclusive amplitudes for differ- 
ent processes in convolution form |^ - universal. The large- momentum behavior of these 
functions can be analyzed either using operator product expansion techniques or, equiv- 
alently, by evolution equations analogous to DGLAP equations [pllKsLKGl in deep-inelastic 
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scattering. Following the second approach, one takes derivatives with respect to of Eq. 
to arrive at evolution equations of the generic form [21 



(91nQ2 



[dy]V {x,,yi,a,{Q^))^'^'^\x,,Q 



(13) 



with distinct kernels K (xj, j/j, q;s(Q^)) for each process at hand, which, to leading order in 
tts, are computable from the single-gluon-exchange kernel. 

To solve the evolution equation, ^^^^ for hadron H has to be expressed as an orthogonal 
expansion in terms of appropriate functions which constitute an eigenfunction basis of the 
particular gluon-exchange kernel, i.e.. 



n=0 



r^7F(^(/i )) 



(14) 



where ^^^-^ is the renormalization group asymptotic distribution amplitude (see below) be- 
ing proportional to the weight w{xi) of the particular orthogonal basis, and denotes 
the corresponding eigenfunctions. The coefficients -B^^'' of this expansion are associated 
with matrix elements of composite lowest-twist operators with definite anomalous dimen- 
sions (after diagonalization of the evolution kernel) taken between the vacuum and the 
external hadron. They represent the non-perturbative input (integration constants of the 
renormalization group equation) in Eq. (|l^ and have to be determined at some initial scale 
of evolution /i^ by non-perturbative techniques. The exponential factor in Eq. (|l^ takes 
care of momentum evolution according to the renormalization group and is governed by the 
quark anomalous dimension 



dZ^ da. 



7f 



_/f 

da. 



(15) 



which in the axial gauge is [87 



7g 



-- + 0(a 

TT 



rmsJ 



(16) 



2. Meson evolution equation 



The advantage of using an eigenfunctions decomposition is that the evolution equation 
can be solved by diagonalization. Consider, for example, the meson evolution equation. 



Using the evolution "time" parameter p5|j21 



as{k 



In 



asifi^ 



In 



and the relation 



d 



as(Q2) "^In/iVA^jcD 



(17) 



Q2 



(18) 
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the meson evolution equation reads 



X1X2 



d_ ; 



Po 



(3q Jo 



V{xi,yi) ^ixi,Q^ 



(19) 



where = xiX2^{xi,Q'^). Expanding in terms of eigenfunctions $„,, or xiX2^n, one 

has 



9 2 2 

from which one readily obtains 

^n(^Xi^ Q ) ^n(^Xi^ ^ ) e 



(20) 



<l>„(xj,/i^) exp 



-7n In 



tts(/i^ 

«s(Q2) 



XiX2$„(xi,/i^) In 



-7n 



A2 



(21) 



Then the equation to be solved becomes 



a^ia:2 -5 In] ^n{xi) = -r [dy] V{xi,yi) $„(?/») 



/5o 



/3o -'0 



(22) 



where V{xi, Ui) is both real and symmetric, and {$„(xi)} form a set of orthogonal functions 
with respect to the weight w{xi) = X1X2: 



[dx] w{Xi) ^„{Xi) $m(Xi) = Kn 6nr 



(23) 



Within this basis of eigenfunctions, the pion distribution amplitude has a convergent expan- 
sion for all Q^, viz.: 



<l>^-\x„Q^)=w{x,) gi^W In 



Q 



2 \ 



n=0 



A2 



(24) 



with expansion coefficients 



In 



-In 



A2 



/ [t/x]«;(x,)l>W(x,)l>W(x„g2) . 



(25) 



To determine the anomalous dimensions 7„ and the corresponding eigenfunctions, it is 
convenient to introduce the relative coordinate — l<C = ^i~2;2<l and express the kernel 
V in terms of a matrix Ujn via a monomial basis 



nc)io = ^(c)ElOf^." 



(26) 



where w{Q = X1X2 = t(1 — C^)- The matrix Ujn turns out to be triangular ||21 
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j > n 



(27) 



25, 



h\h2 



(n + l)(n + 2) 



n+l 1 

-4E 



fc=2 



k 



J = n 



(28) 



-1) 
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hih2 



+ 



4(j + 1) 



(n + l)(n + 2) (n-j)(n + l) 



j <n . 



(29) 



Hence, the eigenvalues are just f/„„ and the eigenf unctions are polynomials of finite order. 
The only polynomials in the interval [—1, 1] orthogonal with respect to w{Q = ^(1 — C^) are 
the Gegenbauer polynomial^ C^/^(C) (see, e.g., ^^), normalized by 



where 



Then, on account of 



{2 + n){l+n) 
2(2n + 3) 



the associated anomalous dimensions become 



In 



9i 

/3o 



hih2 



1 2S 

(n + l)(n + 2) 



> {n even) 



(30) 



(31) 



(32) 



(33) 



where, for the pion, 6^^^^ = 1 (see Fig. ^. The corresponding eigenfunctions are [^>:^''[Xij 



C3/2(xi) = C3/2(xi - xa)) SO that 



^^^\xi,Q^) = xiX2Cl/^ixi - X2) In 



A2 



(34) 



with coefficients going like (cf. Eq. (|25|) ) 



In 



-7n 



A2 



2(2n + 3) 

{2 + n){l + n) J-i 

y/2{2n + 3) 1 



{2 + n){l + n) ViV, 



(01^(0)7+75^^/2 Ud+) ^(0)|7r)(Q.) , 



(35) 



^The Gegenbauer polynomials correspond to the conformally invariant operator product expansion 
for two spin-i operators |23,^|^J^ 
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where in the last step we switched from the relative variable ( back to x, and the matrix 
element of the local operator is evaluated with UV cutoff Q^. Here = D — D with 
D = dfj, — igY^^^it^A"^ and D = + igY^f^it^A^^, and = in the light-cone gauge. 
From this equation one sees that the expansion coefficients B^^'^ are matrix elements of local 
operators and decrease like provided Q^) < Kxj as — ^ for some e > pi 

Thus to leading order the pion distribution amplitude becomes {P^ = 1) 



V2{2n + 3) 



■Cn^x - 1] 



^o(2 + n)(l + n) v^, 



(36) 

For asymptotic values of only the leading logarithm in Eq. (|35| ) with the least anomalous 
dimension 7o = survives (Cq^^ = 1), so that 



f^x{l-x) 



(37) 



where /^r = 93 MeV is the pion weak decay constant. From this, one infers the normalization 
condition (sum rule) 



2VN. 



(38) 



which, given the shape of ^^'^\xi, Q^), normalizes it for every value of Q^. 

The above treatment can be readily used to provide information on the asymptotic 
behavior of the pion form factor as well. The electromagnetic pion form factor in the space- 
like region is defined by 



(39) 



where J^™ = X]f ^i<iilij.<ii is the electromagnetic current, and the index f denotes different 
flavors. To relate the pion form factor to the Fock state wave functions, discussed above, 
the /i = + component of this equation is taken. 

With the spectrum of eigenfunctions and corresponding anomalous dimensions in hand, 
we can write Eq. (|^) in the form {n even because of C-parity) 



167rCFas(Q^ 



ei 



A2 



dxi Xi C^'^{xi — X2) 



+ (1^2) 



(40) 



Recalling that for tt"*' the charges of the struck quarks sum to one, and that for even Gegen- 



bauer polynomials, \q dx\ x\ C'^'^{x\ — X2) 
form factor becomes in even) 



1 fi 
2 



Jo dxx (xi X2) Cll'^{xx - X2) 



2' 



the pion 



47rCFas(g' 



-In 



A2 



(41) 
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Figure 6. Anomalous dimensions of thrcc-quark operators up to polynomial order A-I = 9 in comparison 
with those of two-quark (mcson-liko) operators. Curves are shown for antisymmetric {Sn = —1) and 
symmetric (Sn = 1) eigenfunctions of the diagonalized nucleon evolution equation. The solid line is an 
empirical fit which is compatible with power-low behavior. 



23 



In the limit of asymptotic values of momentum transfer, only the zeroth order term Cq^"^ 
contributes with anomalous dimension 70 = 0, so that 



47rCFas(Q 



2^ 



hm F.(Q^) = --^-l^^ ' . (42) 



This result can be couched in the final form 



Jm FM) = 15%<^ (43) 

by virtue of the decay process 7r+ — > l^vi, il = /i^,e+) which fixes Bq""^ = in terms 
of the pion decay constant /yr, independent of the momentum variance. Unfortunately, the 



asymptotic prediction is not supported by the existing data Indeed, evaluating the 

above expression, say, at = 6 GeV^ for Aqcd = 0.2 GeV, one finds Q^F^(Q^) = 0.12 
which is about a factor of three below the experimental value (and well below the error bars). 
The conclusion is that additional contributions of Gegenbauer terms in the expansion of the 
pion distribution amplitude have to be included, or that one has to take into account higher 
twist contributions (power corrections). In order to increase the value of the pion form 
factor at moderate Q^-values, a broader pion distribution amplitude, relative to X1X2 = 
j{l — x^), is needed with at least i?2'^'*/-^o'^^ > 0. This goal can be achieved using QCD 



sum rules |29,56,93]. The literature on this subject is vast and the selection of references 



, p5|P^ , P7|P^ , P5|jl 00|| is by no means complete. We apologize to those authors whose works 



have not been included. 



3. Nucleon evolution equation 



Let us now turn to the nucleon case. The evolution equation is 



X1X2X3 



^$(x.,g^) + |^$(x.,Q^) 



^ j\dy]V{x,,y,)~^y,,Q') 
Po -Jo 



(44) 



where $ = X1X2X3 $, Cp is defined below Eq. (^), and 



V{Xi, Vi) = 2X1X2X3 '^iVi - ^i) ^iVk 



Xj 



'hihj 



+ 



Xi -\- Xj Ui X{ 



(45) 



Note that V{xi,yi) = V{yi,Xi) is the sum over single-gluon interactions between quark 
pairs {i,j}, and the subtraction prescription A$(?/j, Q^) = ^{yi, Q"^) — ^{xi, Q^) ensures IR 
finiteness at Xi = yi, i.e., V{xi, yi) is not a function but a distribution. For antiparallel spins 
^hjij = 1 for parallel spins it equals 0. Note also that if no gluon is exchanged, each 
Xk in the initial and final wave function is the same because no longitudinal momentum is 
introduced by {q^ = 0). 

We propose to solve the evolution equation by employing factorization of the depen- 
dence on longitudinal momentum from that on the external (large) momentum scale (cf. 
Eq. (p!4|)). The latter is renormalization group controlled according to 
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d ~ 



with solutions 



In 



A2 



(46) 



(47) 



This allows us to write the full nucleon distribution amplitude in the form 



Q"^) ~ X1X2X3 ^n{Xi) In 



n=0 



A2 



(48) 



where ^n{xi) are appropriate but not tabulated polynomials, and the expansion coefficients 
Bn encode the non-perturbative input of the bound-states dynamics at the factorization 
(renormalization) scale. Their determination will concern us in the next section. 

From the factorized form of ^n{xi, Q^) in Eq. (^71) , it follows that the evolution equation 
for the x-dependence reduces to the characteristic equation 



X1X2X3 



3CV 
2J0 



-In 



Hxi) 



(3q Jo w{xi) 



(49) 



where w{xi) = X1X2X3 = —xi — X3)x3 is the weight function of the orthogonal basis and 
Cb = {Nc + l)/2Nc = 2/3 the Casimir operator of the adjoint representation of SU{3)c. To 
proceed, it is convenient to conceive of the kernel V{xi, Ui) as being an operator expanded 
over the polynomial basis I^T] 0:3) = |fc /), (recall that because of momentum conservation, 
only two out of three Xi variables are linearly independent) i.e., to write 



V 



V{xi,yi) 



and convert Eq. ( ^91) into the algebraic equation 

^n{Xi) = Jn^n{Xi 



3CV ^Cb V 
2 (3o (3o 2w{xi 



(50) 



(51) 



In this way, the action of the operator V can be completely determined by a matrix, namely: 



V\kl > 
2w{xi) 



I i+j<M 



(52) 



The corresponding eigenvalues are then determined by the roots rjn of the characteristic 
polynomial that diagonalizes the matrix U: 



V^n{Xi) = -r]nW{Xi)^n{Xi) , 

SO that the anomalous dimensions for order M are given by 

1 



7„(M) = - ( -Cf + 2r]n{M)CB 



(53) 



(54) 
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where the orthogonahzation prescription 

1 . . 1 

{dx]w (Xi) ^,riiXi)^niXi) = — 5mn (55) 
I\m 

has been employed with A^^ being appropriate normahzation constants (see, Table p. 
The explicit form of the matrix U was derived by Lepage and Brodsky |^ and is repro- 



duced in Appendix A. Within the basis \kl), the matrix U can be diagonalized to provide 
eigenfunctions, which are polynomials of degree M = + / = 0, 1, 2, 3 . . ., with M + 1 eigen- 
functions for each M. This was done in |^ by diagonalizing the (M + 1) x (M + 1) matrix 
Uij^ki with i+j = k + l = M and results up to M = 2 were obtained. In our approach, re- 
ported in several meetings P5| , |101| , pi?^ and worked out in [ |103| , |104| ], we make use of another 
method that is based on symmetrized Appell polynomials |91| . 

Appell polynomials are special hypergeometric functions (see Appendix A) of the form 

:^^^^(5,2,2;a;i,X3) = J>nn(a;i, Xs) , (56) 

which constitute an orthogonal polynomial set on the triangle T = T(xi,X3) with Xi > 0, 
X3 > 0, Xi + X3 < 1. They provide a suitable basis for solving the eigenvalue equations for 
the nucleon because within this basis V is block diagonal for different polynomial orders. 
Moreover, introducing a "symmetrized" basis of such polynomials according to [|103|| 



mn('^l5 -^s) 2 mni^Xi, 2^3) i ^ nmi^Xi, 2^3)] 
k-\-l<m+n 

= E Zr\kl) (57) 

k,l=0 

(where + refers to m > n and — to m < n), V commutes with the permutation operator 
P13 = [321] and thus becomes block diagonal within each sector of (definite) permutation- 
symmetry class of eigenfunctions for fixed order M. As a result, the kernel V can be 
analytically diagonalized up to order seven. This is related to the fact that the characteristic 
polynomial of matrices with rank four can be solved analytically. Beyond that order, its 
roots have to be determined numerically. We present here results up to order M = 4. 
Still higher-order eigenfunctions (up to M = 9) were obtained by Bergmann |104|. The 



differential equation defining the symmetrized Appell polynomials is given in Appendix A. 

To appreciate the usefulness of this type of approach, some historical remarks ought to 
be made before we proceed. The solution of the eigenvalue equations for the nucleon beyond 
leading order is a long-standing problem. As already mentioned, results up to order two 



were obtained by Lepage and Brodsky |21]. The anomalous dimensions they computed were 
subsequently confirmed by Peskin ||105|| who considered composite three-quark operators 
containing derivatives and having baryon quantum numbers. Such operators interpolate 
between the nucleon (or the A resonance) and the vacuum at leading twist-three. Their 
anomalous dimension was extracted from the divergent parts of their matrix elements by 
diagonalizing the renormalization matrix of these operators (see in this context also ||106|| ). 
Then Tesima [ p. 07)1 analyzed the light-cone behavior of the Bethe-Salpeter wave functions 
of three-quark bound states with the use of the operator product expansion in terms of a 
conformally invariant operator basis. He presented results for twist-three operators up to 
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Table 1. Orthogonal eigcnfunctions ^ni^i 



kl '^kl ^1^3 



of the nucleoli evolution equation 



up to polynomial order M = 4 in terms of the coefficient matrix a^^ (a^^ = Sn a^p. - no summation 
over n implied) and the corresponding anomalous dimensions 7„ defined in the text. The numerical 
results for n>12 have been obtained with a much higher numerical accuracy than shown in this 
table. 
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order two which in general coincide with Peskin's results and also with those obtained by 
Lepage and Brodsky. In addition, he obtained results for higher orders, but published only 
the eigenvalues and eigenfunctions corresponding to order three. However, his higher-order 



results are not supported by our calculation. Moreover, it was shown by Ohrndorf ||108|| that 
coUinear conformal covariance alone does not suffice to fix trilinear quark operators with 
derivatives uniquely. 

We shall now give a brief description of our approach and present the main results. All 
nucleon eigenfunctions can be represented as linear combinations of (symmetrized) Appell 
polynomials of the same order M: 

m+n=M 

^kix,)= 4n^^n(5,2,2;xi,X3) (58) 

m,n=0 

since Appell polynomials of the same order are not orthogonal to each other. For instance, 
for M = 1, one finds 

l>i(x,) = ^-^1^(5, 2, 2; xi, X3) - ^ff (5, 2, 2; Xi, X3) 

= xi - X3 (59) 

and 

^2ix^) = J-ff (5, 2, 2; xi, X3) + ^iP(5, 2, 2; Xi, X3) 

= -2 + 3(xi + X3) , (60) 

where the notations and conventions of are adopted.! Before diagonalizing, it is conve- 
nient to rearrange J-'mn, which belongs to a definite symmetry class 5'„ = ±1 within order 
M, in the form of an (arbitrary) vector: 

^mn{Xi,X3) I > Tq{Xi,X^) . (61) 

Then Hilbert-Schmidt orthogonalization yields a basis 

\^',) = i:Zli\kl) (62) 

k,l 



with 



I [dx] w{xi) P^P^, oc 5qq' , (63) 

•J 



so that 

V\P,) 1 



2E/a«.M<lu>. (64) 



In particular, the coefficient B2 has the opposite sign relative to |21| 
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Note that the construction of polynomials depending on two variables via the Hilbert- 
Schmidt method has no unique solution, but depends on the order in which the orthog- 
onalization is performed. Since beyond order M = 3 neither the eigenvalues nor the normal- 
ization factors are rational numbers, one has to find which representation is more convenient 
for calculations. 

The last step in determining the eigenf unctions and eigenvalues of V is to define the 
matrix 

Mq'q= f [dx]w{xi,{l- Xi- X3),X3):F'g,{xi,X3)V{xi,X3)yg{xi,X3) (65) 

and calculate the roots of the characteristic polynomial 

r{r]) = det [Mq>q - 7]Ig'g] . (66) 

Consequently, in terms of the eigenvectors = (m^, . . . m^'^ of M.q'q, the eigenf unctions 
of the evolution equation are given by 

4>g(xi,X3) oc ^m^'j^g,(xi,X3) 

= T.4i\kl)- (67) 

k,l 

For every order M, there are M + 1 eigenfunctions of the same order with an excess of 
symmetric terms by one for even orders. The total number of eigenfunctions up to order M 
is '^max(^) = ^{M + 1)(M + 2) and the corresponding (M + 1) eigenvalues are obtained 
by diagonalizing the (M + 1) x (M + 1) matrix. A compendium of the results up to order 
M = 4, meaning a total of 15 eigenfunctions and associated anomalous dimensions, is given 
in Table |I[ The precision of orthogonality is at least 10~^. To get acquaintance with the use 
of Table [l|, we write out explicitly one of the eigenfunctions contributing to order M = 3: 

$9 = floo + '^10 (^1 ~ + ClnXiXs + (^1 ~ ^3) + dli^iXs {xi - X3) + (xf - xfj 

( 2 2\ 59 + V4801 , , 601 -74801 / o s\ / n 
= (xi - X3) - 3 [x^ - XgJ + xxX3 (xi - X3) + - X3J . (68) 

It turns out that the eigenfunctions {$fc} of the nucleon evolution equation satisfy a 
commutative algebra subject to the triangular condition \0{}z) — 0{V)\ < 0{ra) < 0{k) + 
0{l): 



^k{xi)^i{xi) = E Fj:i^^{xi) (69) 

m=0 

with structure coefficients given by 

Ffc7 = Nm [ [dx] XiX3il - Xi - X3)^m{Xi)^k{Xi)^liXi) , (70) 

Jo 

0{k) being defined by 
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0{k) 



r 

1 

2 
3 
4 



= 
l<k<2 
3<k<5 
Q<k<9 
k = 10,11 



(71) 



The structure coefficients are symmetric, i.e., FJ^ 



. Furthermore. 



^kk 



The utihty 



,1 - FU:. Furthermore, F^, = f . 
of this algebra derives from the fact that once the structure coefficients have been computed, 
they can be used to express any function f{xi,X3) in terms of the nucleon eigenfunctions. 
For example, one can calculate the Husini function 



h 



[dx] 



(72) 



Ik 



which gives the probability for finding ^'^^^ in a particular solution The values of F^ 
up to 0{k) = 11 are tabulated in ||104| . 

The obtained results for the exponents (anomalous dimensions) governing the scaling 
behavior of the nucleon distribution amplitude (cf. Eqs. (^Tf ), (48D) are shown in a series 
of four figures associated with successively increasing order MB As outlined above, the 
eigenfunctions correspond to trilinear quark operators with definite anomalous dimension. 
Such operators can be renormalized in a multiplicative way. But as the order M increases, 
so increases also the number of derivatives in these operators entailing a strong amount of 
mixing owing to the fact that they carry the same quantum numbers. As a consequence, 
the anomalous dimensions are degenerate within fixed symmetry classes and hence a "mul- 
tiplet" structure emerges. Because all 7„ are positive fractional numbers increasing with n 
(i.e., order M), higher terms in the eigenfunctions decomposition are gradually suppressed. 
Fig. H displays the anomalous dimensions 7„ up to order M = 9, distinguishing between 
symmetric and antisymmetric eigenfunctions. The trend line of this pattern seems to follow 
the empirical law (solid line) 7„ = 0.37(9(n)°'^^^. The illustration of the spectrum up to 
order 20 is displayed in Fig. 0. 

The next two figures. Fig. ^ and Fig. |^, address the large-order behavior of the anomalous 
dimensions. Also here, both symmetry classes under the permutation P13 are shown in 
unison: open circles for values belonging to S'„ = 1 and black dots for those belonging to 
5'n = — 1. As the order of eigenfunctions increases, a different picture for the large-order 
behavior of the spectrum of anomalous dimensions develops, namely, one of logarithmic 
increase. Indeed, in the limit of very high polynomial orders, symmetric eigenfunctions tend 
to degenerate with their antisymmetric partners as one observes from Figs. |^ and ^, where, 
with increasing order, the dots tend to enter the centers of the circles. At the same time, 
the eigenvalues 7^ increase logarithmically. Hence, to order 20, the spectrum of anomalous 
dimensions as a whole is better described by a logarithmic ffi (solid line) of the form 7^ = 
1.25[logio (C(n) + 1.37)]^-3^ The dotted curve shows the power-law ffi 7„ = 0.52O{nf-^^^ 
which is no more supported by the higher-order values. A more accurate prediction for the 



^This part was done in collaboration with Michael Bergmann. 
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Figure 7. Pattern of anomalous dimensions of three-quark operators up to M = 20. Values associ- 
ated with symmetric eigenfunctions are denoted by open circles. Those corresponding to antisymmetric 
eigenf unctions are marked by black dots. 
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Figure 8. Pattern of anomalous dimensions of three-quark operators up to M = 37. Values associated 
with symmetric eigenfunctions are denoted by open circles. Those associated with antisymmetric eigen- 
functions are marked by black dots. The solid line is an empirical fit already compatible with logarithmic 
behavior. 



asymptotic behavior of 7„ is provided by Fig. ^ where the upper envelope of the spectrum 
is well-described by 7„ = [log]^Q(2.13(9(r2) + 1.4)]^'^^. These findings, presented here for 
the first time, confirm expectations based on exponentiation assumptions of leading-loop 
contributions [|6^ , |109| , p!T0|Jlll|| . Physically, the logarithmic rise is due to enhanced emission 
of soft gluons, reflecting the fact that the probability for finding bare quarks decreases. 
However, in order to establish this trend, still higher orders have to be computed, at least 
up to M = 10^ (three points on a logarithmic plot). 



D. LARGE DISTANCE PART: NON-PERTURBATIVE ASPECTS 



Given the eigenfunctions and associated anomalous dimensions of the nucleon evolution 
equation, the only unknown quantities entering the description of the nucleon distribution 
amplitude are the expansion coefficients i.e., the set of matrix elements of tri-local 
operators interpolating between the one-nucleon state and the vacuum. The QCD sum 
rules method is presently the best approach for extracting these quantities, albeit there are 
conceptual limitations. We do not discuss the pion here since it has been studied extensively 
elsewhere p7|,PMP2| . 



Specifically, we consider the following gauge-invariant proton wave function 

(0| [C(^i, zM) Ua{z^)T [C{Z2^ zM) ci^(z3)|P)^ e'^'^ (73) 



where a, 6, c and 7 are color and spinor labels, respectively, A^(x) = Z]i=i 

the Lie-algebra valued gluon fields {f being the generators of SU{?>)c)- The operator 
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_i I 1 I 1 ! I I ■ I ■ I • t I 

20 40 60 80 100 120 140 

polynomial order 

Figure 9. Pattern of anomalous dimensions of three-quark operators up to M = 150. For technical 
reasons, antisymmetric eigenf unctions are included only up to order M = 37. This has no substantial 
influence on the large-order structure. The obstruction at order M = 45 is caused by a change in the 
numerical precision. The solid line is a logarithmic fit to the upper envelope of the spectrum up to order 
M = 90. 
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C{zi,Zj\A) = V exp 



(74) 



denotes a path-dependent phase factor ("connector" ||113|| ) in which the expansion of the 
exponential along a contour C joining the points Zi and Zj is controlled by the path-ordering 
prescription V. For a short-distance expansion, the connector can be evaluated along straight 



lines with lengths oc 



\Z±i — Z±j 



1/2 



as Q — > cxD. Working in the light 



cone frame, one can use the axial gauge A'^ = by virtue of which the connector at short 
distances can be replaced by unity since all inter-quark separations along the light cone are 
short. A covariant gauge in which the connector can also be neglected is provided by using 
the gauge parameter a = —3 ||113| , |114|| . By virtue of the light-cone gauge, the contour factors 
reduce to unity and the path-dependent matrix element simplifies to a tri-local quantity of 
leading twist-three which can be written in terms of three Lorentz invariant functions of 
positive parity V (vector ^ l^C), A (axial vector ^ 7^756*), and T (tensor ^ cr^^C) ||115 



{0K{zi)ul{z2)d;{z,)\p,+)- 



,abc 



fN[{^C)a/3{l5N)^V{Zi-p) + {^j5C)a(3N^A{z, ■ p) 

(75) 



-i(^t,uPuC)al3ilf,l5N)jT{z, -p)] . 



Here \P, +) is the proton state with momentum P and positive helicity, N denotes the proton 
spinor, and C = ^7072 is the charge-conjugation matrix. The "proton decay constant" /n is 
a dimensionful quantity determining the value of the nucleon distribution amplitude at the 
origin |3gj31|]. 

To perform the actual calculations of the nucleon distribution amplitudes, it is convenient 
to define the functions V, A, and T in the space of longitudinal momentum fractions: 



F{xi,X2,X3) 



A d{z, ■ P) 
11 — — — exp 

i=l 



2tt 



P) 



F{z, ■ P) 



(76) 



Note that these functions depend implicitly on the factorization scale which is supposed 
to serve as the starting point of evolution. Furthermore, renormalization group controlled 



evolution fixes the asymptotic limits of these functions to be lim^ 



V = lim„ 



T 



0as = 120x1X2X3, and lim^2^oo A = 0, the latter because of the Pauli principle. The next 
step is to obtain solutions for these functions at finite evolution scales by determining the 
coefficients Bn and then to represent them as a series expansion over a finite number of their 
eigenf unctions. 

Before we proceed, let us first discuss the symmetry properties of these functions 
in the infinite momentum frame. Because of the identity of the two u-quarks in the 
nucleon, the functions V and T transform under the scalar symmetric representation 
and the function A under the scalar antisymmetric representation of the permutation 
group 5*3. Then one has (in obvious shorthand notation ||116|| : P123 = [123], 



[213]y(l,2,3) = \/(l,2,3) 
[213]T(1,2,3) =T(1,2,3) 
[213]A(1,2,3) = -A(1,2,3). 



etc.) 



(77) 
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On the other hand, the requirement that the total isospin of the three-quark bound state 
should equal 1/2 yields one more relation: 

2T(1, 2, 3) = V{1, 3, 2) - 3, 2) + K(2, 3, 1) - A(2, 3, 1) . (78) 

Thus the system of the functions V , A, and T is redundant and the whole information they 



contain can be expressed by a single function, termed $Ar ||116| , |30 



^^{xi) = V{xi)-A{x,) . (79) 

This function is mixed symmetric, i.e., it transforms under the 2-dimensional (matrix) rep- 
resentation of 5*3. Knowing $jv it is sufficient to determine all other nucleon distribution 
amplitudes. For the time being, there is actually no method of computing hadron distribu- 
tion amplitudes as a whole^ Hence, we must content ourselves with descriptions in terms 
of truncated eigenfunctions series demanding that they (1) comply with theoretical con- 
straints, set for example by QCD sum rules or lattice calculations, (2) are renormalization 
group controlled, i.e., satisfy the nucleon evolution equation at every order of truncation, 
and (3) eventually match the available experimental data. 

In the following exposition the nucleon distribution amplitudes will be determined via 
their moments 

^(nmana) _ \dx]x1^ X^^ X^^ F {xi, X2, X^) , (80) 
JO 

where F{xi) stands for one of the amplitudes V, T, A, or linear combinations of them. These 
moments are related to the covariant derivatives of the invariant functions with respect to the 
light-cone positions Zi ■ p. The latter are the Fourier conjugate variables to the longitudinal 
momenta x,-. Hence we have 



3 



F{xi 



= {Z ■ P)'^l+"2+r^3 / [dx]x1'x':^'x'i'' F{Xi) 

Jo 

= {z ■ P^'^l+^z+'iS p{nin2n3) (81) 

where z is an arbitrary auxiliary vector (z^ = 0) and the caret, " , serves to distinguish 
operators from scalar functions. Taking matrix elements between the proton state and the 
vacuum, we get 

^Q|^(„l"2n3)^0)|P) = h{z ■ p)^'i+"2+n3+ljY^^(nin2n3) ^ (g2) 



where the moments on the rhs are defined in Eq. (|80D . More precisely, we consider matrix 
elements (moments) of the following operators: 

y(n,n,ns)^Q^ ^ . /})"i«(0) [C 7^ (^Z ■ D)"^m(0) {zZ ■ DY'^,d^{id) , (83) 



lOA 

promising approach, applied so far only to the pion, was recently presented in 
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(84) 



(85) 



The determination of the moments of these operators derives from correlators of the generic 
form 



I^''™''^\q,z) =i I rf^xe'«-"(0|T(F("^"2"^)(0)4'"^(x))|0)(z-7)^^, 



■ q) '^i+'^2+"3+ni+3j(nin2n3,m)^^2-j 



(86) 



where z is again a hght-hke auxihary vector and the factor {z-'y)^y serves to project out the 
leading-twist structure of the correlator. The computation of the Wilson coefficients on the 
quark side of the correlator amounts to the perturbative evaluation of diagrams involving 
local quark/gluon condensates PD[. It yields the theoretical side of the sum rule. The 
hadronic (phenomenological) side of the sum rule is obtained by saturating the correlator 
by the lowest-mass baryon state(s) via a dispersion relation. Reconciliation of the two sides 
of the sum rule with respect to the Borel parameter within a continuous interval as large as 
possible determines the variation of permissible values for a particular moment. We do not 
derive sum rules here. For more details we refer to the original works |30J3^J33J33] and the 



review article |3^. We merely use them as theoretical constraints imposed on the moments 
of the nucleon distribution amplitude. However, these constraints cannot fix the shape of the 
nucleon distribution amplitude uniquely and one has to impose additional global constraints 
as opposed to the moment sum rules which constitute local constraints.0 

The moments of the nucleon distribution amplitude $Ar in terms of the expansion coef- 
ficients Bn are displayed in Table |^. 

Because of longitudinal momentum conservation, xi + X2 + = 1, not all the moments 
at a given order M = rii + n2 + n-^ are linearly independent. This implies ||3^ 



(jii ,712,713) 

N 



^(rM+l,n2,?i3) ^(ni,?i2+l,n3) ^ 



(ni ,712,13+1) 



(87) 



For instance, at order M = 3 there are 20 moments out of which only 10 are strict. Which 
combinations are actually taken, depends on the choice of the polynomial basis in which 
the eigenfunctions are finally expressed. We use throughout this report the powers of the 
monomial X1X3, i.e., the basis \kl). In terms of this basis, the moments of $Ar read 



(ni?i2n3) 
iV 



1 rl—Xi 

dxi / dx3 
Jo 



EE(-ir 

i=0 j=0 




ni+i-j n-i+j 



^n{Xi,X3] 



(88) 



It was first shown in ||49|j50|JT0T[| and then outlined in [ p.03|| that it is possible to derive a 
closed-form relation between the expansion coefficients Bn and the strict moments of $Ar, 
defined by 



^^One can always add some oscillating function which vanishes at the points fixed by the local 
constraints but which contributes outside. 
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Table 2. Analytical expressions for the moments of the nucleon distribution amplitude 
^(nin2n3) ^ ^[fc] -^^ terms of the expansion coefficients B„ up to order M = ni + n2 + ns = 3. 



Moments 



Bo 

7B0+B1+B2 
21 

7B0-2B2 
21 

7B0—B1+B2 
21 

108 Bp +27 Bi+27 ^2+9 B3--B4--B5 
756 

18BQ-9B2+B3+Br, 
126 

108 -Bo -27 Bi+27 B2+9 -B3+B4--B5 
756 

72 -Bii +9 £?i - 9 -B2 - 3 -B3 +B4 -3 Br 
750 

36/j()+9B2-3/ji+2-B5 
378 

72 Bo-9 -51-9^2-3 -63-54-3.65 
756 

87120 Bq +29040 .Bi +29040 -B2 + 1 7424 B3 - 1936 B4 - 1936 .65 

1219680 , , 

-33 B6-33 V97 B6-33 B7+33 \/97 B7+I46 B8+2 V4801 Bs+146 B9-2 V4801 Bg 

1219680 

165 Bp -110 B2+22 B3+22 B5+B6+B7 
2310 

87120 Bo-29040 Bi +29040 B2+17424 B3+1936 B4-1936 B5 

1219680 , , 

-33 B6-33 V97B6-33 B7+33 V97B7-146 B8-2 V4801 B8-146 Bq+2 V4801 B9 

1219680 

5940 Bo+1320 Bi+88 B4-44O B5+I5 Bs+S VWf Be+15 B7-3 ^97 B7-I2 Bs-12 B9 

166320 

130680 Bo+14520 Bi +43560 B2-8712 B3-968 B4+10648 B5 

3659040 , , 

-231 B6+33\/97Bs-231B7-33 a/97 B7-174Bs-6V4801Bs-174B9+6 V4801 B9 

3659040 

495 Bn+55 Bi-165 B2-II B3+II B4-II B5-3BB-3B7+BS+B9 

13860 

130680 Bo-14520 Bi +43560 B2-8712 B3+968 B4+10648 B5 

3659040 , , 

-231 Bfi+33 V97Bs-231 B7-33 V97B7+174 Bs+6 V4801 Bs+174B9-6 V4801 B9 

3659040 

495B0-55 B1-I65B2-IIB3-IIB4-IIB5-3B6-3B7-B8-B9 

13860 

5940 Bn-1320 Bi-88 B4-440 B5+I5 B6+3 \/97 Bfi+15 B7-3 ^97 B7+I2 B8+12 B9 

166320 
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(iOj) 
N 



[dx] X\xl xi <l>Ar(Xfc,/i^ 



(89) 



to arrive at 



120 



ln(/iVA^^^) J 



1 -In 



OO 



(90) 



i,j=0 



where the normahzation constants Nn, the matrix coefficients a^-, and the anomalous dimen- 
sions 7„ (up to M = 4) are those tabulated in Table |I]. This expression enables the analytical 
calculation of the coefficients to any desired order of polynomial expansion. The utility of 
Eq. ( pUD is twofold: (1) As repeatedly stated, the moments of hadron distribution amplitudes 
are not accurately determined. Thus, without such an explicit relation between expansion 
coefficients and strict moments, one has to perform a simultaneous and self-consistent fit 
to the moment sum rules constraints, a procedure which obviously becomes increasingly 
tedious as the moment-order grows. (2) As already outlined in the previous chapter, the 
Hilbert-Schmidt orthogonalization procedure of polynomials with more than one variable 
is not unique, and one is well-advised to look for higher-order eigenfunctions which are as 
simple as possible, absorbing numerical coefficients into the normalization constants. This 
non- uniqueness entails that one can compare expansion coefficients 5„, obtained in different 
approaches, only if they are normalized. Without the knowledge of the normalization con- 
stants Nn, the values of Bn are of no significance or practical usefulness. However, having 
derived Eq. (|90|) , the knowledge of the specific normalization used in different approaches 
becomes superfluous. The coefficients i?„ can be self-consistently computed on the basis of 
the strict moments alone, which are universal quantities, modulo an overall normalization 
through the value of ^j^''^''. In the following. 



[dx] $Ar (xi, fi' 



1 



(91) 



is used. Hence, evaluating the strict moments of the nucleon distribution amplitude up 
to a given order, one can determine the corresponding expansion coefficients of the same 
order and vice versa. All these advantages will become successively apparent through the 
applications of the formalism in subsequent chapters. The link between i?„ and the strict 
moments of $7v is exemplified below (see also |41,|5l|): 

1260 



B.ifi'] 



120 
420 

120 
756 



$(100) _ $(001) 



2$r^ - 3^T^ 



N 



3$(ooi) 



120 
34020 

120 

1944 

120 



2$ 



(000) 
N 



(100) 
N 



7$ 



(100) 

N 

(001) 
N 



N 



■N 



(200) 
N 



+ 3*. 



(002) 



2$r^ - 7$!^°°) 



N 



7$(ooi) 



14 



N 



(200) 
N 



14 



N 



(002) 
N 



(92) 



Note that Bq is fixed to unity by the normalization of (cf. Eq. (pl|)). Furthermore, 



recall once again that the notations of |3^ are used. The next-to- leading order expansion 
coefficients (M = 3) are: 
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4620 485 + 11V97 



120 

31 
+3— 



97 



97. 



(101) 
N 



(000) 

17 + \/97 



6 f$r^ + ^T') 



41 + V97 
4 



($(200) ^ ^(002)) 



16 



(300) 



_^31-y97 ^^(,01) _^ $(102)) 



4620 485 



iV 

llV97 
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+3 



97 



(000) 



N 



31W97 (101) 



17 
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N 
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($(300) _^ ^(003)) 



_^3lW97 (^goi) _^ ^(102)) 



27720 33607 - 247V480T 
^20 4801 

601 + VisoT 



(100) 
N 
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(200) 
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59- V480T 
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$(300) _ $(003)\ 
$(201) _ $(102)\ 
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($(201) _ ^(102)) 
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The values of the moments of $Ar are extracted from the correlator in Eq. ( pq) for ni + n2 + 
na < 3 and m = 1 at some self-consistently determined normalization point /i = /ip of order 
1 GeV (see, e.g., p0|j35[] ) at which a short-distance operator product expansion can be safely 
performed and quark-hadron duality is supposed to be valid. Table ^ shows the range of 
values obtained by COZ |Q up to order M = 3 together with those calculated by KS ||3^ 
for the first- and second-order moments. 

This exposition completes our discussion of the expansion coefficients and the brief re- 
view of the conceptual essentials underlying their non-perturbative determination. Before 
moving on in the next chapter to the actual models for the nucleon distribution amplitude, 
we briefly discuss now the momentum evolution of the normalized expansion coefficients 
Bn (recall Eq. ^) against the order of expansion in terms of eigenf unctions. The result is 
illustrated in Fig. |1^. One observes that the ratio ^"^^^^j decreases both with increasing poly- 
nomial order and increasing momentum transfer so that a truncation at low orders seems, 
in principle, justifiable. However, in order that the truncation of the infinite eigenfunctions 
expansion, given by Eq. (0), is computationally useful, we have to ensure dominance of 
the lowest-order contributions. Given this premise, the guiding principle is to minimize the 
influence of higher-order terms which we have discarded. This procedure is in some sense 
analogous to the optimization of the renormalization-scheme dependence of physical quan- 
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Table 3. Numerical values of the moments M = ni + n2 
distribution amplitude in comparison with those of previous 



constraints derived by Chernyak, Ogloblin, and Zhitnisky [38 



+ ns < 3 of the heterotic nucleon 
models versus the QCD sum-rule 
and King and Sachrajda [p7[. 





COZ-SR 


KS-SR 


I (nira2n3) 
CD ^ ^ / 

^N/het 


^N/COZ 


^N/CZ 


^(nin2n3) 
^N/GS 


I (nin2ra3 
^N/KS 


(000) 


1 




1 


1 


1 


1 


1 


1 


(100) 


0.54- 


-0.62 


0.46—0.59 


0.5721 


0.5790 


0.630 


0.6269 


0.550 


(010) 


0.18- 


-0.20 


0.18—0.21 


0.1837 


0.1920 


0.150 


0.1371 


0.210 


(001) 


0.20- 


-0.25 


0.22—0.26 


0.2442 


0.2290 


0.220 


0.2359 


0.240 


(200) 


0.32- 


-0.42 


0.27—0.37 


0.3380 


0.3690 


0.397 


0.2879 


0.350 


(020) 


0.065- 


-0.088 


0.08—0.09 


0.0661 


0.0680 


0.0233 


0.0321 


0.090 


(002) 


0.09- 


-0.12 


0.10—0.12 


0.1696 


0.0890 


0.080 


0.0079 


0.120 


(110) 


0.08- 


-0.10 


0.08—0.10 


0.1386 


0.0970 


0.110 


0.1080 


0.100 


(101) 


0.09- 


-0.11 


0.09—0.11 


0.0955 


0.1130 


0.123 


0.2309 


0.100 


(Oil) 


-0.03- 


-0.03 


unreliable 


-0.0210 


0.0270 


0.017 


-0.0029 


0.020 


(300) 


0.21- 


-0.25 




0.2101 


0.2445 


0.257 


0.1281 


0.2333 


(030) 


0.028- 


-0.04 




0.0392 


0.0381 


0.0013 


0.0169 


0.0573 


(003) 


0.048- 


-0.056 




0.1392 


0.0485 


0.0413 


-0.0515 


0.0813 


(210) 


0.041- 


-0.049 




0.0789 


0.0587 


0.068 


0.0463 


0.0593 


(201) 


0.044- 


-0.055 




0.0490 


0.0658 


0.0713 


0.1135 


0.0573 


(120) 


0.027- 


-0.037 




0.0504 


0.0243 


0.0253 


0.0278 


0.030 


(102) 


0.037- 


-0.043 




0.0372 


0.0331 


0.0353 


0.0836 


0.0320 


(021) 


-0.004- 


-0.007 




-0.0235 


0.0056 


-0.003 


-0.0127 


0.0027 


(012) 


-0.005- 


-0.008 




-0.0068 


0.0073 


0.003 


-0.0241 


0.0067 
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Figure 10. Evolution behavior of the normahzed expansion coefhcients Bn versus the tradeoff with the 
order of expansion in eigenfunctions. 



titles computed in a perturbative scheme [principle of minimum sensitivity ||1 1 7|| ) . Without 
enough understanding of the underlying non-perturbative dynamics, it remains a challenge 
to develop a method of computing hadron distribution amplitudes as a whole. For the time 
being, we must content ourselves with an effective description in which the local low-order 
moment constraints are supplemented by additional global constraints to impose restrictions 
on the shape of the nucleon distribution amplitude as a whole. This may be achieved my 
means of a hierarchical (with respect to moment order) x^-criterion (see next section) or by 
demanding, for instance, smoothness of the nucleon distribution amplitude |^ . In this way. 



one can enforce the dominance of the lowest-order contributions and minimize the influence 
of the disregarded higher-order terms. Since the extraction of reliable estimates of higher- 
order moments from QCD sum rules is severely limited |^ (the accuracy of the moment 
values one can extract decreases as the moment order increases), there is actually no other 
pragmatic alternative to this type of approach. So, if one is satisfied with a given accuracy (of 
observables calculated with these nucleon distribution amplitudes) relative to existing data, 
reconstructing an analytic representation of the nucleon distribution amplitude in terms of 
higher and higher eigenfunctions is unnecessary - maybe even irrelevant. 



E. NUCLEON DISTRIBUTION AMPLITUDES 

It is evident from the discussion in the previous chapter that the method of QCD sum 
rules enables the determination of the basic characteristics of the nucleon distribution ampli- 
tude, namely: (1) its value at the origin, /n, and (2) a model representation over longitudinal 
momenta via a finite number of moments. However, the reliability of this procedure depends 
crucially on the particular way of using the QCD sum-rule constraints. The basic assump- 
tions in truncating the infinite series of eigenfunctions should be: 

• A truncated representation at relatively low orders provides a reasonable approxima- 
tion of the true nucleon distribution amplitude, given that the normalized expansion 
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coefficients decrease with the (moment) order, albeit slowly. 



The intrinsic inaccuracy of the sum rules for higher-order moments is taken into ac- 
count. 

The evolution equation is always satisfied. 



The first premise says nothing about the convergence of the eigenfunctions series. We 
know that this series converges for very large values of to the asymptotic solution 0as, 
but only logarithmically. Nevertheless, it seems reasonable to believe that the characteristic 
properties of the nucleon distribution amplitude belong to the entire series and that they 
do not ffist show up in higher orders. On the contrary, the polynomial structure of the 
eigenfunctions (comprising not only symmetric but also antisymmetric Appell polynomials) 
will inevitably introduce oscillations ("wiggles") at every step of truncation because of the 
normalization condition imposed on $Ar (cf. Eq. (0)). These wiggles are unphysical and 
should be washed out by destructive interference at subsequent steps of truncation. Hence, 
it is evident that the trial nucleon distribution amplitude at every step of truncation should 
be chosen to have a shape as smooth as the sum-rule constraints on its moments allow. This 
procedure has close resemblance to the Tamm-Dancoff method | |118 | of truncating time- 
ordered products in quantum field theory in seeking for a ground-state solution.lli In this 
method, as in our approach, the parameters of the low order solutions are effective, meaning 
that they incorporate by construction crucial high-order effects. 

Following this strategy, we introduce a x^-criterion which not only parameterizes de- 
viations from the sum rules laterally, i.e., within some fixed moment order, but - in ad- 
dition - which weighs the quality of sum rules vertically, i.e., according to their order 
P,P|J^,|5^P5|. This can be realized by defining 



xl 



2 I 2 
Xk,{a) + Xk,{b) 



{nik 



(94) 



with 



Xfc,(a) 



mm 



(95) 



or 



|^max| 



(mfc [k = 1, 2, ... , 18) denoting collectively the moments), where Nk 
whether lies on the left- or on the right-hand side of the corresponding sum-rule interval 

iXtot = T.k xl) and 



Xfc,(fe) 



100, 1 < A; < 3 
10, 4 < < 9 
1, lQ<k<li 



(96) 



Let us consider a concrete example in order to make this criterion and its use more trans- 
parent. Suppose a trial distribution amplitude corresponds to Xtot = 132.85. This means 
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I wish to thank Prof. D. Shirkov for discussions on this point. 
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Table 4. Theoretical parameters to classify the nucleon distribution amplitudes discussed in 



the text. The samples shown refer to the COZ sum rules. 



Model 




Bi 




B2 




^3 


B4 




B5 




R 




Symbol 


Het 


3 


4437 


1 


5710 


4 


5937 


29.3125 


-0 


1250 


-1.89 


.104 


33.48 


• 


Het' 


4 


3025 


1 


5920 


1 


9675 


-19.6580 


3 


3531 


24.44 


.448 


30.63 


• 


coz°p^ 


3 


5268 


1 


4000 


2 


8736 


-4.5227 





8002 


9.13 


.465 


4.49 


■ 


COZ^'P 


3 


2185 


1 


4562 


2 


8300 


-17.3400 





4700 


5.83 


.4881 


21.29 


+ 


COZ 


3 


6750 


1 


4840 


2 


8980 


-6.6150 


1 


0260 


10.16 


.474 


24.64 


□ 


cz 


4 


3050 


1 


9250 


2 


2470 


-3.4650 





0180 


13.40 


.487 


250.07 


♦ 




3 


5818 


1 


4702 


4 


8831 


31.9906 





4313 


-0.93 


.0675 


36.27 





KS/COZ°P^ 


3 


4242 


1 


3644 


3 


0844 


-3.2656 


1 


2750 


9.47 


.453 


5.66 







3 


5935 


1 


4184 


2 


7864 


-13.3802 


2 


0594 


13.82 


.482 


40.38 





KS 


3 


2550 


1 


2950 


3 


9690 


0.9450 


1 


0260 


2.47 


.412 


116.35 





QgOpt 


3 


9501 


1 


5273 


-4 


8174 


3.4435 


8 


7534 


80.87 


.095 


54.95 


▲ 


Qgmin 


3 


9258 


1 


4598 


-4 


6816 


1.1898 


8 


0123 


80.19 


.035 


54.11 


▼ 


GS 


4 


1045 


2 


0605 


-4 


7173 


5.0202 


9 


3014 


78.87 


.097 


270.82 


A 


Samples 





3 


3125 


1 


4644 


3 


1438 


-1.0000 





8750 


7.67 


.441 


4.63 


+ 


1 


3 


2651 


1 


4032 


3 


5466 


2.8685 


1 


7954 


8.94 


.405 


5.11 


+ 


2 


3 


4026 


1 


4917 


3 


0629 


7.3430 





6719 


8.75 


.385 


16.07 


+ 


3 


3 


7225 


1 


5030 


3 


6592 


10.7265 


1 


5154 


9.29 


.355 


17.78 


+ 


4 


3 


8407 


1 


4968 


3 


2142 


14.4093 





8757 


10.49 


.325 


19.41 


+ 


5 


3 


6544 


1 


4000 


3 


0993 


15.5614 


-0 


1329 


6.35 


.305 


18.15 


+ 


6 


3 


8607 


1 


4000 


3 


2375 


19.8571 


-0 


1635 


6.32 


.255 


20.57 


+ 


7 


3 


9783 


1 


4000 


3 


2706 


22.4194 


-0 


4805 


5.29 


.225 


21.69 


+ 


8 


4 


1547 


1 


4000 


3 


3756 


26.1305 


-0 


5855 


5.02 


.175 


23.52 


+ 


9 


3 


4044 


1 


5387 


4 


3094 


25.5625 





0625 


.01 


.153 


30.80 


+ 



that it breaks one first-order sum rule, three second-order sum rules, and two third-order 
sum rules, whereas the total breaking is 85%. 

This "hierarchical" treatment of the sum rules (1) accounts for the higher stability of 
the lower- level moments relative to the higher ones and (2) does not overestimate the 
significance of the unverified constraints for the third-order moments. Consequently, our 
parameterization in Eq. (|96D favors solutions that satisfy best the lowest-order moments and 
gives (arbitrary) penalty points to those solutions which violate them. This relegates third- 
order terms to contribute only marginally, i.e., to merely refine (if at all) the shape of the 
nucleon distribution amplitude determined in the first step solely on the basis of the second- 
order moment constraints. It goes without saying that this procedure intrinsically suppresses 
solutions with unphysical oscillations. A physical analog of truncating the representation 
of the nucleon distribution amplitude this way is perhaps provided by a holographic image 
which is not destroyed when cut into smaller pieces (corresponding here to a lower order of 
truncation) but becomes rather less sharp |^T]]. In contrast, a series truncation according 
to a simple (i.e., lateral) x^-criterion corresponds to a conventional image which is really 
destroyed when cut into smaller pieces. 
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Table 5. Theoretical parameters to classify the nucleon distribution amplitudes discussed in 



the text. 


The sampl 


3S shown refer to a 


combined use 


of the COZ and KS 


sum 


rules. 




Model 


Bi 


B2 


Bs 




B5 


^[deg] 


R 


x' 


Symbol 


Samples 





3.7520 


1.3845 


2.7000 


-10.5000 


1.8500 


14.24 


.480 


28.77 





1 


3.7065 


1.3298 


2.9000 


- 8.0500 


1.1500 


10.15 


.475 


18.28 





2 


3.4075 


1.4191 


3.3813 


- 7.6500 


1.4750 


7.82 


.465 


9.02 





3 


3.6695 


1.3186 


2.8375 


- 5.8875 


1.1125 


10.62 


.470 


16.90 





4 


3.4120 


1.3906 


3.2375 


- 5.5875 


1.3625 


8.51 


.460 


6.76 





5 


3.4242 


1.3644 


3.0844 


- 3.2656 


1.2750 


9.47 


.453 


5.66 





6 


3.3500 


1.3710 


3.1192 


- 0.9556 


1.2995 


9.48 


.440 


5.75 





7 


3.3501 


1.4045 


3.3004 


0.2513 


1.4203 


8.95 


.430 


5.86 





8 


3.3536 


1.4327 


3.4787 


4.8241 


1.5609 


8.97 


.395 


15.11 





9 


3.3500 


1.3303 


3.1262 


8.8918 


1.3173 


10.75 


.360 


15.19 





10 


2.9067 


1.3664 


4.0326 


12.0701 


1.0180 


2.46 


.300 


26.37 





11 


2.9300 


1.3899 


4.4263 


16.4589 


0.8764 


-0.019 


.250 


28.91 





12 


3.1760 


1.4491 


4.6009 


22.9637 


0.5213 


-1.016 


.180 


31.06 





13 


3.2912 


1.4545 


4.6802 


25.5161 


0.4903 


-1.014 


.150 


32.96 





14 


3.4017 


1.4595 


4.7558 


27.9607 


0.4594 


-1.014 


.120 


34.20 





15 


3.5078 


1.4638 


4.8284 


30.3070 


0.4310 


-1.012 


.090 


35.39 






All told, let us present the results. They are gathered in Tables ^ and ^ in correspondence 



with Figs. 11 and 12. We can see the utility of this procedure of fitting the sum rules as 



follows: the parameter space of the decomposition coefficients ( "Appell" coefficients) that 
project on to the nucleon eigenf unctions, derived in the previous chapter, is systematically 
scanned seeking for evolving morphologies of geometric features such as the peak pattern 
of the amplitudes V, A, and T. Quantification of a feature or region involves its isolation, 
extraction, and stability. Next, one tracks features corresponding to particular solutions to 
determine how they change their form and how they compare with features in other regions. 
By testing the effects of various combinations of symmetric versus antisymmetric components 
of the distribution amplitudes, the whole parameter space can be systematically explored 
with respect to local minima of and a complete picture of their distribution pattern can 
be pieced together. Using for the first- and second-order moments either the COZ or the 
KS sum-rule constraints in conjunction with those of COZ for the third-order moments, a 
simple scaling relation between the ratio R = and the expansion coefficient 

emerges as one progresses through the generated solutions. That smooth scaling behavior 
is inherent in the particular set of sum rules used in the fit and does not rely on additional 
assumptions. 



Fig. 12 illustrates the x quality of the COZ samples. Depending on the degree of 



matching with the corresponding sum rule, a solution found this way appears as a local 
minimum of the criterion. We have plotted in the (-B4, R) plane interpolating solutions 
to the COZ sum rules (+ labels) and such to a combined set of KS/COZ sum rules (o labels). 
The latter set is not quite consistent because the typical sum-rule parameters, such as Borel 
intervals, energy thresholds, etc., are treated differently in the two approaches. Nevertheless, 
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Figure 11. Pattern of nucleon distribution amplitudes matching the COZ sum-rule constraints (crosses) 
and such to a combined set of KS/COZ sum rules (open circles). A striking scaling relation between the 
ratio R of the magnetic nucleon form factors and the coefficient B4 that projects on to the corresponding 
eigenfunction is revealed. The various models of nucleon distribution amplitudes discussed in the text are 
compiled in Table ^. Interpolating samples of distribution amplitudes with respect to the COZ sum rules are 
denoted by crosses (see Table ^); those referring to the combined set of the KS/COZ sum rules, by circles 
(see Table 1^). The inset at the lower left expands the vertical scale between 0.455 and 0.495 (corresponding 
to the B4 interval [—10,0]) to exhibit the close agreement between the fiducial orbit and a variety of 
proposed amplitudes with third-order Appell polynomials, not included in the fit. Significantly, those model 
amplitudes which appear as isolated points are exactly the ones which possess unphysical features, namely, 
either in the form of spurious oscillations, or because they violate the renormalization group equation, 
leading this way to a wrong QCD evolution behavior of form factors (see [^). The curves are fits to the 
local minima of the COZ sum rules (solid line) and the KS/COZ sum rules (dotted line). They constitute 
a fiducial orbit with respect to beginning in the heterotic region (small R and large positive -B4) and 
terminating past the COZ cluster (large R and large negative B4). 



it is instructive to examine what changes the inclusion of the KS sum rules may cause. As 
it turns out, there is no significant difference between the two treatments (see Fig. |ll]), and 
this insensitivity justifies the whole approach and promotes the orbit structure to a key 
element of the analysis. 

The solutions to the sum rules agree within less than 1% with the empirical fits : 

R = 0.437338 - 0.006016^4 - 0.000176B| (97) 

(only COZ constraints) or 

R = 0.431303 - O.OO752S4 - 0.000241S| + 3.851221 x lO'^Bl (98) 

(combined KS and COZ constraints), represented by the solid and dotted curves, respec- 
tively. For simplicity, we use the term "orbit" to refer to both characteristic curves. 

Furthermore, we find that the magnitude of the ratio R is restricted at both ends of the 
orbit. For the first discussed case, the upper bound is 0.4881 and for the second one, 0.482. 



The corresponding lower bounds are 0.104 and 0.0675. In the course of this analysis |51,52 
optimized versions - with respect to the sum rules - of all previous models P(]|J53|P7|J55 
have been determined. These amplitudes, denoted by the superscript "opt", are shown in 
Fig. |13[ [The assignments of models to symbols are given in Table ^.] 
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Figure 12. Two-dimensional view on solutions to the QCD sum rules as successive local minima in the 
plane spanned hy R = /G^ and B4 with respect to the "hierarchical" x^-criterion discussed in the 

text. Referring to the previous figure, only the COZ branch of the "fiducial orbit" is shown. 

As mentioned previously, the lower part of the orbit is associated with the heterotic 



solution which gives the smallest possible ratio still compatible with the sum-rule con- 
straints. This solution, although degenerated with respect to R, is distinct from the GS one. 
As we shall see in the next chapters, the predictions extracted from the heterotic model 
are dramatically different compared to those following from the GS model. The upper re- 
gion of the orbit controls COZ-type amplitudes and contains a cluster of solutions densely 
populating the orbit in the i?-interval 0.455 -7- 0.495 (see the inset in Fig. |TT|). This cluster 
contains the amplitude denoted COZ°p* (see Table ||), which is associated with the abso- 
lute minimum of ^^^^ plays the role of an attractor for all other solutions with similar 
features. Strictly speaking, also the CZ model and the KS one, although in the vicinity 
of the orbit, are actually isolated points because they correspond to much larger values 
(cf. Fig. |l^) in correspondence with Table ^) than the proper solutions on the orbit. The 



heterotic amplitude and the original COZ amplitude correspond to similar values. It 



is remarkable that the heterotic solution matches the KS |^ sum-rule constraints better 
than the original COZ amplitude. This is worth noting because the KS results have been 
independently verified in . 

Analogously, the solution with the smallest value on the KS/COZ branch of the orbit 
is the amplitude denoted KS/COZ"^^ (see Table |). GS-type amplitudes constitute an 
isolated region in the (i?4, R) plane and are separated from the characteristic orbit by a 
large barrier. There is yet another type of solutions (we termed Het') past the upper end 
of the orbit - first discussed in [^. Het' belongs to a secondary-branch of the main orbit 



which can be reached by coercing the amplitudes V and T to have a particular symmetry 
pattern of maxima and minima, corresponding to the reversed combination as compared to 
the heterotic solution. Remarkably, each of the discussed categories of solutions has unique 
geometric characteristics specific for its class. One has to exercise a certain amount of care 
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R = 0.465 







Figure 13. Profiles of tlie optimized versions of model distribution amplitudes, commonly discussed in 
the literature, shown together with the heterotic one. 

to be sure that these are indeed the only regions in the parameter space contributing to the 
sum rules at the desired level of accuracy. 

Having completed the classification of distribution amplitudes rendering the degree of 
Appell polynomials M < 2, let us now turn to models with functional representations 
which attempt to include third-order eigenfunctions, in connection with additional cutoff- 
parameters p]| , |119| ]. The inset in Fig. ^ shows how such models |^D| (marked by stars) 



and ||119|] (marked by light upside-down triangles) group around the optimum amplitudes 



COZ°P* and KS/ COZ°^^, thus establishing the scaling relation between R and in a much 
more general context. This result suggests that the inclusion of higher-order eigenfunctions 
(Appell polynomials) in the nucleon distribution amplitude is indeed a marginal effect, as 
argued above. If so, the antisymmetric content of the nucleon distribution amplitude is suf- 
ficiently accounted for by the eigenfunction $4(0:4) since higher antisymmetric components 
are offset by this term. On the other hand, those model amplitudes proposed in |40|jll9|| 
which appear as isolated points scattered toward the GS "island" are unacceptable on phys- 
ical grounds, either because they exhibit unrealistic large oscillations in the longitudinal 



momentum fractions |40] or because they entail a wrong evolution behavior for the nucleon 



form factors [ p.l9| |. In addition, as was pointed out in |1T[], models which use cutoffs, like 
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Figure 14. A two-dimensional example showing how the profile of the nucleon distribution amplitude 
$iV changes along the orbit. The amplitudes are labeled by the corresponding R value. 



those in [p.l9||, have difficulties in preserving the vahdity of the evolution equation. 



The general picture emerging from the presented analysis is a pattern of nucleon dis- 
tribution amplitudes which develops into several regions of the (i?4, R) plane. The main 
sector is characterized by an orbit which shows a striking scaling behavior between the ratio 
R (a measurable quantity) and the (theoretical) expansion parameter i?4. For a variety 
of amplitudes this result is unique, irrespective of their functional representation. Isolated 
samples in this plane are relegated to spurious solutions owing to their unphysical features. 
The proffies of the distribution amplitudes across the orbit change in an orderly sequence 
of gradations with some mixture of COZ and GS characteristics until the COZ amplitude is 
transmuted into the heterotic solution. The variation with shape of the nucleon distribution 



amplitude with R is shown graphically in Fig. |14. 

Despite differing functional representations of the nucleon distribution amplitudes (use 
of higher-order eigenfunctions with or without cutoffs), the results plotted in Fig. |Tl| show 
that all nucleon distribution amplitudes complying with the sum-rule constraints actually 
collapse to the fiducial orbit, thereby providing convincing evidence that the scaling relation 
between the ratio R and the expansion parameter i?4 reflects a genuine non-perturbative 
feature of the true nucleon distribution amplitude. To underline the hybridity character of 



the heterotic distribution amplitude, we show in Fig. 15 the invariant functions V ^ A, and 
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T associated with it in comparison with their counterparts for the optimized versions of 
the COZ and the GS models. To leading order M = 2, the algebraic expressions for these 
functions in terms of the expansion coefficients Bn are p5| : 

(5o + ^2 - 5^3 - 555) + ^ (fil - 3^2 + 11^3 + 54 + 2155) (Xi + X2) 



V{Xi) = (f)as{x 

- (5i + 54) X3 - (453 + 1455) X1X2 + I (I2S3 - 454 - 2855) (xj + 



X2 



1 
3 



+^ (24^3 + 4^4 + 1455) xl 



(99) 



A{xi) = 0as(a;i) 



- i-Bi - 3B2 + 3B3 - 54 - 7S5) (xi - X2) 



+- (-12^3 + 4^4 + 28^5) (xl - xl 
b ^ 



(100) 



T{Xi) = (f)as{Xi) 



{Bo + B2- 553 - 555) + (-3^2 + 753 + 755) 0:3 

(101) 



+ (4^3 + 1455) X1X2 + (8^3 + ^55) {xl + 



where the appropriate coefficients for every model have to be inserted.li^l While Vuet has a 
symmetry pattern similar to that of Vgs (one main maximum in the central region) , Tjjet is 
characterized by two maxima, much the same as Tqoz- The inverse heterotic combination 
is realized by the "mirror" solution (see Table 

These remarks can be put on more mathematical grounds by considering a classifica- 
tion scheme of the various nucleon distribution amplitudes based on the observation that 
the optimized versions of the COZ-type and GS-type amplitudes, we derived, are almost 
orthogonal to each other with respect to the weight w{xi) = (f)as{xi)/120. Their normalized 
inner product (COZ°p*|GS'°p*) yields 0.1607, which corresponds to an angle of 80.8°. Thus 
these functions form a quasi-orthogonal basis that can be used to continuously parameterize 
the nucleon distribution amplitudes in terms of a "hybridity" angle i), defined by 

where the index i denotes collectively any of the nucleon distribution amplitudes listed in 
Tables ^ and |^. Here the bracket denotes the normalized inner product 
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Note that at the central point xt = 1/3, 



120 

V {xi) = T {xi) = {xi) = — (27 - 12^3 - 2B5) 



to leading order M = 2. 
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Figure 15. Invariant distribution amplitudes V, A, and T associated with the optimized versions of 
the COZ and GS models in comparison with the heterotic model. 
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where 



{i,j)= f [dx]w{xi)<^^{xk)<^j{xk) (104) 



and 



(105) 



The hybridity angle parameterizes the minghng of geometrical characteristics attributed to 
COZ-hke and GS-like amplitudes and provides a quantitative measure for their presence 



in any solution conforming with the sum-rule constraints (see Fig. 16). The superimposed 



dashed line in Fig. 0(c) is a fit given by i^/[deg] = 8.5693 + 0.016054+0.00735|-0.00067B, 



The dashed line in Fig. ^6Ka) represents the empirical fit i? = 0.436415 — 0.005374S4 — 
0.000197i?|. An improvement of this fit was presented above in connection with the orbit 



structure of the nucleon distribution amplitudes (cf. Fig. |TlD. The mixing of different geo- 
metrical characteristics is particularly relevant for the heterotic solution, which generically 
amalgamates features of both types of amplitudes (cf. Fig. |T3]). In this role, the heterotic 
model has the special virtue of simultaneously fitting the twin hopes for making reliable pre- 
dictions with respect to the experimental data while being in agreement with the sum-rule 
constraints. The other models can be tuned to fit some aspects of the data, but never all 
aspects simultaneously. 

To conclude this chapter, we remark that fixing the value of the ratio R by experiment, 
one could use the presented classification scheme to identify the corresponding nucleon 
distribution amplitude determined by theory. 

F. DELTA DISTRIBUTION AMPLITUDES 

The heterotic concept developed for the nucleon distribution amplitude can be extended 
to derive an optimal distribution amplitude for the A"*" (1232) resonance as well [^]. The 



quark structure of this particle and the A^ — A transition form factor were investigated within 
the method of QCD sum rules by Farrar et al. "FZOZ" (denoted by the acronyms of 
the authors), and independently by Carlson and Poor (CP) [|^. Both groups derived sum- 
rule constraints on the moments of the A distribution amplitude and used them to construct 
model distribution amplitudes in terms of low-order eigenf unctions (M = 2). Note that only 
symmetric eigenf unctions under permutations xi X3 can contribute to the A distribution 
amplitude, so that Bi = B4 = 0. Thus the general form of the distribution amplitude for 
M = 2 is 



+ ^Bt) (xl + xl) + (-3B^ + 7B^ + W^) x. 



14 
3 

+ (4^3^ + 14^5^) xxx^ + (^o^^ + B^ - hB^ - 55 



(106) 
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Figure 16. Classification scheme of nucleon distribution amplitudes complying with the constraints set 
by QCD sum rules (Table I), (a) The ratio R = as a function of the expansion coefficient B4. 

The positions of the minima are indicated (see Fig. |l^). (b) Distribution of local minima of (on 
a logarithmic scale) plotted vs B4. (c) Pattern of nucleon distribution amplitudes parameterized by the 



hybridity angle i?, defined in Eq. (102), vs B4. The dashed curves are the fits described in the text. 
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Table 6. Numerical values of the moments ni + n2 + 71-3 < 3 of model distribution amplitudes 
for the A"'"-isobar in comparison with the sum-rule constraints of Carlson and Poor (CP) |45| and 



those of Farrar, Zhang, Ogloblin, and Zhitnitsky (FZOZ) |44|. The numbers in parentheses are 
those pubhshed by FZOZ. 



Moments 


Sum rules 




Models 






(^171271-3) 


TAiFZOZ) 


CP 


FZOZ 


heterotic 


FZOZ°P* 


(000) 


1 


1 


1 


1 


1 


(100) 


0.31—0.35 


0.350 


0.325 (0.32) 


0.321 


0.325 


(001) 


0.35—0.40 


0.300 


0.350 (0.36) 


0.357 


0.350 


(200) 


0.14—0.16 


0.160 


0.150 


0.140 


0.156 


(002) 


0.15—0.18 


0.123 


0.160 


0.151 


0.154 


(110) 


0.07—0.1 


0.101 


0.080 (0.07) 


0.078 


0.071 


(101) 


0.09—0.13 


0.089 


0.095 (0.1) 


0.103 


0.098 


(300) 


0.06—0.09 


0.085 


0.083 (0.085) 


0.073 


0.090 


(003) 


0.06—0.10 


0.060 


0.085 (0.081) 


0.071 


0.078 


(210) 


0.025—0.04 


0.039 


0.030 (0.025) 


0.027 


0.026 


(201) 


0.04—0.06 


0.035 


0.037 (0.04) 


0.040 


0.040 


(102) 


0.035—0.06 


0.031 


0.037 (0.039) 


0.040 


0.038 


VAiCP) 


(001) 


0.33—0.37 


0.350 


0.325 


0.321 


0.325 


(002) 


0.14—0.18 


0.160 


0.150 


0.140 


0.156 


(101) 


0.072—0.12 


0.095 


0.088 


0.091 


0.085 



Note also that Ta(xi, X2, 0:3) = [$a(132) + $a(231)] and that in addition Ta{xi, X2, X3) = 
<I'a(231) because $a(132) = $a(231). The shapes of these two model distribution ampli- 
tudes look quite different, and also the predictions one can extract from them are different 
(see for details in subsequent chapters). While the CP amplitude has almost asymptotic 
profile, the FZOZ one exhibits two maxima. The corresponding sum rule estimates are 
compiled in Table ^. The same table contains also the moments of the proposed model 
distribution amplitudes. 

Following again similar ideas of "heteroticity" , as in the nucleon case, one can attempt 
to combine the CP sum rules with those of FZOZ in order to obtain an amplitude that is 
capable to satisfy both sets simultaneously. The incentive is to obtain an amplitude that is 
optimized with respect to the sum-rule constraints as well as to improve its phenomenological 
capabilities. Both hopes can be satisfied by an amplitude termed again "Het", derived in 
and shown in Fig. This amplitude satisfies all FZOZ moment constraints while 



providing the best possible compliance with those of CP (only one of their sum rules is 
violated). Concerning its shape, one sees that it is, nevertheless, close to the CP amplitude, 
though its maximum is somehow shifted towards the center of phase space. 

There is yet another specific A amplitude denoted "FZOZ°p*, which was derived 



by demanding that the nucleon-Delta transition form factor, which involves in the nucleon 
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Figure 17. Profiles of model distribution amplitudes, defined in Table 0, for the A~''-isobar, derived 
from QCD sum rules. 

channel the COZ°p* distribution amphtude, is at least positive.0 Unfortunately, it turns out 
that this amplitude cannot reproduce any data (see below), though it is perfectly acceptable 
from the point of view of the sum rules. Thus we reiterate that optimum agreement with the 
moment sum rules - provided by COZ°p* - does not necessarily entail best phenomenological 
behavior. The expansion coefficients for all considered distribution amplitudes are listed in 
Table 0. 



^^All COZ-like distribution amplitudes yield N-A transition form factors compatible with zero. 



Table 7. Expansion coefficients for A-isobar distribution amplitudes derived from QCD sum 
rules within a truncated expansion of eigenfunctions. Note that = 1 by normalization and that 
only eigenfunctions symmetric under xi X3 contribute, i.e., = B^ = 0. As in the nucleon 
case, the notation of |35| is adopted. 



Expansion coefficients 


CP 


FZOZ 


Models 

heterotic 


FZOZ°P* 


Bt 


0.350 


-0.175 


-0.2499 


-0.175 


Bt 


0.4095 


1.071 


0.3297 


1.4117 


Bt 


0.1755 


-0.486 


-1.6205 


-1.620 
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IV. ELECTROWEAK FORM FACTORS 



The technical apparatus discussed in the previous chapters can now be used to calculate 
several electroweak form factors of the nucleon. 

The standard differential inclusive cross section for electron-nucleon scattering in terms 
of directly measured quantities is 



(Pa 



dVLdE' 4^2 sin^ ^ 



2W,{Q\v) sin2^ + iy2(Q',^)cos2^ 



(107) 



where E' = E — Q'^/2M^ and Wi and W2 are Lorentz- invariant structure functions of the 
target nucleon expressed in terms oiQ^, v = (5^/2Mjv = {E — E') is the energy loss of the in- 
cident lepton, and = (P' — P)'^ = —AEE'siii^ f ^ {Q^ = —q^ > 0) is the four-momentum 
transfer. If the nucleon (with mass Mn) were a point-like particle like the electron, say, we 
would have P^ = P'^ = and P ■ g = Mni^ so that u = Q^/2Mn; furthermore, Wi{Q^) 
and W2{Q'^) would reduce respectively to Q'^/AMn and 1. Then Eq. (|107|) would simply be: 



a 



AalE'' 



dQdE' 



2O 
cos - 



2Ml 



■ 2O 
sm - 



5 



2M, 



N . 



Integrating over dE' one gets the familiar formula 

da 

dn " 




Mott 




tan 



(108) 



(109) 



where the Mott differential cross section for a singly-charged point-like target, assumed 
to have an infinite mass, is cos^ ^ / AE'^ sin'^ |. However, because of the internal (quark) 
structure of the nucleon, this differential cross section does not describe the elastic scattering 
of an electron off the nucleon. Hence, one has to assume that the incident electron is scattered 
elastically by a nucleon with a "diffuse" structure which can be revealed by increasing the 
momentum transfer. Then we can write Eq. ( |108| ) in the form 



d^a AalE 



12 



dQdE' 



-6 



V — 



2M, 



X 



G|(g2) + (gV4M|,)G^(Q2 



1 + gv4M2, 

where GKQ^) and G^((5^) are the electric and magnetic ("Sachs") form factors of the 
nucleon, respectively, parameterizing in some sense our ignorance about its non-perturbative 
(binding) structure. A phenomenologically successful parameterization is provided by the 
empirical "dipole" form factor, viz: 



9 



cos 



+ 



2M2, 



e' 



sm 



(110) 



GUO) 



r<2 



(0) 



1 



O.riGeV 



(111) 



Obviously, the finite-size ("diffuse") nucleon causes the elastic scattering cross section to 
decrease rapidly as grows - in contrast to the behavior of a point-like nucleon described 
by Eq. (|108|) . Hence, as long as there is a single hadron in the final state, there is an 
extremely rapid drop-off of the cross section for » 0.7 GeV^. Clearly, the challenge is 
to calculate the nucleon's form factors within QCD, employing for the intact nucleon wave 
functions of the sort we discussed in the previous sections. 
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A. NUCLEON FORM FACTORS 



From Lorentz covariance, charge conservation, and invariance under space reflections 
there exist two electromagnetic form factors for a spin-1/2 particle. Then, in the Sachs 
parameterization, we have 



2M. 



N 



u{P) , (112) 



where r = {l + ^) , = (P^ + P^), and i?^ = f (7^. f 4- 4 Pit.) with normal- 
izations as G'e(O) = Cat (e^ = l,e„ = —1); G'm(O) = /i being the total magnetic mo- 
ment in units ejlM^: Hp = 2.7928, fin = —1.9131, and with the threshold condition 
G'e(4M^) = G'm(4M^). The connection to the Dirac parameterization is established by the 
on-shell relation 

u'{P')tRXP) = ^\P') i^iM^a^^q' + q\) u{P) (113) 



with 

Fi = T \ Ge H t-ptGm Dirac form factor 

F2 = T {Gm — Ge) Pauli form factor 

Ge = -^1 7-pr-^2 electric form factor 

Gm = Fi + F2 magnetic form factor (114) 

so that the on-shell nucleon current in the Dirac parameterization is given by 



{p'\j;-mp) = u'{p') 



Fi{r)l, + WT^F2{q')a,,q^ 



2M. 



N 



u{P) . (115) 



Recalling Eq. (^, the nucleon helicity-conserving form factor can be written in convolu- 
tion form as 

G'm(Q') = J^\dx]Jjdy]<^*{y„Ql)Tn{x„y„Q\a,{fi')) ^x„Ql) . (116) 

In the evaluation of this expression one may encounter singularities owing to the fact that 
the gluon virtualities should enter the arguments of the strong coupling constant in Th, i.e., 
as(/i^ = Q'^Xiyj). These virtualities depend on the longitudinal momentum fractions, and 
as a result, as becomes singular in the endpoint region. Several options have been suggested 
to avert this singularity: 



• The simplest choice is to assume that one can use two collective scales, one for the 
hard and another for the soft gluon propagator, and then take the geometric average 
of both couplings in the form factor. In this way the dependence of as on the frac- 
tional momenta is avoided and as can be taken outside the integral, amounting to a 
rescaled Q^-dependent pre-factor ("peak approximation"). It seems reasonable to fix 
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the gluon scales by the gluon virtuahties from that Feynman graphs which contribute 
most for each particular nucleon distribution amplitude by substituting the longitu- 
dinal momenta at the position of the main maximum of the corresponding nucleon 
distribution amplitude. Then one has as{Q'^) = [cts(Q^'?^hard)«s(Q^'ysoft)]^''^, where 
typically fhard = Xifji and fsoft = xji/j. This technique was used in several works, like 
[0,0,^,^,0, and the results are quite reasonable. 



A refined version of this approximation was proposed by Stefanis in . Recalling that 
there are in total 8 topologically distinct diagrams which give non-zero contributions 
to the form factor, one can use the virtuahties of the hard and soft gluon propagators 
in each individual diagram to calculate the corresponding as values. This procedure 
improves the predictions, though this improvement depends on the specific nucleon 
distribution amplitude used; for example, the effect when using the CZ-distribution 
amplitude is negligible 



In both cases considered above, one can only change the absolute value of the form 
factor but not its slope. This can only be improved if as can be retained inside the 
integral. To this end, one can either cutoff as explicitly - typically between 0.5 to 07 
- or saturate its momentum dependence by introducing, say, an effective gluon mass. 
Such an analysis was performed in ||120|| , having recourse to a modified expression for 



as, proposed by Cornwall ||121|| (see also ||122| 
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<Q") = . , r/^o . .X ,.0 1 ■ (117) 



/?oln[(Q2 + 4m2) /A^^cD 



Here stands for the effective gluon mass with values in the range (500 ± 200) MeV. 
Technically, is an IR-cutoff serving implicitly to regularize that gluon propaga- 
tor which becomes soft (compared to the large external momentum) in the endpoint 
region. It is evident that such a nonzero gluon mass will "freeze" the value of the 
running coupling constant for < Am^. Since a dynamically generated (gluon) mass 
is not a constant but vanishes at large momentum, the asymptotic behavior is not 
affected. Strictly speaking, one should use to saturate as a scale-dependent gluon 
mass which is governed by the renormalization group equation with an appropriate 
anomalous dimension [[L21|| . Such a saturation procedure was adopted in and will 



be presented below. Physically, the IR regularization of as connects to our previous 
discussion concerning the existence of a fundamental IR scale in the non-perturbative 
regime. This scale my be thought of as being the average transverse momentum of vac- 
uum partons and is of order 300 MeV to 600 MeV, but cannot be reliably computed 
at present. Within the framework of quark/gluon condensates, this scale separates 
hard from soft momentum flows of gluon (or quark) propagators which go into the 
corresponding condensates. It is also worth noting that according to Cornwall, rrig 
and Aqcd are interrelated by the consistency relation mg/AgcD ~ 1-5 — 2.0. It was 
warned in that the treatment in ||120|| violates this relation, handling rrig as an 



additional fit parameter for improving the slope of the form factors. 

Another, theoretically more appealing, possibility is provided by including in the cal- 
culation of the form factor gluon radiative corrections in the form of Sudakov-type 
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damping exponentials ||7^. This procedure is technically more complicated though 
theoretically more sound than, more or less plausible, cutoff or saturation prescrip- 
tions. The "Sudakov option" will be discussed in detail in the second main part of 
this report. 

In the remainder of this chapter we shall present results obtained within the peak approxi- 
mation. In this case, the nucleon form factor can be cast in the form [^| 



Q^GUQ') = ^ [47ra.(g^)J |/n|^ [ [dx] j [dy] 
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i=l 



i=8 



(118) 



where the amplitudes Ti{xj,yj) = ^i{xj)Tii{xj,yj)^i{yj) represent convolutions of Th with 
the appropriate distribution amplitudes evaluated term by term for each contributing dia- 
gram (marked by the index "i"). Inputting 

^N{Xi) = 0as(Xi)[(5o + B2- 5^3 " S^g) + (5i + B,){xi - X3) + {SB^ + 7Bs + 7B,)X2 

4 14 4 14 

+ (453 + UB,)xiXs + (SBs - -E4 + —B,)xl + (SB-, + -B^ + —B,)xl] (119) 

and integrating over variables Xi and yi, the proton and neutron form factors acquire the 
following form 



1 
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(120) 
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121) 



Analytic expressions for and 

|; such up to order M = 3 



where and J" are functions of the coefficients B„ 
up to order M = 2 {n = 0, 1, 2, . . . , 5) were published in 
(n = 0, 1, 2, . . . , 9) are given in Appendix B. 

The results for the main models in use are shown in Fig. |18] in comparison with existing 
data (for the proton) and ||123|| (for the neutron). Note that the Q^-evolution of the coef- 
ficients Bn has been neglected and that the average value as{Q'^) = [<ys{Q'^Xiyi)as{Q'^Xjyj)Y^'^ 
has been used to account for the different virtualities of the involved gluon propagators. For 
instance, for the heterotic amplitude we use ^^(Q^) = [as(<5^ x 0.427)a;s(Q^ x 0.178)]^/^. 
Here and below the values Aqcd = 180 MeV and = 5.0 x 10-^ GeV^ fS^H are 
taken. 

Since there are no direct data on beyond, say, 5 GeV^, it is virtually impossible 
to extract experimental values of the neutron form factors in a model-independent way. 
Therefore, we show in Fig. ^ the theoretical predictions along with GJ^-data, extracted 
under the proviso of two extreme, albeit reasonable assumptions: (1) assuming that the 
electric neutron form factor is at least of the order of the magnetic one (open circles), and 
(2) assuming GJ^ = (black dots). 

The main conclusion to be inferred from these results is twofold: (1) The heterotic model 
yields the best relative agreement with the experimental data, though the prediction for the 
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Figure 18. Comparison with available data of the magnetic form factor of the proton and the neutron 
for the heterotic and other model distribution amplitudes, labeled by the acronyms of the corresponding 
authors. 
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neutron form factor is too small to account for the observed form factor. (2) The COZ°p* 
amplitude - designed to yield optimum agreement with the COZ sum rules - fails in both 
cases to describe the experimental data. 

Another place to test these results is in the data for the elastic cross sections cXp and ct„. 
For small scattering angles, where the terms oc tan'^{6/2) can be neglected, there are two 
main possibilities for the ratio Cn/crp. If the Dirac form factor is zero, or small compared 
to the Pauli form factor F2 [ |125| |, then cr„ should be due only to the higher-order term F^. 
At large the ratio would become 

(modulo logarithmic corrections due to anomalous dimensions) and would decrease with 
increasing due to the extra power of of the Pauli form factor. Alternatively, if 
is comparable to F2, then cr„ would eventually be due to at large Q^. Then the ratio 
(Jn/up would be given by some constant determined by the nucleon wave functions 

^ ^ (91. 
Op ^Cf - 

In the above expressions, the wave-function characteristics are parameterized by the (dimen- 
sionful) coefficients C,, which are functions of the expansion coefficients i?„ and the "proton 
decay constant" |/n| = (5.0 ±0.3) x lO'^ GeV^ Oi^. 

From this analysis it follows that in the domain where the helicity-flipping parts of 
the form factors can be ignored, Onjop should be within the range 0.238 and 0.01. If the 
results of the combined KS/COZ sum rules are taken seriously, then this range is somewhat 
shifted: 0.232 -h 0.005. Comparing with available data |123|, we see that the measured 



(123) 



value of (Jn/op enters the estimated range already at ~ 8 GeV^/c^ (see Fig. |T9|). Here 
two concluding remarks are in order: (1) The present accuracy of QCD sum rules (modulo 
technicalities) is sufficient to provide bounds on Onjo-p within the observed region. (2) The 
available data in the range ~ (8-^10) GeV^/c^ are well below the calculated upper bound 
and still decreasing. This indicates that distribution amplitudes which give IGJJjI/G^^ ~ 0.5 
may be in contradiction to experiment because they yield a Dirac form factor F° which 
starts to overestimate the data already at ~ 8 GeV^/c^. In particular, there is no 
room for contributions due to the Pauli form factor. On the contrary, models which give 
a small value of |GmI/^m can explain the data only under the assumption that in this Q"^ 
region the Pauli contribution is still dominant. For instance, the heterotic model requires 
at = 10 GeV /c^ a Pauli contribution to the neutron form factor approximately three 
times as large as that from the Dirac form factor. The principal result from this discussion 
is that in the intermediate domain, an/cTp should be within the range 0.238 and 0.001. 

We close this section by presenting results for the axial- vector and isoscalar axial- vector 
form factors of the nucleon. These form factors were analyzed with perturbative QCD by 
Carlson and Poor in [ |12(j| ] and [ |127| ]. The axial- vector currents in question are: 



d-f^-f5U , Ajf) = ^ (u7^75U - rf7^75C?) , (124) 
= I (w7/.75W + d-f^l5d) . (125) 
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Figure 19. Bounds on the elastic cross sections (Jn/cTp obtained from QCD sum rules (shaded area). 



The data are taken from |123] 



The first two currents are related by an isospin rotation. Then the axial-vector form factor 
qa and the pseudoscalar form factor gp are defined by the matrix element with nucleon 
states as follows 



{n{p\>:)\A^p{p){PA)) = uy{p') gAmi'r + 9p{Q')q'r mp) 



(126) 



where q = P' — P, and = —q^. The axial- vector form factor can be isolated by taking 
the A'^ component and working in a frame with g+ = 0. Then 



{n\A+\p) = 2P+gA, 



(127) 



and 



y(3) 



{p\A+^'>\p) = 2P+Gf 
{p\A+^^^\p) = 2P+Gf . 

The first two form factors are related by isospin invariance to give 



.(3) 



1 



(128) 



(129) 



while the last one is the non-leading isoscalar form factor. Both form factors can be cast in 
convolution form to read 



9a{Q' 



\dx\ f\dy]^yi, Q')TH5{xi, y^, Q'Mxi, 
Jo 



(130) 



1 rl 

[dx] / 
Jo 



131) 
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Assuming, as previously, constant arguments of the strong coupling constants and expanding 
the nucleon distribution amplitudes in terms of their eigenfunctions within the basis of Appell 
polynomials, the integrals over the fractional momenta can be carried out to arrive at the 
following expressions 



47r\2 
27/ 



(132) 



Q'Gf{Q') = \a,m]'ls , (133) 



27/ 

where I a and Is are functions of the non-perturbative expansion coefficients B^. Carlson 
and Poor truncate the eigenfunctions decomposition at leading order M = 2. A more general 
expression for qa which includes the next-to-leading order eigenfunctions (M = 3) is given 
in Appendix B. 

Let us now present some results for the heterotic model. The calculation for the nucleon 
axial form factor Qa^Q'^) according to |l26l yields at ~ lo GeV^ Q^Qa^Q"^) = 0.90 GeV^ 
for Aqcd = 100 MeV, and Q^c/a(Q^) = 1.44 GeV^ for Aqcd = 180 MeV. These results 
compare fairly well with the value Q^QAiQ"^) ~ 1-5 GeV^ extrapolated from the data ||128 



Also the ratio QAiQ"^) IGy^iQ"^) ~ 1.19, in the region where the calculations can still be 
trusted, is consistent with the (extrapolated) experimental value Qa^Q"^) IG\a{Q'^) ~ 1-35. As 
for the isoscalar nucleon form factor we find at ~ lo GeV^ Q^Gf\Q'^) = 0.83 GeV^ 
for Aqcd = 100 MeV and Q^G^a\Q^) = 1-34 GeV^ for Aqcd = 180 MeV. Assuming isospin 
invariance, we combine these results with those for gA to obtain Ga^/Ga^ ~ 1.85, where 
Ga is the isovector axial- vector nucleon form factor. If a dipole form 



with G[^\o) = 0.38 from SU{Q), is used to describe the dependence of G^"* ||127|| , then, 
in the high-Q^ region, the heterotic model yields Mas = (1-15 — 1.22) GeV for Aqcd = 
100 MeV and Mas = (1-27 - 1.37) GeV for Aqcd = 180 MeV. These values comply with 
the experimental value Ma = (1.032 ± 0.036) GeV. More systematic results are compiled in 
Table |. 

B. NUCLEON-DELTA TRANSITION FORM FACTOR 

The electromagnetic N — A transition is another exclusive process which offers the pos- 
sibility to test the quality of the model distribution amplitudes for the nucleon and the A 
resonance, obtained in the previous chapters. This process involves only interactions medi- 
ated by vector bosons so that, neglecting quark masses, quark helicity is conserved. Hence, 
analogously to the elastic electron-nucleon scattering, one can define a helicity-conserving 
form factor Gj^"^^ which according to perturbative QCD should fall off like l/Q^ when 
becomes large. On the experimental side, there are no exclusive data for above 3 GeV^, 
though there are inclusive results, obtained in SLAG experiment E133 spanning the range 
2.5-10.0 GeV^ |8|. 
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Table 8. Values of the axial-vector form factor qa and the isoscalar axial-vector form factor 
G^'' at = lOGeV^/c^ for Aqcd = 180 MeV. The estimated empirical dipole masses Ma, Mas 



are also shown for some selected model distribution amplitudes for the nucleon. 



Model 


0'^ a A 

^ tf/i 


Ma 


^ ^ A 


Ma^ 


1 ^ A 


QaIG\, 


Het 


1.4374 


1.0347 


1.3379 


1.3700 


1.862 


1.195 




4.1185 


1.3462 


1.6104 


1.4348 


0.782 


1.531 


Qgopt 


2.2419 


1.2236 


2.4376 


1.5915 


2.175 


1.111 


cz 


3.5349 


1.2957 


1.2152 


1.3373 


0.688 


1.529 


coz 


4.2747 


1.3588 


1.6404 


1.4414 


0.768 


1.536 


KS 


6.1508 


1.4882 


2.8808 


1.6593 


0.937 


1.489 


GS 


2.2130 


1.1526 


2.2436 


1.5588 


2.028 


1.119 



The form factor Gf^"*^ was calculated within perturbative QCD by Carlson ||129|| . Start- 
ing from the convolution form 



[dx] / [dy]^Aiyi,Qy)Tnixi,yi,Q)^Nixi,Q^) , 



(135) 



10 Jo 

the result after integrating over longitudinal momenta can be expressed in the following form 
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1 2 



V2 



/(A) 



(136) 



where it is assumed that the argument of as depends on only, and I^^^ is a function 
of the expansion coefficients which project the A distribution amplitude on to the 
eigenfunctions of the interaction kernel: 



(A) j^{N) 
j 



(137) 



The hard part of the process derives from Th, encountered in the calculation of the nucleon 
form factor, so that $a can be expanded over the same eigenfunctions basis, i.e., the Appell 
polynomials. The difference here is that, as stated in Subsection [111 l*] , only symmetric terms 
under Xi ^ x^, contribute. The coefficients Eij in Eq. ( |1371 ) are calculable with perturbative 
QCD. Their values for i, j = 0, 1, . . . , 5 are tabulated in | |129| |. The non-perturbative expan- 
sion coefficients -B^^-* can be taken from Table ^ those for B^'^ from Table |^ Several options 
are possible, depending on the favored choice among the various distribution amplitudes. In 
Fig. |2^ we show predictions (solid lines) for 7]9A+, calculated with the heterotic distribu- 
tion amplitude for the nucleon and various A distribution amplitudes, in comparison with 
existing experimental data (see |^8| , |130|| ). The short horizontal lines on the right margin are 
predictions for the absolute value of at = 15 GeV^ from other nucleon distribution 
amplitudes, listed in Table ^ In all considered cases the CP value |/a| = 11-5 x 10~^ GeV^ 
has been used, which is within the spread of the FZOZ estimate. We emphasize that the 
sign of predicted by CZ PU 



COZ [^, and GS comes out negative for all A distri- 
bution amplitudes discussed here, accept for the amplitude FZOZ°p* which yields for COZ°p* 
a small positive form factor by construction. In contrast, convolution of FZOZ°p* with the 
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Figure 20. Transition form factor 'jpA^ , calculated with the heterotic distribution amplitude for the 
nucleon |^ , ^ and various distribution amplitudes in comparison with existing experimental data (see 
[130|). The dashed-dotted line shows a calculation |4^] which takes into account a dynamical gluon mass to 
saturate at low Q^. The dashed line illustrates the effect of evolution of the coefficients B^. Predictions 
from previous nucleon distribution amplitudes are indicated on the right margin. The open circles denote 
the data from Stoler's analysis. Note that for all curves the value |/a| = 11.5 X lO"^ Q^y^ from |4| 
has been used. 



heterotic nucleon distribution amplitude yields a form factor which overshoots the data by 
orders of magnitude. Only the heterotic nucleon distribution amplitude and the KS ||3^ one 
yield a positive sign for and only the heterotic A distribution amplitude comes close to 
the data. 

This agreement can be significantly improved. In order to account for (unknown) con- 
finement effects at low Q^, we saturate ctg in the one- loop approximation by introducing an 
effective gluon mass ||121|| : as(Q^) ^ cts(Q^ + 4mg(Q^)). In contrast to other approaches of 
this type [[12CI|J119|| , we use a dynamical, i.e., scale-dependent gluon mass derived by Cornwall 



m„ 




-1 -12/11 



A2 



;i38) 



Due to the positivity of the anomalous dimension of the mass operator, this gluon mass van- 
ishes asymptotically. This soft behavior at short distances leaves the validity of form-factor 
evolution at large momentum transfer virtually unaffected (as it should, if the dynamical 
mass generation is to be consistent with the renormalization group). In the fit represented 
by the dashed-dotted line (Fig. [2^) we use = 380 MeV, which matches Cornwall's con- 
sistency relation, mentioned previously. Referring to the same figure, we see that including 
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Figure 21. Transition form factor 7pA^, calculated with the FZOZ (top figure) and the CP distribution 
amplitude (lower figure) for the A resonance, and various nucleon distribution amplitudes in comparison 
with existing experimental data (see |l3Ct] ). The dashed-dotted line shows the calculation with the COZ 
distribution amplitude; the dashed line that one with the COZ°P* distribution amplitude. Predictions from 
previous nucleon distribution amplitudes are again indicated on the right margin. The open circles denote 
the data from Stoler's [Q analysis. Note that for all curves the value |/a| = 11.5 x 10"^ GeV^ from |4| 
has been used. 



65 



the perturbative evolution of the coefficients S^, it is sufficient to provide a good fit to 
the data above ~ 3 GeV^/c^ (dashed hue) (see, the recent data analysis in [46,|7||). At 
lower values, additional non-perturbative parameters have to be introduced, e.g., parton 
transverse momenta, effective parton masses, higher-twist contributions, quark clustering, 
etc. to account for the limitations of the leading-order formalism. In a recent analysis ||131| 
Belyaev and Radyushkin calculated on the basis of local quark-hadron duality the purely 
non-perturbative soft contribution to the 7*p — > A form factor and showed that their result 
can account for the observed behavior of the form factor at values as low as 3 GeV^. 
They also argued that this contribution remains leading up to the highest laboratory Q'^ 
values without any need for the contribution due to hard-gluon exchange (embedded in Th). 
It is clear that experimental efforts to measure beyond 10 GeV^/c^ would be extremely 
helpful to discriminate among these options. 

To complete our discussion, we show in Fig. ^ the form-factor predictions for ob- 
tained from other nucleon distribution amplitudes in convolution with the model A distribu- 
tion amplitudes of CP (top figure) and FZOZ (lower figure). One may interpret these figures 
as indicating that none of the shown amplitudes neither for the nucleon nor for the A can 
account for the observed behavior of the form factor. COZ-type distribution amplitudes for 
the nucleon exhibit strong cancellation of symmetric vs antisymmetric contributions | |132| , 
while GS-type distribution amplitudes yield unrealistically large negative contributions. 



V. CHARMONIUM DECAYS 



The last phenomenological application we consider here within the standard convolution 
scheme here deals with exclusive decays of charmonium states to pp and AA. Such decays 
are sensitive to the nucleon (A) distribution amplitude and hence may serve to discriminate 
among proposed models. Calculations of charmonium decays have been performed by many 
authors p6 i|l 3 3|,p!3^ , |l 3 5|] within the QCD convolution framework. We follow [^. We con- 
sider ffist the ^Pj states with J = 1,2. The branching ratio for the decay of the Xci state 
( J^*-^ = 1"*"^) into "pp is given bylll 



BR 



0.75 



^Pi ^ all; ln(M/A) 



167r2 


/n 


729 


M2 



Ml 



(139) 



where M ^ 2mc ~ 3 GeV and A = 0.4 GeV (the last value from [|136|| ). The non- 
perturbative content of Eq. ( |139|) is due to /n and the decay amplitude for the process 
^Pi — > pp, denoted Mi, which involves $iv. Inserting expression ( |119|) , in connection with 
Table ^, the decay amplitude Mi for each model is computed with an elaborate integration 
routine which properly takes account of contributions near singularities [ p,04|| . The results 
for different model distribution amplitudes are compiled in Table |^. 

The analogous expression to ( |139D for the Xc2 state {J^'" = 2~^~^) has the form 



15 



I wish to thank Victor Chernyak for bringing to my attention that the factor (ttcks)^ was missing 
in [^Pp3|. 
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Table 9. Exclusive charmonium decays into pp for a variety of nucleon 

The numbers in parentheses are those g 

381 



The data are taken from [ 137 | 

amplitudes termed Mq were independently verified by Bolz 



distribution amplitudes, 
iven in IM]. The decay 



Models 


M{^Pi pp) = Ml 


71 J" / .'"i 7^ — \ 71 'T 

M('^P2 pp) = M2 


M{^Si pp) = Mo 


heterotic 


99849.6 


515491.2 


13726.8 


LZ 


28olU.4 (13.60X113"^) 


246052.8 (2.87x10 j 


7545.6 (0.72x10 j 


COZ 


53625.6 (0.88x10^) 


298123.2 (3.4x10^) 


8758.8 (0.79x10^) 


coz°p* 


55137.6 


289728.0 


8499.6 


GS 


26366.4 


232632.0 


928.8 (0.7xl03) 


Qgopt 


19915.2 


230659.2 


986.4 


Qgmin 


17193.6 


210729.6 


964.8 


KS 


94723.2 (1.35x10^) 


416937.6 (4.84x10^) 


11484.0 (1.15x10^) 


rm asympt. 


20086.6 (0.2x10^) 


43099.2 (0.43x10^) 


1517.4 (0.14x10^) 


Observables 


^^[^) in % 


BR(^)in% 


r(35i ^ pp) [eV] 


heterotic 


0.22 xlO-2 


0.89 X 10-2 


138.37 


CZ 


0.017x10-2 


0.20 X 10-2 


41.81 


COZ 


0.06 xlO-2 


0.30 X 10-2 


56.34 


coz°p* 


0.066 X 10-2 


0.28 X 10-2 


53.07 


GS 


0.014x10-2 


0.18 xlO-2 


0.63 


GS°P* 


0.086 X 10-3 


0.18 xlO-2 


0.71 


Qgmin 


0.066 X 10-3 


0.15 xlO-2 


0.68 


KS 


0.198 X 10-2 


0.59 X 10-2 


96.84 


asympt. 


0.086 X 10-3 


0.01 xlO-2 


1.69 


E760 


(0.78 ± 0.10 ±0.11)xl0-2 


(0.91 ±0.08 ±0.14) X 10-2 


180 ± 16 ± 26 
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BR 



pp 



3P2 



all 



0.85(7ras 



,16_ 
729 



M2 



(140) 



which is Eq. (20) of with an obvious minor correction. The results for the branching 
ratio of this process, shown in Table have been calculated with as{mc) = 0.210 ± 0.028 



(see third paper of ||136|| ). 

Similar considerations apply also to the charmonium decay into pp of the level ^5*1 with 
J^*" = 1 . The partial width of J/ip (i.e., Xco) into pp is defined by 



TCSi^pp) = {7ia,y 



1280 \U\^ 


/n 


2437r M 


M2 



where determines the value of the ^S'l-state wave function at the origin, viz, 



;i4i) 



(142) 



where is the J/ ip polarization vector. The value of can be extracted from the leptonic 
width r(3^i ^ e+e" 



= (5.36 ± 0.29) keV ||139|| via the Van Royen-Weisskopf formula: 



e e 



(143) 



The result is |/^| = 409 MeV with mj/^ = 3096.93 MeV. (Note that M ^ Mj/^.) Then, 
using the same values of parameters as before, we obtain the results shown in Table One 
sees from this table that the agreement between the predictions of the heterotic model and 



the recent high-precision data of the E760 experiment at Fermilab [|1371| is quite good. The 
only exception is the branching ratio for the process ^Pi pp which comes out too small, 
though the heterotic amplitude yields the largest value, compared to all other models. The 
experimental result can be reproduced using as — 0.32. 

In a completely analogous way, one can calculate the corresponding exclusive decays 
of (helicity-conserving) charmonium states in AA. The branching ratios follow by re- 
placing /n by /a/v^- The results of this calculation are summarized in Table |T0[ As 
regards the branching ratio of the ^5*1 state, the heterotic distribution amplitude yields 



BR{^Si AA/^Si all) 
all considered decays as{mc 



0.30 X 10-2%, where Ttot 
= 0.210 ± 0.028 as before. 



85.5t^:8 keV is used. Note that in 



VI. MODIFIED CONVOLUTION SCHEME 

In casting the hadron form factors in convolution form (cf. Eq. ( |116|) ), we tacitly assumed 
that the fc^-dependence of the quark and gluon propagators in Th can be ignored. This is 
tantamount to factorizing the fc^-dependence into the distribution amplitudes which are the 
wave functions integrated over up to the factorization scale. Thus, in the limit Q"^ ^ oo, 
the only gluon radiative corrections remaining uncanceled are those giving rise to wave- 
function renormalization. Indeed, recalling Eqs. (p^), (^8]), one sees that asymptotically 
the most likely configurations are those in which the valence quarks share longitudinal 
momentum in a uniform way, i.e., Xi = 1/2 for the pion and Xi = 1/3 for the nucleon. 



68 



Table 10. Exclusive charmonium decays in AA for the models discussed in the text. 



Amplitudes 


CP 


FZOZ 


Heterotic 




11418.23 


21585.99 


11651.26 




30924.07 


48233.23 


26277.40 




1480.67 


1882.31 


1134.94 


Observables 




0.089 X 10-3% 


0.3196 X 10-3% 


0.093 X 10-3% 




0.100 X 10-3% 


0.106 X 10-3% 


0.072 X 10-3% 


r(35i ^ AA) 


0.439 X 10-2 keV 


0.709 X 10-2 keV 


0.258 X 10-2 keV 



When confinement sets in, a quark is not able to venture too far from the anti-quark (in 
the pion) or the other two valence quarks (in the nucleon), and this poses constraints on 
the off-shellness of the involved propagators (see Ihs of Fig. ^2|). However, in the end-point 
region, the parton transverse momenta in Th cannot be a priori ignored, since, say, for 

the pion, (^k_|_j + ^'±j^ ^ x^x'jQ'^. As a result, the transverse distance between the quark 
and the anti-quark becomes large compared to 1/Q and the corresponding gluon line is no 
more part of the hard-scattering process but should be counted to its soft part. This is in 
contrast to the perturbative QCD treatment in which the struck quark connects to the other 
valence quarks via highly off-shell gluon propagators, meaning that the transverse inter- 
quark distances are rather small, viz. of order 1/Q and that all partons share comparable 
fractions of longitudinal momentum. In other words, the mechanism of hard-gluon exchange 
becomes inapplicable and should be replaced by that of Feynman ||140|| . According to the 
latter mechanism, almost all of the hadron's momentum is carried off by a single parton, the 
others being "wee" . This picture is consistent with a configuration in which only the struck 
quark is within an impact distance l/Q of the electron while all other partons have rather 
random positions in the transverse direction, building a soft "cloud" with transverse size 
^ l/Q [ p-41|l - Once the elastic scattering has happened, rearrangements are necessary to 
change quarks and gluons into an unscathed hadron. This conversion procedure (visualized 
on the rhs of Fig. |2^) is controlled by the overlap of the initial and final state wave functions 
and cannot be computed within pQCD. 



A. MODIFIED FACTORIZATION 

The physical basis of the modified convolution scheme is to dissect the process in such a 
way as that for transverse distances large compared to 1/Q (the latter being the playground 
of the hard-scattering mechanism) but still small relative to the true confinement regime 
- characterized by 1/Aqcd ^ the hadron wave function is modified to exhibit the effect of 
Sudakov enhancements explicitly up to the transverse scale retained in Th (see Fig. p3|) . 
Going over to the transverse configuration space, the modified wave function in the axial 
gauge A'^ = reads 
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Figure 22. Mechanisms for momentum transfer during elastic scattering. The Ihs shows hard-gluon 
exchange within pQCD. The blob S containing soft gluon lines (and an analogous one with soft 
quark-anti-quark lines not shown here) spoils factorization but is power-suppressed, i.e., non- leading. The 
(rhs) shows the Fcynman mechanism using, for purposes of illustration, quark and gluon lines. The leading 
quark is denoted by a heavy line, while all other lines represent wee quarks and soft gluons. 



The Sudakov exponential factor re-sums contributions from two-particle reducible diagrams 
(giving rise to double logarithms), whereas two-particle irreducible diagrams (giving rise to 
single logarithms) are absorbed into the hard scattering amplitude Th [76]. Hence, e~'^ can 
be conceived of as being a finite IR renormalization factor to the hadron's wave function 



25| , encoding the exponentiation of the probability for no-emission of soft gluons. This 
factor renormalizes the wave function in addition to the conventional renormalization factor 
Z2 due to UV divergences, encountered before. While Z2 contains single logarithms in 
leading order, the Sudakov renormalization factor is dominated by double logarithms. The 
leading double logarithms derive from those momentum regions where soft gluons (all four- 
momentum components small) and coUinear gluons to the external quark lines overlap. 
These contributions are numerically dominated by the term 

\ A v^Aqcd ln(l/6,AQ0D) J ' 

where is one of the fractions Xi,Xi = 1 — Xi (for incoming valence quarks) or x[,x'i = 1 — x'^ 
(for outgoing valence quarks), and (3q is the first-order term of the Gell-Mann and Low func- 
tion encountered before. The single logarithm in Eq. (eq:doublelog) stems from the running 
of the coupling constant, and the double logarithm contains the exponentiated higher-order 
corrections - required by the renormalization group - rendered finite by the inherent IR- 
cutoff l/bi- This issue marks the crucial difference between the modified convolution scheme 
and previous approaches dealing with isolated quarks where such IR-cutoff parameters had 
to be introduced as external regulators. Here IR regularization is provided in situ without 
additional assumptions. For small transverse distances (or equivalently, 1/bi ^ ^iQ), glu- 
onic radiative corrections are treated as being part of Th and are excluded from the Sudakov 
form factor. Consequently, for < ^/2/biQ the Sudakov functions s{^i,bi,Q) are set equal 
to zero. On the other hand, as bi increases, decreases, reaching zero at fojAqcD = 1- 



The full expression for the Sudakov exponent was calculated by Botts and Sterman [176 . 
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Figure 23. Modified factorization of soft-gluon contributions tantamount to a finite renormalization of 
the nucleon wave function. 
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It is given b; 



1=1 



s{xi,bi,Q) + I \ — 7g(^(/i^)) 
'i/bi /i 



+ E 

1=1 



•j2 



s{x'i,bi,Q) + / —lq{g{li )) 

Jl/bi /i 



wherein the Sudakov functions s{^i,hi,Q) are given by 



AW 



Am 



1^-11- 



■qi 



hi{2hi) + 1 /n(2g;) + 1 



16/5f 



In 



32/3? 



In2(2g0-ln'(2&, 



(146) 



(147) 



Here = xi,xi or X;, (/ = 1, 2, 3) and the variables q and b are defined as follows 

qi = ln[6Q/(V2AQCD)] , (148) 



ln[l/6iAQCD]. 



The coefficients A^'^^ and Pi are 



67 1 



10 



27 



(149) 



33 - 2nf 
12 



4 



/3o, /32 



153 - 19nf 
24 ' 



(150) 



where nf is the number of quark fiavors and 7 = 0.5772 ... is the Euler-Mascheroni constant. 

The Sudakov function, s{^i,bi,Q), in Eq. (|147|) takes into account leading and next-to- 
leading order gluonic radiative corrections of the form shown in Fig. and accounts for 
renormalization group evolution from the IR-scale 1/bi to the renormalization scale firen via 
the quark anomalous dimension. The quantities bi {I = 1, 2, 3) are IR cutoff parameters, 
naturally related to, but not uniquely determined by the mutual separations of the three 
quarks ||72| . A physical perspective on the choice of the IR cutoff is provided by the following 
analogy to ordinary QED. One expects that because of the color neutrality of a hadron, its 
quark distribution cannot be resolved by gluons with a wavelength much larger than a 



^^A corrected form of this expression was derived by Bolz |142(| ; see also [ 143 |. 
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characteristic quark separation scale. Hence, long-wavelength gluons probe the color singlet 
proton and radiation is damped. On the other hand, radiative corrections with wavelengths 
between the IR cutoff and an upper limit (related to the physical momentum Q) yield to 
suppression. It is understood that still softer gluonic corrections are already taken care of 
in the hadron wave function, whereas harder gluons are considered as being part of Th. 

In the pion case, there is only one transverse scale, notably, the quark-anti-quark separa- 
tion b, so that suppression is automatically accomplished. Indeed, when it happens that one 
Sudakov function s{^,b,Q) = (meaning that the corresponding exponential is set equal to 
unity) the other (negative) Sudakov function in the exponent, s(l — ^,b,Q), diverges, thus 
providing sufficient suppression. This IR-protecting behavior due to the Sudakov form fac- 
tor, makes it possible to choose a renormalization scale which depends on the initial and final 



longitudinal momenta. It was shown in Ref. ||144||, within the collinear approximation, that 



such a subtraction point minimizes the contributions of the next-to-leading order corrections 
to the pion form factor, ensuring dominance of the leading-order perturbative contribution. 
In the standard convolution scheme, such a choice leads to singularities in the running cou- 
pling constant as, unless it is saturated by additional IR regulators; e.g., a dynamical gluon 
mass. Here the modification of as becomes superfluous because the Sudakov effect inhibits 
soft-gluon emission, thus effectively screening the ag-singularities. Physically, the Sudakov 
suppression enhances the dominance of quark configurations with a small color dipole mo- 
ment and this enhancement becomes more pronounced as increases. Concerning the 
nucleon, the situation is more complicated because several transverse scales are involved 
and therefore a careful IR regularization is required. Now the Sudakov function comprises 
six terms (^/ = xi or ^/ = x'l for / = 1,2,3), each depending on its own transverse scale 
and with corresponding QCD-evolution contributions, driven by the anomalous dimensions 
associated with quark self energy. Recall that the IR cutoff to be used is the factoriza- 
tion scale, below which operator product expansion becomes a poor approximation. It is 
also understood that genuine non-perturbative momenta below the IR cutoff are implicitly 
accounted for in the nucleon wave function. 

Different choices of the IR cutoff have been used in the literature: Thus, Li chooses 
bi = bi (this choice is termed hereafter the "L" prescription), whereas Hyer |p.45|| in his 
analysis of the proton-antiproton annihilation into two photons and of the time-like proton 
form factor as well as Sotiropoulos and Sterman ||146|| take &i = &2, &2 = &3 = &3 (this 



choice is denoted the "H-SS" prescription). Still another possibility, proposed in [^], is 
to use for reasons that will be explained below as IR cutoff the maximum of the three 
inter-quark separations, i.e., to set 

b = max{6i, 62, ^3} = &i = &2 = &3 • (151) 

This choice, designated by "MAX" , is analogous to that in the meson case, wherein the 
quark-anti-quark distance naturally provides a secure IR cutoff. The specific features of each 
particular cutoff choice will be discussed in detail later. 

The integrals in Eq. ( |146| ) arise from the application of the renormalization group equa- 
tion. The evolution from one scale value to another is governed by the anomalous dimensions 
of the involved operators. The integrals combine the effects of the application of the renor- 
malization group equation on the wave functions and on the hard scattering amplitude. 
The range of validity of Eq. (|147|) for the Sudakov functions is limited to not too small bi 
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Figure 24. The exponential of the Sudakov function s(^;, bi, Q) vs ^; and ^/Aqcd for Q 
In the hatched area the Sudakov function is set equal to zero according to Li's requirement |79 



30A, 



QCD- 



values. Whenever 1/6; is large relative to the hard (gluon) scale $,iQ, the gluonic corrections 
are to be considered as higher-order corrections to Th and hence are not included in the 
Sudakov factor but are absorbed into Th (hierarchy of scales). For that reason, Li sets 
any Sudakov function s{^i,bi,Q) equal to zero whenever < \/2/{Qbi). Moreover, Li holds 
the Sudakov factor e~^^ equal to unity whenever it exceeds this value, which is the case 
in the small 6/-region. Actually, the full expression Eq. ( 146 ) shows in this region a small 
enhancement resulting from the interplay of the next-to-leading logarithmic contributions 
to the Sudakov exponents and the integrals over the anomalous dimensions. 

The IR cutoffs 1/6/ in the Sudakov exponents mark the interface between the purely non- 
perturbative soft momenta, which are implicitly accounted for in the proton wave function, 
and the contributions from soft gluons, incorporated in a perturbative way in the Sudakov 
factors. Obviously, the IR cutoff serves at the same time as the gliding factorization scale 
/ip to be used in the evolution of the wave function. For that reason, Li |79 
Sotiropoulos and Sterman ||146|| take /ip = min{l/6/}. The "MAX' 
adopted in [^, naturally complies with the choice of the evolution scale proposed in |l79| , |146|| 



as well as 
prescription ( |151[ ), 



In Fig. we display the exponential of the Sudakov function exp[—s{C,i,bi,Q)] for Q 
30 Aqcd by imposing Li's requirement W^: s{^i,bi,Q) = whenever < ^/2/Qbi. 



B. MODIFIED NUCLEON FORM FACTOR 

The crucial element to incorporate radiative corrections in the calculation of form fac- 
tors within the modified convolution scheme is to retain the explicit fc^-dependence in the 
convolution of the wave functions with the hard scattering amplitude. This new convolution 
formula can formally be derived by using the methods described in detail by Botts and 
Sterman |7^. According to Li |79], the convolution formula for the proton form factor can 
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be written in the form 



G'm(Q^) = ^ / / [d^k±][£k'j^]^Tn^{x,x',k±,k'^,Q,iJ,)Yj{x,x',k±,k'^,lxp) . 

6 JO J 

(152) 

One recognizes that this formula appears as an intermediate step in the derivation of the 



standard hard-scattering formula for the proton form factor Note, however, that the 
notation adopted in the present exposition is slightly different from that used by Li. Making 
use of the symmetry properties of the proton wave function under permutation, discussed 
previously, the contributions from the diagrams involved in the calculation of the proton 
form factor at Born level can be arranged into two reduced hard-scattering amplitudes of 
the form 

r =2 (4^^s(/i))^ 

3 ^ [(1 - a;i)(l - x[)Q^ + (k^i - k'^,)^] [x^x'.Q^ + (k^s - k'^,)^] ' ^ ^ 
T =2 (4vr«s(/i))^ . . 

3 ^ [xix;g2 + (k^^ _ kYj2] [a^^x'^g^ + (k^^ - ki^)] ' ^ ' 

In the hard scattering amplitudes only the fc_L-dependence of the gluon propagators is in- 
cluded, whereas that of the quark propagators has been neglected. There is, in principle, no 
problem to include in the calculation all fc^-dependence of Th. In the case of the pion form 
factor this has been exphcitly demonstrated by Li |79[. He found an additional suppression 



of the final result of about 10%. However, in the proton case, such terms would lead to an 
1 1-dimensional integration which cannot be carried out to sufficient accuracy with present 
day computers.0 

The functions Yj in Eq. ( |152| ) are short-hand notations for linear combinations of products 
of the initial- and final-state wave functions "ifijk ,^X'j'k'^ weighted by a;i-dependent factors 
arising from the fermion propagators, namely: 



Yl = ^{4^^23^123 + 4^^32^132 + ^231^231 + ^321^321 

+2^*'31^132 + 2^t:i2^231 + 2^^21^123 + 2^r23^32l} , (155) 

% = ^ {3^r32^132 - ^^31^231 " ^231^132 " ^r32^23l} 

^ {4^*'2,^32l + ^^23^123 + 2^*'21^123 + 2^^23^321} (156) 



X3X1 



^^Explorative studies by using, for instance, the A;_l -dependence of only a certain subset of quark 
propagators in addition to those of the gluon propagators (leading to a 9-dimensional integration) 
show that the inclusion of such terms yields an additional suppression of the perturbative result. 
Since we are mainly interested in estimating rather the maximum perturbative contribution, we 
dispense with such terms. 
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for the proton [|78[ and by 



-'2 



1 

=-=^{-2^^23^123 - 2^^32^132 + ^231^231 + ^321^321 
-^231^132 - ^^32^231 - ^321^123 " ^m^32l} , 

^ {^231^231 + ^231^132 + ^t32^23l} 

+ J=r {2^321^321 - ^r23^123 + ^321^123 + ^^23^321} 



(157) 



XsXi 



(158) 



for the neutron |S^. The subscripts on refer to the order of momentum arguments, for 
example \l/i23(x, k^) = \l/(xi, k_Li; X2, k_L2; X3, k_L3). Note that, in general, the wave function 
depends on the factorization scale fip. We make the following convenient ansatz for the wave 
function: 



^123(3;, kj 



1 



/N(/iF)$iv(a;, /ip = l/bi) n{x, kj 



(159) 



The A;_L-dependence of the wave function is contained in the function Q which is normal- 
ized according to 



[d^k_L]l^i23(a;, kj 



1 . 



(160) 



Due to Eqs. (|9T|) and ( |160|) , /n is the value of the wave function at the origin of the configura- 
tion space. Recall that in contrast to the pion decay constant f^^, which has zero anomalous 
dimension, /n exhibits evolution behavior driven by the leading anomalous dimension 70 
according to 



/n(/^f) = /N(/io) ( 



2/3/3o 



;i6i) 



At the starting point of evolution it has the value obtained from QCD sum rules ||38| , |37| : 
/^(/io) = (5.0 ± 0.3) X 10-3 GeV2. 

In Eq. (|159|) \E' represents the soft part of the proton wave function that results by 
removing the perturbative part and absorbing it into the hard-scattering amplitude Th. 
The perturbative tail of the full wave function behaves as l/k\ for large fc_L, whereas the 
soft part vanishes as l/k\ or faster. The non-perturbative or intrinsic /c^-dependence of 
the soft wave function, being related to confinement, is parameterized in our analysis as a 
simple Gaussian distribution according to 



123 



[x, kj 



'167r 



2\2 



X1X2X3 



exp 



i=l 



(162) 



This parameterization of the intrinsic fc^-dependence of the wave function, which is due 



to Brodsky, Huang, and Lepage ||147|| , seems to be more favorable than the standard form of 
factorizing x- and fc^-dependences. At least for the case of the pion wave function, this has 
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recently been effected by Zhitnitslcy on the basis of QCD sum rules. He found that a 



factorizing wave function is in conflict with some general theoretical constraints with which 
any reasonable wave function should comply. Zhitnitsky's QCD sum-rule analysis of the 
pion wave function seems to indicate that the fc^-distribution may also show a double-hump 
structure (like the distribution over longitudinal momenta in the CZ pion distribution am- 
plitude), which means that small and large values of k± are favored relative to intermediate 
values. It is likely that the proton wave function may exhibit a similar behavior, though 
this kind of analysis has yet to be done. 

In Eq. ( |162| ) the parameter a controls the root mean square transverse momentum 
(r.m.s.), {k^y^, and the r.m.s. transverse radius of the proton valence Fock state. From 
the known charge radius of the proton, we expect the r.m.s. transverse momentum to be 
larger than about 250 MeV. The actual value of (^i)^''^ may be much larger than 250 MeV, 
even as large as 600 MeV. Indeed, Sotiropoulos and Sterman [|146|| show that application of 
the modified convolution scheme to proton-proton elastic scattering leads to an approximate 
t~^-behavior of the differential cross section at moderate \t\. The behavior da/dt ~ t^^^, 
predicted by dimensional counting, appears only at very large |t|. At precisely which value of 
\t\ the transition from the to the t~^^ behavior occurs, depends on the transverse size of 
the valence Fock state of the proton. Since the ISR [ p.48| | and the FNAL [ p.49[ | data are rather 
compatible with a t~^-behavior of the differential cross section, Sotiropoulos and Sterman 
conclude that the transverse size of the proton is small, perhaps < 0.3 fm. Correspondingly, 
the r.m.s. transverse momentum is then larger than 600 MeV. It is worth noting that such a 



large value is supported by the findings of the EMC group ||150|| in a study of the transverse 
momentum distribution in semi-inclusive deep inelastic fip scattering. A phenomenologi- 
cally successful approach to the standard convolution scheme, in which baryons are viewed 
as bound states of a quark and an effective diquark particle, also uses a value of this size 
for pi| ^ |i5^ ^ p3^ . There is a second constraint on the wave function, namely the 



probability for finding three valence quarks in the proton: 

P ,4 f'u i 2('l'i23(x))' + 'l'l32(x)'l'23l(a^) 

"sq = (tto) / [dx\ < 1 . (163) 

3 JO X1X2X3 



Similarly to Sotiropoulos and Sterman |p.46||, we write the valence quark component of the 



proton state with positive helicity in the form 

|P,+) = Ijdx] J [d\^]{ X1>123^™ + ^213^ 



aia2a3 



- (^132 + ^23l)A^+T_"^}eaia,a3 , (164) 



where we assume the proton to be moving rapidly in the 3-direction. Hence, the ratio of 
transverse to longitudinal momenta of the quarks is small. The 3-quark state with helicities 
Ai, A2, A3 and colors ai, 02, 03 is given by 

^\\Z% = l^ai; Xi, k_Li, Xl)\Ua2] X2, k_L2, A2) l^a; X3, k_L3, A3) . (165) 

y/XiX2X3 

Since the orbital angular momentum is assumed to be zero, the proton helicity is the sum 
of the quark helicities. The quark states are then normalized as follows 
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X. 



k'^„ KlQa^x,, \) = 2x,{27r)%,,A'xS{< - x,)(5(k'^, - k^J . (166) 



In the numerical analysis of [[78| , p0| , to be presented below, we make use of two different 
values of the r.m.s. transverse momentum: (1) One option is to use that value which is 
obtained by requiring P^^ = 1 for a given wave function. [This corresponds to the minimum 
value of the r.m.s. transverse momentum.] (2) As another option for the r.m.s. transverse 
momentum, we consider the rather large value of 600 MeV. In the latter case, the probability 
for the valence quark Fock state is much smaller than unity and depends on the structure 
of the particular wave function. 



In order to include the Sudakov corrections, it is advantageous to re-express Eq. ( |152| ) in 



terms of the variables b_Lj, which are canonically conjugate to k_Lj and span the transverse 
configuration space. Then 

GuiQ') = idx][dx'] j -^-^E ^;(^'^''b±,Q,/^)^;(a;,x',b^,/XF)e-^^ , (167) 

where the Fourier transform of a function /(k^) = /(k^i, kx2) is defined by 

/(b_L) = -^^^ j (i2k_Lic/2k_L2exp{-ib^i- k^i - ib_L2- k_L2}/(k_L) . (168) 

Since the hard scattering amplitudes depend only on the differences of initial- and final- 
state transverse momenta, there are only two independent Fourier-conjugate vectors b_Li (= 
and b_L2 (= b^g). They are, respectively, the transverse separation vectors between 
quarks 1 and 3 and between quarks 2 and 3. Accordingly, the transverse separation of quark 
1 and quark 2 is given by 

b±3 = b±2 - bii . (169) 



[Note that Sotiropoulos and Sterman ||146|| define the transverse separations in a cyclic way 
which results in the interchange b_|_x < — > ~b_|_2, as compared to the definition used here.] 

The fact that there are only two independent transverse separation vectors is a conse- 
quence of the approximation made in the treatment of the hard scattering amplitudes (given 
by Eqs. ( |153D and ( |154| )) which disregards the /c^-dependence of the quark propagators. This 



approximation is justified by the enormous technical simplification it entails, given that the 
thereby introduced errors are very small. Then by virtue of rotational invariance of the sys- 



tem with respect to the longitudinal axis, the form factor (|167|) can be expressed in terms 
of a seven-dimensional integral instead of an eleven-dimensional one. Physically, the rela- 
tions h^i = h±[, h±2 = bx2 mean that the physical probe (i.e., the photon) mediates only 
such transitions from the initial- to the final-proton state, which have the same transverse 
configurations of the quarks. 

The Fourier-transformed hard scattering amplitudes appearing in Eq. ( p.67|) read 



Ti = -CFa,{tn)as{tu)Ko ( {1 - Xi){l ~ x[)Qh ) ( \/x2x'2Qb2 ) , (170) 



T2 = -CFasit2l)as{t22)Ko[^JxlX[Qbl] Ko[^Jx2X2Qb2] , (171) 
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where Kq is the modified Bessel function of order 0, and 6/ denotes the length of the cor- 
responding vector. We have chosen the renormahzation scale in such a way that each hard 
gluon carries its own individual momentum scale tji which in turn appears as the argument 
of the corresponding ctg. The tji are defined as the maximum scale of either the longitudinal 
momentum or the inverse transverse separation, associated with each of the gluons: 



tii= max 
t2i= max 



^{l-x^){l-x[)Q, 1/W 



\Jxix\Q, l/hi 
ti2= t22 = max 



X2X'2Q, 1/&2 



(172) 



One may think of other choices. However, they are not expected to lead to significantly 
different predictions for the form factor |7^ . 

The quantities Yj contain the same combinations of initial- and final-state wave func- 
tions as those in Eq. ( |155| ) and Eq. (|156|) , the only difference being that now the products 
^*/ /^/\E'ijfc are replaced by corresponding products of Fourier-transformed wave functions: 



^'*yfc/(a;',b_L,/iF)^ijfc(a:, b_L,AiF)- Using Eq. 
the wave function reads 



and Eq. 



the Fourier transform of 



^123(2;, b_L,/XF) 



1 



h{l^F)^123ix, flF)Qi23{x, bj 



where the Fourier-transform of the fc^-dependent part is given by 

1 



fiio.sfx, bj 



(47r) exp 



4a2 



xiXsbj + X2X362 + xiX2bl 



(173) 



(174) 



The exponentials e"'^^ in Eq. ( |167| ) are the Sudakov factors, discussed previously, which 
encapsulate the effects of gluonic radiative corrections. This makes it apparent that Eq. (|167|) 
is not simply the Fourier transform of Eq. (|15!2|) but an expression comprising additional 



physical input. Thus, with hindsight, Eq. ( |167| ) may be termed the "modified convolution 
formula" in which hard-gluon re-scattering can still be isolated (factored out). 



1. Screening of ag singularities 

Using as argument in ctg a renormahzation scale, which is independent of the longitu- 
dinal momentum fractions, may induce large contributions from higher orders in the end- 
point region. Indeed, for the pion form factor this has been explicitly shown, at least at 
next-to-leading order ||144|| . Surely such large higher-order contributions would render the 



leading-order calculation useless. For this reason, one may be tempted to use as a more 



appropriate choice the renormahzation scale yx2x'2Q, for such a scale would eliminate the 
large logarithms arising from the higher-order contributions. Unfortunately, this is achieved 
at the expense that as becomes singular in the endpoint regions and the convolution form 
of the form factors diverges. To render the form factors finite, additional external param- 



eters, like an effective gluon mass ||121|| or a cutoff prescription have to be introduced, as 
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we discussed before. One of the crucial advantages of the modified convolution scheme, is 
that there is no need for external regulators because the Sudakov factor may suppress the 
singularities of the "bare" (one-loop) as inherently. Indeed, in the pion case, it was shown 



77| that the transverse quark-anti-quark separation is tantamount to an IR regulator which 
suffices to suppress all singularities arising from the soft region. 

Concerning the proton form factor, we shall effect in the following that the screening 
of the ag-singularities by the Sudakov factor depends sensitively on the choice of the IR 
regularization prescription in transverse configuration space. 

Ultimately, the suppression of the as-singularities relies on the fact that whenever one of 
the Qfs couplings tends to infinity (owing to the limit tji —* Aqcd), the Sudakov factor e~'^^' 



rapidly decreases to zero. As it can be observed from Fig. this is not the case in the region 
determined by < \/2Aqcd/Q simultaneously 6/Aqcd — ^ 1, where exp[—s{^i,bi,Q)] is 
fixed to unity. In the pion case thjs does not matter, since the other exp[— s(l — ^, b, Q)] — *• 
faster than any power of ln[l/(6AQCD)] and, consequently, the Sudakov factor drops to 
zero ensuring IR protection. The treatment of the proton form factor is technically more 
subtle. In that case, e"'^^ does not necessarily vanish fast enough to enforce the suppression 
of the Ofs-singularities. This can be illustrated by the following string-like configurations: 
xi < a/2Aqcd/Q and 6iAqcd 1 and 62AQCD — ^sAqcd — 1/2. Using the "L" prescription, 
it is obvious that s{xi, bi,Q) =0 and the other two Sudakov functions are finite. Therefore, 
there is no suppression of the cts-singularities owing to the limit 6iAqcd 1- Sufficient 
suppression of the ag-singularities is provided only if all three bi are coerced to be equal. 
Allowing for the three bi to be different, the Sudakov factor provides suppression only through 
the contributions of the anomalous dimensions which are only logarithmic, and hence cannot 
provide sufficient suppression. Similar observations hold for the "H-SS" prescription. Since 
the "L" and "H-SS" prescriptions allow for different bi values in the Sudakov functions, the 
integrand in ( |167|) has singularities behaving as 

^ ln-« I _i — I (175) 
V6/Aqcd/ 

for 6;Aqcd — 1 and xi hold fixed. The maximum degree of divergence is given by 

^=^(t + 27max - 2) + 1 , (176) 



Po V3 

where the first term 4/3 comes from the evolution of /n, Eq. ( [L61| ), and the constant 7max is 
related to the anomalous dimension driving the evolution behavior of the proton distribution 
amplitude (see Eq. (|48|) and Table |lT]). Here 7max is the maximum value of the {%} within 



a given polynomial order of the eigenfunctions expansion of the distribution amplitude. We 
reiterate that the 7^ are positive fractional numbers increasing with n. Thus, the singular 
behavior of the integrand becomes worse as the expansion in terms of eigenfunctions extends 
to higher and higher orders. The term —2 in Eq. (|176|) stems from the integrations over 
the anomalous dimensions in the Sudakov factor e"'^^ (see Eq. ( |146D ). Finally, the term 1 
originates from that as{tjk) which becomes singular in Eq. ( |167|) (cf. Eq. (|172|) ). Which one 



of the tts couplings becomes actually singular depends on the prescription imposed on the 
IR cutoff parameters 6/. The integral in Eq. (|167|) does not exist if 7max > |. As Table ^ 
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Table 11. Expansion coefficients for selected nucleon distribution amplitudes, taken from |38] 
and [ |5H^ , and used in the modified convolution scheme. The notation of is adopted. The 
{jn} are related to those given in Table ^ as follows: 7^^^^*^ ^ = -^{^ + In), where | = 70 in this 
table. The associated r.m.s transverse momentum and the oscillator parameter for each model 
wave function, normalized via = 1, are shown. 



n 


7n 


B„(COZ"P) 


B„(COZ) 


S,,(COZ°Pt) 


Sn(Het) 


B„(low) 


1 


20/9 


3.2185 


3.6750 


3.5268 


3.4437 


4.1547 


2 


24/9 


1.4562 


1.4840 


1.4000 


1.5710 


1.4000 


3 


32/9 


2.8300 


2.8980 


2.8736 


4.5937 


3.3756 


4 


40/9 


-17.3400 


-6.6150 


-4.5227 


29.3125 


26.1305 


5 


42/9 


0.4700 


1.0260 


0.8002 


-0.1250 


-0.5855 


a [GeV- 


[MeV] 


271 
0.9893 


271 
0.9939 


267 
1.0069 


317 
0.8537 


299 
0.9217 



reveals, this happens already for proton distribution amplitudes which include (Appell) 
polynomials of order 1, i.e., for all distribution amplitudes except for the asymptotic one. 
Thus application of the "L" and "H-SS" prescriptions to choose the IR cutoff parameters 
bi in the proton form factor entails the modified convolution scheme to be invalid. In view 
of these results, Li's numerical analysis of the proton form factor |^ is seriously flawed. 
Nevertheless, his general approach is a decisive step towards a deeper understanding of the 
electromagnetic baryon form factors. 

A simple recipe to bypass the singular behavior of the integrand is to ignore completely 
the evolution of the distribution amplitude or to "freeze" it at any (arbitrary) value larger 
than Aqcd- Hyer ||145| | suggests, for example, to take for the factorization scale fip 



ma.x{l/bi}. In this case, 7max appears in Eq. ( | 176| ) only if all three bi tend to 1/Aqcd at 
once. But then at least one of the exp[—s{C,i,bi,Q)] drops to faster than any power of 
In (^1/6;Aqcd)- Apparently, Hyer's choice of the factorization scale avoids singularities of 
the form of Eq. (|175|) , but seems physically implausible. Since he gives numerical results for 
the proton form factor in the time-like region only, his results cannot be compared directly 
with those presented here. 

Another option, and actually the one favored here, is to use a common IR cutoff not 



only for the evolution of the wave function but also in the Sudakov exponent ||78| , |80| , |25 



Indeed, for a common cutoff b, the Sudakov factors always compensate the ag-singularities, 
i.e., if, for a given /, it happens that one Sudakov function is in the dangerous region of 
phase space, ^/ < -\/2Aqcd/Q, ^Aqcd 1, at least one of the other two Sudakov functions 
lies in the region > \/2Aqcd/<5, ^Aqcd ^ 1 7^ and therefore provides sufficient 
suppression, as outlined above. In particular, we favor b = max{6;} as the optimum choice 
("MAX" prescription), since it does not only lead to a regular integral but also to a non- 
singular integrand of the form factor as well. The Sudakov factor e~'^^ subject to the "L" 
and "MAX" prescriptions is plotted for a specific quark configuration in Fig. This figure 
makes it apparent that the Sudakov factor in connection with the "MAX" prescription is 
unencumbered by singularities in the soft regions. As a consequence of the regularizing 
power of the "MAX" prescription, the perturbative contribution to the proton form factor 
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^2 • Aqcd 

Figure 25. The Sudakov factor 6"*^^ vs 6iAqcd and 62AQCD evaluated at Q = 30AqcD: and 
Xi = x\ = 0.9, X2 = X3 = X2 = x'^ = 0.05 assuming a linear quark configuration (i.e., h±i and h±2 
are parallel to each other). The upper and lower figures correspond to the "L" and "MAX" prescriptions, 
respectively. 
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Figure 26. The proton magnetic form factor as a function of 6cAqcd- The curves shown are obtained 
at Q = SOAqcd for the COZ distribution amplitude. The solid line corresponds to the "MAX" prescription 
including evolution. The dotted (dashed, dashed-dotted) line represents results using the "H-SS" ("MAX", 
"L") prescription, by ignoring the intrinsic dependence, and also evolution. 



Eq. (167) saturates, i.e., the results become insensitive to the inclusion of the contributions 
of the soft regions. A saturation as strong as possible should be regarded as a prerequisite 
for the self-consistency of the perturbative contribution. 

To demonstrate the amount of saturation, we calculate the proton form factor through 
Eq. ( [167| ), employing a cutoff procedure to the fe^-integrations at a maximum value be- In 
Fig. ^the dependence of on be for the three choices, labeled "L", "H-SS", and "MAX" is 
shown, using, for reasons of an easier comparison with previous works, the COZ distribution 
amplitude and ignoring evolution. El 

As one sees from this figure, the "MAX" prescription empowers saturation, i.e., the soft 
region 6cAqcd > 0.7 does not contribute to the form factor substantially. In fact, already 
50 % of the result are obtained from the regions with 6cAqcd < 0.48, while as increases to a 
value of 0.95 at &cAqcd ~ 0.48. This indicates that a sizeable fraction of the contributions 
to the form factor is indeed accumulated in the perturbative region. 



^^Evolution has been dispensed with to avoid a concomitant singularity in Q^G^ as 6cAqcd 1 
when imposing the "L" and "H-SS" prescriptions. 
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Unfortunately, this saturation is achieved at the expense of a rather strong damping of 
the perturbative contribution to the proton (and neutron) form factors. Using the other two 
prescriptions ("L" and "H-SS") and ignoring evolution, larger results for G^j, and can 
be obtained, but there is no indication for saturation: The additional contributions to the 
form factor gained this way are accumulated exclusively in the soft regions, i.e., for values 
of &cAqcd near 1. Hence, this would be a rather deceptive improvement since the criticism 
of Isgur and Llewellyn-Smith and of Radyushkin [T^JT^ would again apply. 

These findings are in evident contradiction to Li's results (figure 5 in |^9[) for which 
an acceptable saturation was claimed. On the other hand, the saturation behavior of the 
proton form factor calculated by Hyer ||145|| in the time-like region is confirmed. A saturating 
behavior of the form factor should be regarded as a stringent test for the self-consistent 
applicability of pQCD. Therefore, calculations, which accumulate large contributions from 
soft regions (large be), can hardly be considered as theoretically legitimate, even if they seem 
to be phenomenologically successful. 

The role of the evolution effect subject to the "MAX" prescription is also effected in 



Fig. It is evident from this figure that the effect of evolution is large, though finite, owing 
to the strong suppression provided by the Sudakov factor (recall that the factorization scale 
is hf = 1/b). The significant feature of the evolution effect is that it tends to neutralize the 
influence of the IR cutoff. Thus one may trade larger values of the proton (neutron) form 
factor for a slightly worse saturation. 

One may criticize the use of the one-loop formula for as in the soft region. The "true" as 
coupling may likely differ from that in the soft region. However, the saturation behavior due 
to the "MAX" prescription (see Fig. ^6|) reveals that there is practically no contribution to 
the form factor from soft regions (say for bAqcn > 0.8). The Sudakov factor together with 
the intrinsic transverse-momentum dependence of the wave function suppresses this region 
completely. Therefore, it is irrelevant what description of as one is using in that region, and 
hence more involved as parameterizations can be avoided. 

We now comment on the connection between the modified convolution scheme, employed 
here, and approaches (e.g., ||120|| ) which make use of an effective gluon mass [121] to regular- 
ize the running coupling constant and gluon propagators at small values of Q^. The purpose 
of the following discussion is to show the utility of the k_i_ (or transverse configuration space) 
representation in compensating the divergent contributions of the soft region of phase space. 
According to Eq. (|152|) or equivalently Eq. (|167|) , all values of have been integrated out, 
in particular those corresponding to large distances, i.e., such of order I/Aqcd, which are 
not governed by perturbative QCD. As discussed above, once ~ Aq^j-,, the external 
probe (the virtual photon) no longer "sees" the substructure of the nucleon, it rather sees 
the "diffuse" nucleon as a whole. This means that the strength of the dynamical coupling 
constant ceases to increase and "freezes" at some scale < = 4mg (recall Eq. (|117|) ), 
below which the long-range forces become saturated [^. The scale can be fixed by 
the requirement that it should lead to the experimentally observed form factors. Then the 
leading-order coupling constant is replaced by the modified expression | p,21| , p,22|| given by 
Eq. (|117|) , where the arguments of as are = Q^^i^l or Q^^i^l with ^/ = (1 — ^/). In this 
case, the IR-protected as remains small enough that perturbation theory still applies. Fur- 
thermore, when the longitudinal momentum fractions ^/ tend to zero, the gluon virtuality 
remains finite and equal to M^, and well above the Landau ghost pole. On the other hand, 
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Table 12. Simulation of Sudakov suppression and intrinsic transverse momentum by an 
IR-cutoff parameter M in the running coupling constant at the one-loop order. The values 
of Q^G^ are calculated with the COZ distribution amplitude and the r.m.s. transverse momen- 
tum (/ci)^/2 = 271 MeV. 



g2 [GeV2] 


Q^Glj [GeV^] 


M [MeV] 


as(9^) 


4.0 


0.3593 


330.9 


1.147 


6.0 


0.3823 


343.8 


1.079 


8.0 


0.3962 


353.4 


1.035 


10.0 


0.4042 


361.6 


1.001 


12.0 


0.4094 


368.5 


0.974 


14.0 


0.4134 


374.6 


0.953 


16.0 


0.4164 


380.0 


0.934 


18.0 


0.4186 


384.9 


0.919 


20.0 


0.4195 


389.6 


0.904 


22.0 


0.4197 


394.1 


0.891 


24.0 


0.4193 


398.4 


0.879 


26.0 


0.4185 


402.5 


0.868 


28.0 


0.4176 


406.4 


0.857 


30.0 


0.4165 


410.2 


0.848 



when becomes large enough, the form factor becomes insensitive to the large distances 
or small momenta, characterized by M. Thus, Eq. ( |117|) simulates the effect of the Sudakov 
factor in suppressing contributions from large values of |b_L|. However, in contrast to the 
Sudakov damping factor, which has universal character, the IR-cutoff parameter M has to 
be adjusted to the external momentum. This becomes evident from Table |T2|, where the 
results for the Sudakov-suppressed magnetic proton form factor are simulated in the range 
of momentum transfers 4 to 30 GeV^, using again as reference model the COZ distribution 

amplitude. Evolution is taken into account via the factor —r^ r , where /ig = 1 

GeV. Indeed, while the form factor begins to scale at 14 GeV^, the values of the ad hoc 
parameter M still increase with Q^, in order to provide sufficient IR-protection. Such a 
behavior of M can hardly be associated with the renormalization group controlled evolution 
of a dynamical gluon mass with a definite anomalous dimension [[121| . 

This exercise demonstrates that IR regularization by such techniques is conceptually 
inferior to the inclusion of the explicit dependence of within the modified convolution 
scheme. Perhaps the most significant feature of the latter approach is that it remains valid 
in any regime of momentum transfer without recourse to ad hoc assumptions. 



2. Numerical form-factor analysis 
In this section we present numerical results for the proton and the neutron form fac- 



tors obtained in [ 78, 80, 25 1. In these calculations the "MAX" prescription is employed, with 



85 



> 

ID 

o 



a 



0.6 - 



0.4 - 



0.2 - 












20 



40 



60 



80 



100 



[GeV^] 



Figure 27. The influence of the intrinsic transverse momentum on the proton magnetic form factor 
within the modified convolution scheme. The curves shown are obtained for the COZ distribution amplitude 
by imposing the "MAX" prescription, and including evolution. The solid line represents the result without 
A;_L-dependence, whereas the dashed and dotted hues are obtained with (/c^)^/^ = 271 MeV and 600 MeV, 
respectively. 



evolution included, using Aqcd =180 MeV and /^o = 1 GeV. Before proceeding with the pre- 
sentation of the results, let us first investigate the effect of including the intrinsic transverse 
momentum in form-factor calculations. The fc_L-dependence of the nucleon wave function 
effectively introduces a confinement scale in the formalism, the importance of which may be 
appreciated by looking at Fig. 

This figure shows results, obtained for the COZ distribution amplitude without and 
with A;_|_-dependence, using two different values of {k'j_)^^'^. To describe the intrinsic k±- 
dependence, one can use Eq. ( |162| ) or, after Fourier-transforming to the transverse config- 
uration space, Eq. (|174|) . Notice that in Li's approach the Gaussian in Eq. (|174|) has been 



replaced by unity. The oscillator parameter a is determined in such a way that either the 
normalization of the wave function is unity (resulting into (^i)^^^ = 271 MeV for the 
COZ distribution amplitude), or by inputting the value of the r.m.s. transverse momentum. 
In the second case, a value of 600 MeV is used, which implies Psq = 0.042. As can be seen 
from this and the analogous figure for the neutron (Fig. the predictions for the form 
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Figure 28. The influence of the intrinsic transverse momentum on the neutron magnetic form factor 
within the modified convolution scheme. The curves shown are obtained for the COZ distribution amplitude 
by imposing the "MAX" prescription, and including evolution. The solid line represents the result without 
A;_L-dependence, whereas the dashed and dotted fines are obtained with (/c^)^/^ = 271 MeV and 600 MeV, 
respectively. 



factor are quite different for the three considered cases. 

In fact, the intrinsic fc^-dependence of the nucleon wave function leads to further sup- 
pression of the perturbative contribution, which becomes substantial if the r.m.s. transverse 
momentum is large. On the other hand, this suppression is accompanied by an increasing 
amount of saturation because the Gaussian distribution, Eq. ( |174|) , also provides suppres- 
sion; predominantly of contributions from the soft regions, viz., the large 6-regions. In 
contrast to the Sudakov factor, however, this suppression is Q-independent (there is no 
evolution). The interplay of the two effects, Sudakov suppression and intrinsic transverse 
momentum, leads to a different Q-behavior of the form factor, depending on the value of the 
r.m.s transverse momentum, as can be seen from Figs. ^ The Q-dependence beyond 
10 GeV^ is rather weak, being approximately compatible with dimensional counting (mod- 
ulo renormalization group generated logarithmic corrections). For very large values of Q^, 
say, beyond 1000 GeV^, the three curves have already approached each other within 10% 
accuracy. This happens when the Sudakov factor dominates the Gaussian /c^-distribution 
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Figure 29. The proton magnetic form factor vs Q^. Data are taken from |154,3j]. The data are 
represented by black dots, whereas those for Ff are indicated by open circles. The theoretical results are 
obtained using the "MAX" prescription including evolution and normalizing the wave functions to unity. The 
shadowed strip indicates the range of predictions derived from the set of distribution amplitudes determined 
in 1^,^ in the context of QCD sum rules (see text). The solid (dashed, dotted) line corresponds to the 
COZ (heterotic, optimized COZ) distribution amplitude. 



and selects those configurations with small inter-quark separations. In this region, which one 
may consider as the pure perturbative region, the results for the form factor are independent 
of the confinement scale introduced by the r.m.s. transverse momentum. 

The penalty of the additional suppression of the perturbative contribution caused by 
Eq. (|174|) is mitigated by the advantage of promoting the perturbative contribution to be- 
come self-consistent by the in situ IR protection due to the Sudakov form factor. This is 
clearly indicated in the enhanced amount of saturation with increasing r.m.s. transverse mo- 
mentum. Adapting the criterion of self-consistency, originally suggested by Li and Sterman 



77| for the pion case, namely that 50% of the results should accumulate at moderate values 
of the coupling constant, say, < 0.5, we find for the nucleon form factors self-consistency 
for ~ 7 GeV^ (still referring to the COZ distribution amplitude). 

The final results are shown in Fig. ^ (proton) and in Fig. |30| (neutron) for different 
nucleon distribution amplitudes. 

Though we have not considered the pion form factor in the modified convolution scheme. 



we find it instructive to include it in Fig. RT] which shows the (corrected) calculation by 
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Figure 30. The neutron magnetic form factor vs Q^. The theoretical resuhs are obtained using the 
"MAX" prescription including evolution and normalizing the wave functions to unity. The shadowed strip 
indicates the range of predictions derived from the set of distribution amplitudes determined in [5l| in the 
context of QCD sum rules (see text). The solid (dashed, dotted) line corresponds to the COZ (heterotic, 



optimized COZ) distribution amplitude. The data are taken from [154| 
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Figure 31. The pion form factor in the space-hke region including Sudakov effects and those due to 
the intrinsic A;_L-dependence of the pion wave function evaluated for Aqcd = 200 MeV. The solid line 
shows the result for the Chernyak-Zhitnitsky amplitude ^cz = 30xiX2{xi — 2:2)^, and the dotted line that 
for the asymptotic solution to the evolution equation, <I>as = 6x1X2. The data are taken from [p2|. 



Jakob and Kroll [81 



A set of 45 nucleoli distribution amplitudes [§T],^ have been investigated, distribution 
amplitudes which respect the QCD sum-rule constraints p8| , p7| with comparable values. 
The results for the various distribution amplitudes - or more precisely wave functions, since 
their intrinsic transverse-momentum dependence is taken into account - obtained under 
the "MAX" prescription with evolution included, form the shaded area shown in Figs. |2^, 
PU| . Note that all wave functions are normalized to unity and that the corresponding r.m.s. 
transverse momenta vary between 267 MeV and 317 MeV (see Table [Tl|). 

The theoretical form-factor predictions span a "band" congruent to the "orbit" of solu- 
tions found previously in [|5T| , |52| and discussed already in detail in this report. The upper 
bound of the "band" corresponds to the distribution amplitude COZ"p, which yields the 
maximum value of the form-factor ratio jCJ^I/G^ = 0.4881 in the standard convolution 
scheme. The lower limit of the "band" is obtained using the distribution amplitude "low" 
(sample 8 in ^2[) (cf. Table ^) with |G^|/G^i = 0.175. Explicitly shown are the results for 
the COZ distribution amplitude, its optimized version COZ°p', and the heterotic distribu- 
tion amplitude We note that the differences among these curves practically disappear 
already at about = 80 GeV^, despite the fact that these amplitudes have distinct geomet- 
rical characteristics. This behavior effects once more that momentum evolution at larger 
values is enough to wash out shape peculiarities due to the truncation of the eigenfunctions 
series. 

Because the true valence Fock state probability is likely much smaller, or invariably the 
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r.m.s. transverse momentum larger than a value of order 300 MeV, the "band" describes 
rather maximal expectations for the (leading-order) perturbative contributions to the form 
factor; at least for proton wave functions of the type we considered. Comparison with the 
experimental data reveals that the theoretical predictions amount, at best, to approximately 
50% of the measured values. This is the benchmark against which we have to discern 
novelties and aberrations. 

Concerning the proton form factor, we note that, since we are calculating only the 
helicity-conserving part of the current matrix element, it is not obvious whether we should 
compare the theoretical predictions with the data for the Sachs form factor or with 
those for the Dirac form factor Ff. Therefore, we exhibit in Fig. EOlboth sets of data ||154 



for comparison. However, since the two sets of data differ by only 10%, our conclusions 
concerning the smallness of the theoretical results remain unaffected. 

The results for the neutron magnetic form factor are shown in Fig. |3^. Unfortunately, 
there is no form-factor data available in the region where the perturbative contribution 
is self-consistent. Yet the trend of the low data ||154]| seems to indicate that the size 
of the perturbative contribution is rather small. Measurements of the neutron magnetic 
form factor beyond 5 GeV^ are extremely important in order to check the validity of the 
theoretical predictions. 



One place to test these results beyond 5 GeV^ is in the data ||123|| for the differential cross 
sections for elastic electron-proton and electron- neutron scattering, ap and respectively. 
For small scattering angles, where the terms oc tan^(^/2) can be neglected, and for large Q^, 
the ratio (7„/ Cp becomes approximately proportional to the square of the ratio of the neutron 
to the proton magnetic form factor, provided the electric form factors are negligible. For the 
electric form factor of the proton this assumption has recently been verified experimentally 
| |123| |. Its neutron counterpart is experimentally still unknown beyond 5 GeV^, but many 
phenomenological models currently in use predict a small G^. We note that the low data 
on | ]154[ are compatible with zero. 

Combining the calculations for the proton with those for the neutron, we can extract 
theoretical predictions for (T„/ ap, using the same set of nucleon model distribution amplitudes 
as before. The results are shown in Fig. ^ (shaded area) in comparison with available data 
|123|| . From this figure we see that the measured values of Onjop enter the estimated range 
already at ~ 8 GeV^. [The corresponding values of the ratio —G\IG\j^, allowed by 
this analysis, range between -0.2 and O.5.] It is important to emphasize that the ratio of 
the magnetic form factors is the first observable for which the modified convolution scheme 
yields predictions which, albeit in a slightly mo del- dependent way, indicate overlap with 
the existing data ||123|| . This tentative agreement occurs at data points corresponding to 



the largest momentum transfers measured, where, incidentally, the presented theoretical 
calculations become self-consistent. 

The various model (nucleon) wave functions considered lead to self-consistency of the 
perturbative contribution, meaning that 50% of the results are accumulated in regions where 
ctg < 0.5, i.e., in the range of between 6 and 10 GeV^. 

Closing this chapter we note that using other values of the r.m.s. transverse momentum 
but the same set of distribution amplitudes, similar "bands" are obtained with about the 
same relative widths. Even for zero r.m.s. transverse momentum (see Fig. |2^), the "band" 
does not overlap with the data. 
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Figure 32. The ratio (Jn/o'p of the differential elastic electron- neutron to electron- proton cross section 
vs at scattering angles of 10°. The shaded area and the model distribution amplitudes for the nucleon 
correspond to those shown in Figs. E9 and pO. The data are taken from ll22 
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Our present knowledge of the proton wave function resides only on QCD sum-rule calcu- 
lations and is thus rather limited - contrary to the pion case. Therefore, we present two sets 
of results, such with a minimal r.m.s. transverse momentum (corresponding to P^q = 1), 
and such with 600 MeV (corresponding to P^^ ^ 0.04) in order to demonstrate the role 
of this parameter. The probability of the "true" valence Fock state wave function, say, at 

= 10 GeV^, is certainly not 1 but smaller, probably close to the value of the other set. 
However, according to Botts and Sterman []7B[ (see also ||2^), one may view the Sudakov 



form factor as part of the wave function, describing the perturbative part of the transverse 
momentum distribution. Accepting this idea and interpreting the case P^q = 1 as the soft 
wave function in the region — 1 GeV^, the product of wave function and Sudakov form 
factor at about = 10 GeV^ leads to Psq = 0.21 and {k±^y/^ ~ 500 MeV. Thus, in order 
to avoid double counting, it is perhaps reasonable to start with the intrinsic /c_L-dependence 
of about {k±Y/^ ~ 300 MeV and Pgq = 1. 

We finally note that the value of the proton wave function at the origin, /n, is also 
burdened by large errors, which in turn induce further uncertainties in the form factor 
predictions. In present lattice theory calculations the value Z^* = (0.29 ± 0.6) x 
10~^ GeV^ was found. Literal use of such a value would lead to form factor results even 
smaller than the predictions shown in Figs. 12^, BD. 



3. Higher-order nucleon distribution amplitudes 

To effect that the depletion of the perturbative contribution is not the consequence of 
truncating the nucleon distribution amplitude at the level of second-order eigenfunctions, we 
show in Fig. |3^ predictions for the proton and the neutron magnetic form factors calculated 
with nucleon distribution amplitudes, determined by Schafer [Q, which incorporate Appell 



polynomials of order three.ll3 The solid lines show the results for the amplitude Sch II (in 
Schafer's notation) which deliberately incorporates third-order terms. It is clearly obvious 
that both form factors G^j and GJJj overshoot the data and have, in addition, a wrong Q"^ 
evolution. Besides, their saturation behavior deteriorates, for they still increase at as 
large as 50 GeV^ and insensitivity to the soft region biAgco 1 sets in at ~ 17 GeV^, 
which is a much larger scale compared to that found for the second-order amplitudes: < 
10 GeV^. As we have discussed in detail in previous publications [p5| , |103| , ^ , p] ,pl|1 , and 
also in previous sections of this report, the moment sum rules are not stringent enough to 
exclude this sort of amplitudes. Additional criteria have to be imposed in order to filter out 
physically meaningful solutions, i.e., those which ensure dominance of the lowest leading- 
order contributions. 

Along similar lines of thought, Schafer |^0[ has determined a collection of amplitudes 



which, although including third-order Appell polynomials, effectively resemble those of poly- 
nomial order two. The dotted lines in Fig. p3| exemplify the form-factor predictions obtained 
with such amplitudes. To be specific, the underlying amplitude (denoted Sch IVb) was deter- 



^^The expansion coefficients i?„ for each of these amplitudes can be readily obtained from the 
corresponding moments via Eqs. (p2t), (p3[). 
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Figure 33. Plots for the proton and neutron magnetic form factor within the modified convolu- 
tion scheme calculated with nucleon distribution amplitudes, determined by Schafer pO| , which include 
third-order eigenfunctions (Appell polynomials). As in the previous figures, the "MAX" prescription and 
the value Aqcd = 180 MeV were used. 
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Figure 34. Ratio of tlic nuclcon magnetic form factors computed with the Schafer amplitudes vs the 
momentum transfer Q^. 



mined under the proviso of a "smoothness" criterion with the purpose of excluding spurious 
oscillations ("wiggles") caused by the third-order (Appell) polynomials. This criterion is 
quite restrictive and, as one sees, the obtained results have the same overall quality as those 
computed with amplitudes truncated after taking into account bilinear combinations of lon- 
gitudinal momentum fractions. This indicates that the truncation at the second polynomial 
order is justifiable, since the shape of the corresponding amplitudes is a characteristic prop- 
erty of the entire series and that the errors (and amount of cutoff dependence involved) are 
of sub-leading importance if properly incorporated. For completeness. Fig. shows also 
form factors, computed with an amplitude (termed Sch III and represented by dashed lines) 
which exhibits unphysical oscillations. The form-factor ratio — GJJj/G^, as a function of the 
momentum transfer Q^, for the three amplitudes Sch II, Sch III, and Sch IVb is shown in 

Fig. pa 



Global pattern of nucleon distribution amplitudes 



It is remarkable that the collective pattern of solutions to the QCD sum rules ||38| , |37 



found within the standard convolution scheme |51,52|, pertains to the inclusion of transverse- 



momentum contributions comprising the Sudakov factor and those due to the intrinsic trans- 
verse momentum (see Fig. |35| ). 

Indeed, the solutions (nucleon distribution amplitudes) arrange themselves again across 
an "orbit" in the (i?4, — G^/G^m) plane which is, however, somewhat shifted as com- 
pared to the original one (Fig. 0). Now the "orbit" is slightly Q^-dependent, though 
its coherent structure is not destroyed but changed in a global way. The solid and 
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Figure 35. The relation between tlie ratio R = of the magnetic nucleon form factors and 

the expansion coefficient of the AppeU polynomial decomposition of the nucleon distribution amplitude, 
including the effect of the Sudakov form factor and intrinsic transverse-momentum dependence. The results 
are obtained at = 30 GeV^, employing the "MAX" prescription. The superimposed solid line is an 
empirical polynomial fit similar to the original one given in |5l|| within the standard convolution scheme. 
The dashed line serves to illustrate the dependence on the momentum scale (Q^ = lO"' GeV^). 
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dashed lines correlate to the calculations within the modified convolution scheme presented 
above at two different scales, viz., at = 30 GeV^ and = 10^ GeV^, respectively. 
The orbit at = 30 GeV^ (solid line) can be characterized by the empirical relation 
-Gm/GIi = 0.426 - 9.91 x 10-^54 - 4.27 x IQ-^Bj + 4.59 x IQ-^Bl which complies with 
that found within the standard convolution scheme. The dashed line in Fig. |3^ represents a 
similar fit for = 10^ GeV^. We observe that as increases the modified orbit approaches 
again the original one. This observation is double-edged: it means not only that the orbit 
structure is stable and can be used in both convolution schemes in order to exclude unphys- 
ical solutions, for example, Sch III, it also supports the view that in the axial gauge, 
operates like a finite wave-function renormalization factor |^5[, as outlined above. This con- 
ception may be used reversely to model hadron distribution amplitudes with an improved, 
i.e., "softened" endpoint behavior relative to those determined by QCD sum rules alone. 



VII. CONCLUSIONS AND OUTLOOK 

We have reviewed a large body of exclusive processes which relates to testing the validity 
of perturbative QCD, in particular factorization theorems, and the universality of non- 
perturbative hadron distribution amplitudes modeled on the basis of QCD sum rules. 

Within the standard convolution scheme, our "heterotic" distribution amplitude for 
the nucleon - which amalgamates the pivotal features of COZ-type with GS-type 



33| , |53|J35| distribution amplitudes - successfully correlates a wide variety of unrelated ob- 
servables to the data without tuning free parameters from case to case and without using 
additional phenomenological constraints ||103|| . This agreement is non-trivial; no other the- 
oretical model matches the available experimental data nearly as well. 

By analyzing the QCD sum-rule constraints on the moments of nucleon distribution am- 



plitudes [p8| , p7| with the aid of a "hierarchical" criterion - that takes into account the 



higher accuracy |31| , |35|1 of lower-order sum rules - we were able to systematically explore the 



characteristic features of distribution amplitudes complying with the sum-rule constraints 



5T| . Our extensive analysis does reveal definite evidence for a global pattern of such distri- 
bution amplitudes across an orbit in the plane spanned by the ratio R = \G^\/Gy^ of the 
magnetic form factors of the nucleon and the expansion coefficient B4, the latter projecting 
on to an antisymmetric eigenfunction of order two of the nucleon evolution equation. This 
orbit is finite, ranging from distribution amplitudes with a ratio value R around 0.1 (het- 
erotic region) up to COZ-type solutions with R < 0.49, where the orbit terminates (fixed 
point). Proposed nucleon distribution amplitudes which fall outside the orbit turn out to 
be unphysical, irrespectively of their functional representation in terms of eigenfunctions. 
Hence, we are inclined to believe that the orbit captures a genuine feature of the true nu- 
cleon distribution amplitude, i.e., of the whole eigenfunctions series, resulting in turn into 
an insensitivity to particular schemes of truncation. We were able to formulate the concept 
of "heteroticity" in mathematical terms by introducing a "hybridity" angle to quantify the 
mingling of COZ- and GS-type distribution amplitudes, thus completing the classification 
scheme of nucleon distribution amplitudes. 

Factorizing the exclusive-scattering amplitude for the nucleon form factors according to 
the modified convolution scheme - which includes gluon radiative corrections in the form of 
Sudakov-like damping exponentials - the orbit retains its structure and changes only slightly 
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as a whole. For momentum-transfer values on the order of 10^ GeV^ it rotates back to the 



one obtained within the standard convolution scheme pSipOf . On the phenomenological 
side, limiting the value of the form-factor ratio by experiment, one could use the orbit to 
extract the theoretically appropriate nucleon distribution amplitude. 

Similar ideas of "heteroticity" could be applied to QCD sum-rule analyses pi| , ^ which 
constrain the moments of the A"*" distribution amplitude. A theoretical model was derived 
P3| , termed again "heterotic", that provides the best agreement with existing data PB|J^ 
on the (electromagnetic) transition form factor G^. 

Overall, the current status of the QCD sum rules program to exclusive high-momentum 
processes produced a rich body of theoretical ideas and experimental tests that could be 
utilized to scrutinize new developments, like the use of nonlocal condensates. 

But there are problems, both technical and conceptual. The inversion of few moments 
cannot fix hadron distribution amplitudes uniquely |]53|j35| , p!7| . On the other hand, the calcu- 
lation of higher-order moment constraints is burdened by intrinsic theoretical uncertainties 



3T| , p!7[| ; e.g., power corrections grow with the moment order thus overwhelming the pertur- 
bative contribution. Clearly, from the theoretical point of view, a QCD sum-rule method 
to calculate hadron distribution amplitudes as a whole would be desirable, but is still to be 
developed. First steps in this direction ]T7|PB|] are promising but have to be further pursued 
to include baryons. 

On the perturbative side, we were able to calculate a total of 55 next-to-leading eigenfunc- 
tions of the nucleon evolution equation, tabulating results up to polynomial order M = 4 
(Tables ^ |^). We investigated the corresponding anomalous dimensions of leading- twist 
three and estimated their large-order behavior in successive steps from M = 9, M = 20, 
M = 38, up to M = 150. This behavior can be best reproduced by a logarithmic fit, thus 
confirming previous expectations based on exponentiation of one-loop results. 

Repeating the form-factor calculations for the nucleon within the modified convolution 
scheme |7^J7^, which incorporates the effect of (soft) gluon radiative corrections by means 
of Sudakov-type form factors |^6[, we showed how the screening of as singularities in the 
endpoint region can be enforced by employing an infrared regularization prescription, based 
on the assumption that quarks at large transverse separations act coherently and thus can- 
not radiate soft gluons ("MAX" prescription). This prescription renders the calculations 
finite and reinstates the validity of the perturbative treatment providing sufficient infrared 
protection, since, even with evolution, the integrand in Eq. (|167|) remains finite []78| , |80| , p5 



A significant feature of this treatment is that the nucleon form factors saturate, i.e., be- 
come insensitive to the contributions from large transverse separations. The other choices 
of infrared-cutoff prescriptions ("L", "H-SS") p^|145| , p!15[] do not lead to saturation. 

However, a heavy price is paid. A reliable saturation and infrared protection of the form 
factors are achieved at the expense of a strong reduction of the perturbative contribution. 
The damping of the nucleon form factors becomes enhanced if the intrinsic transverse mo- 
mentum dependence of the nucleon wave function (see Fig. ^ and Fig. is taken into 
account. This has been done by assuming a non-factorizing x and /c^-dependence of the wave 
function of the Brodsky-Lepage-Huang type ||147|| , and fixing the value of (A;^.)^''^ either via 
the valence quark probability Psq or by inputting the value 600 MeV by hand ||146|| . A 
remarkable finding is that the form factors calculated this way show only a mild dependence 
on the particular model distribution amplitude used. 
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The perturbative contribution to the form factor becomes self-consistent in all cases for 
momentum transfers larger than 6 to 10 GeV^, albeit the actual value of the onset of self- 
consistency depends on the particular wave function and the r.m.s. transverse momentum 
chosen. Self-consistency is defined such that 50% of the result are accumulated in regions 
where is smaller than 0.5. 

Comparing these theoretical results with the data, it turns out that they fall short by 
at least 50%. This is true not only for the COZ distribution amplitude (which we have 
exemplarily used as reference model) but actually for the whole spectrum of amplitudes 
determined in |pT| , |5^ that satisfy existing sum-rule requirements. Depending on the actual 
value of the r.m.s. transverse momentum, the reduction of the perturbative contribution 
may be even stronger than 50%. 

This damping may perhaps be counteracted by taking into account pairing effects in the 
hadron, effects which go beyond the simple impulse approximation. The latter assumes that 
the interaction of the external photon takes place with just one parton at a time. One may 
argue that non-perturbative vacuum fluctuations, e.g., due to instantons [|100|| , give way to 



parton correlations and provide additional finite renormalizing factors I^S] greater than 1 
to bridge the gap to the experimental data at laboratory momentum transfer. Or perhaps 
higher-order corrections of perturbation theory giving rise to a rather large K-factor of the 
order 2, as found for other large-momentum transfer (inclusive) processes [ 155|| , may im- 
prove the theoretical predictions. Furthermore, the failing of the leading-order form-factor 
contribution to reproduce the existing data, may be a signal that soft contributions (non- 
logarithmic contributions) not accounted for so far in the formalism should be included. 
Such contributions include, e.g., improved and/or more complicated wave functions com- 
prising higher-twists, orbital angular momentum, higher Fock components, quark-quark cor- 
relations (diquarks), higher- dimensional quark/gluon condensates, quark masses, etc. Also 
remainders of genuine soft contributions, like vector-meson-dominance terms or the non- 
factorizing overlap of the soft parts of the wave functions (Feynman contributions), may 
still be large at accessible momentum transfers. The rather large value of the Pauli form 
factor F2 around 10 GeV^, as found experimentally (second reference in ||123|| ), indicates 
that sizeable higher-twist contributions still exist in that region of momentum transfer. One 
may suspect similar or even larger higher-twist contributions to the helicity non-flip current 
matrix element controlling Fi and Gm- 

Once all this is accomplished, we shall probably be able to make contact with the data 
in a more meaningful way. 

A real theoretical breakthrough is, however, not a project that simply reconciles form- 
factor calculations with experimental data. The basic challenge is finding a way to describe 
an exclusive process in regimes where the impulse approximation becomes invalid. This hap- 
pens when for some reason the interaction time becomes comparable with that of hadron 
formation. Then, speaking in terms of pictures, the dynamics is not a "snapshot" (the 
essence of the impulse approximation) but rather a "blurred" (i.e., "smeared") field configu- 
ration which extends to a large range of scales linking the initial with the final hadron state 
(scale intrusions). The consequence is that factorization of the amplitude breaks down and 
a short-distance part amenable to perturbation theory can no more be isolated. This is the 
case in the kinematic endpoint region when partons can become very "wee" . The appearance 
of soft propagators in the hard part of the amplitude signals infrared sensitivity and means 
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that scale locality is lost. These problems cannot be resolved insisting to use valence quarks. 
These are not asymptotic fields and hence cannot be part of the true asymptotic dynamics 
(Hamiltonian of the system). One should use instead quark lines which are quark-gluon 
composites and can inherently accommodate soft modes within an effective approach. The 
goal is to transcend perturbation theory by providing all-order expressions in the coupling 
constant for dynamical quantities, like Green and universal vertex functions. Work in this 
direction is in progress P^ , |156|| . 
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APPENDIX A: EIGENFUNCTIONS OF THE NUCLEON EVOLUTION 

EQUATION 



The explicit form of the matrix Uij^ki which relates the nucleon evolution kernel to the 
polynomial basis \kl) is 



V\kl > _ ] 
2w{xi) ~ 2 Jo 

h 1 

rf"" rv*'' 



w{Xi) 
1 



k+l 1 

-3E7 



(/c + l)(/c + 2) ^2 J (/ + l)(/ + 2) 



i+i 1 



2^-7^ r^x' 



i=l j=0 ^ 

1=1 j=0 

i+j<M 



El 

j=0 

:;) 

' k 
J 



-lyxt 



I- J 

I - 1 

(k~3 



rn 



m=2 



k — i 



m=2 



E 



(Al) 



with u'(xj) = X1X2X3 — Xi{l — Xi —X3)x3. The operation of y on a polynomial |A; /) of degree 
M = k + l projects on another polynomial of degree i+j< M. This projection is expressed 
through the matrix U . In the present work U is represented in terms of AppeU polynomials 
which are defined on the triangle T = T{xi, X3) with xi > 0, 0:3 > 0, xi + < 1 as follows 



-r / \ -U r{xi + m)T{x3 + n) 

^mniXO, Xl, Xs; Xi, X3) = W {Xo, Xl, Xs] Xl, X3) — — r 

r(xi) r(x3) 

am+n 



dx^dx'^ 

Here F denotes the F-function with 

r(x.) 



Xi ixi- i)---(xi + i -1) 



(A2) 



(A3) 



(A4) 



and the weight function is defined by 

w(xo,Xi,X3;a:^i,a;3) = x^'~^ x^''~'^ {1 - xi - xs)^"'^'"^^ . 
The polynomials Trnn are of order M — m + n and analytical, provided the Xi a-^e real, i.e., 
3^(Xi) > , 3?(X3) > , 3fJ(xo) > 3?(Xi + Xs) - 1 • (A5) 
The scalar product of an arbitrary function f{xi,X3) on the triangle T with Trnn gives 
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{fl'F mn) 



X 



Xo 



1 - Xi - X3 



_ ^^^m+n+X()-Xl-X3 



(A6) 



where a "bracket" notation, resembling that in Chapter [ill C 3| , is used. If the function 
f{xi,x^) is another Appell polynomial, then the following two relations hold: 



(1) m + n ^ m' + n' 



^•F m'n'\-F mn) 



(A7) 



meaning that Appell polynomials of different degree are mutually orthogonal. 

(2) m + n = m' + n' 



^•F mn\F m'n') 

dxi dX3 



-1 



G 



^rn+Xl-l3,n+X3-l(i 



r(xi)r(x3) dxfdx^ 

m+n+xo-Xi-X3 



Xi - Xs) 



r(xi) r(x3) r(xo + m + n - xi - X3 + 1) 

r(xo + 2m + 2n + 1) 



-1) 



m+n ^ -I m'n' 



dx^dx^ 



(A8) 



meaning that Appell polynomials of the same degree are not a priori orthogonal to 
each other. 

In the nucleon case, M = m + n = 0, 1, 2, 3, . . . equals the total number of derivatives 
in the interpolating three-quark operators between the nucleon and the vacuum. In order 
to ensure that the asymptotic behavior of the nucleon distribution amplitude predicted by 
perturbative QCD is described by the amplitude ^tJyXi) = 120x1X2X3, the weight function 
of the Appell polynomials has to be w{xi) = X1X2X3 with X2 = 1 — Xi — X3 because of 
momentum conservation. Hence, Xi = Xs = 2 and Xo = 0, so that the relevant Appell 
polynomials are 



^i^)(5,2,2;xi,X3) 



WiXi 



k+l<m+n 

E (mn) 
k,l 



1 

m + 1)!(?7, + 1)! dx'idx'^ 
kl) . 



X 



m+l ™n+l 



X q 



'1 - Xi - X3 



im+n+1 



(A9) 



The symmetrized Appell polynomials used in the text are defined by the differential equation 

r(2 + m)r(2 + n) 



J'm„(Xi,X3) = XiX3(l -Xi -X3)- 

Qm+n 



,m 



X 



m+l n\ 



Xn 



'1 - Xi - X3 



,m+n+l 



dxi^dx^ ± dxidx^ 

where the plus sign refers to the case m > n and the minus sign to the case m < n. 



(AlO) 
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APPENDIX B: ELASTIC FORM FACTORS OF THE NUCLEON 

The functions P{Bn) and I"'{Bn), entering, respectively, the calculation of the proton 
and neutron Dirac form factors are: 

P = 1400 So + + 1800 Bo B, + '^^^^ + 1200 B^+ 

^ „ 22000 Si S3 „ „ 18800^3^ 1000 So S4 
6600 So S3 + ^-^^ + 4800 S2 S3 + ^ 3"^" 

2600 Si S4 2200 S2S4 4600 S3 S4 2600 S42 „ „ 

^ H - 1400 So S5- 

27 9 9 243 ° ^ 

3500 Si S5 1100 S2S5 5900 S3 S5 4100 S4S5 7700 Sg^ 

9 3 3 ^ 81 ^ 27 

I35S0S6 855^^1 So Sg 715SiS6 265^^1 Si Sg 

8 8 24 8 

1025 B-2 Bq 625 /TFT S. Sg 3775 S.-^ S(, 1175 S,, Sg 

8 8 24 24 ^ 

35S4Sg 65VwrS4Sg 1025 SgSg 325^/^755^6 

3 ^ 9 48 ^ 16 ^ 

94901 Sg^ 2629 ^/^^JlSg^ I35S0S7 855^/«JrSoS7 

768 ^ 768 8 ^ 8 

715S1S7 265v^SiS7 1025 S2S7 625y^S2S7 

24 8 8 8 

3775 S3 S7 4175^/SyrS3S7 35S4S7 65^/^«rS4S7 

24 ^ 24 ^ 3 9 

1025 S5S7 325VWIS5S7 „ „ ^ 94901 S7^ 

ZZ I jdr JD7 H 

48 16 ^ ^ 768 

2629 ^S72 8309 So Sg 109^SoS8 67525 Si Sg 

768 ^ 44 ^ 44 ^ 1188 ^ 

925v^SiS8 17197 S2S8 m ^B2Bs 120643 S3 Sg 

1188 132 132 396 

1643v^53S8 10939 S4S8 139v/^S4S8 94717 SgSg 

396 891 891 2376 

1717 ^B^Bs 331907 Be Bs 265651 ^SgSg 

2376 63 360 63360 

2b07^BeBs _ 3451 ^W^uh Be B^ _ 331907 S7S8 

63360 63360 63 360 ^ 

265651^^^7^8 2507 ^B^Bs 3A51 -./wtw^ B7 Bg 

63360 63360 63360 ~^ 

2390087 Sg^ 33287 yw^S8^ 8309 S0S9 109^SoS9 

348480 ^ 348480 ^ 44 44 ^ 

67525 Si Sg 925^SiS9 17197 S2S9 m ^B2B^ 

1188 1188 ^ 132 132 ^ 

120643 S3 Sg 1643v^S3Sg 10939 S4Sg 139y^S4S9 

396 396 891 ^ 891 

9 1717 S5 S9 1717 ^> B:, So 331907 S,; Sg 2507 S(, Sg 

2376 ^ 2376 63360 6 

265651 ^SgSg 3451^/wIIy^S6Sg 331907 S7S9 



39 

63360 
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265651 ywyiSySg 2507^ Br Bg _ y/UJrw^ Br Bg 

63360 ^ 63360 63360 ^ 

791 Eg Eg 2390087 Eg^ 33287 ^w^Eg^ 



1188 ^ 348480 348480 



, 2200 Bi^ 2800 Si Bo 
1800 - 1400 Bo Si + - 1800 Sq Sg 

„ 2 -. ^„ „ 22000 Si S3 _ „ „ 17000 Ss^ 
200 S2^ - 1000 So S3 ^ - 2000 S2 S3 

1000 So S4 2600 Si S4 2200 S2S4 4600 S3 S4 2600 S42 

3 27 ~'~ 9 9 ~'~ 243 ^ 

2000 So S5 3500 Si S5 500 S2 S5 6500 S3 S5 _ 4100 S4 S5 _ 

3 ' 9 ' 9 ' 9 81 
7250 Ss^ 725 Bo Be ?>2b^BoBQ 715 B^Bq 265 ^B^Bq 

81 ^ 8 ^ 8 ^ 24 ^ 8 

II5S2S6 195^S2S6 5815 S3 Se 4295vW53 56 

8 8 72 72 

35S4S6 65^B4B6 1705 B5 Be 935 ^B^Bq 

3 9 ^ 144 144 

91813 Se^ 2357 ^FTSe^ 725 Bo Br 325 y/w^Bo Br 

2304 2304 ^ 8 8 ^ 

715S1S7 265V«^SiS7 II5S2S7 195VS^S2S7 

24 8 ^ 8 8 ^ 

5815 S3 S7 4295^ S3 S7 35 B^ Br 65 y/wlB4 Br 

72 72 3 ^ 9 ^ 

1705S5S7^ 935V^S5S7 , ^ ^ 9I8I3S7' 

r I I Dp. Dr r 

144 144 ^ 7 2304 

2357 ^xBr^ 8309 SqSs 109 ^SqSs 69973 Si Sg 

2304 44 44 ^ 1188 ^ 

973^^51^8 17197 SgSg 197 ^^B^B^ 120643 S3 Sg 

1188 132 132 396 

1643v^S3S8 _ 10939 S4S8 _ 139v^S4S8 94717 SgSg 

396 891 891 ~'~ 2376 ~'~ 

1717 ^B^Bs 331907 SeSg 265651 ^/wTSg Sg 

2 376 ^ 63360 ^ 63360 ^ 

3451 ^Af^ ^6^8 331907 S7S8 _ 265651 VwrS7 S8 

63360 ^ 63360 63360 ^ 

2507^ Br Bs 3A51 s/WTw^ Br Bs 4062433 Sg^ 

63360 63360 ^ 580800 ^ 

169699 y^Sg^ 8309 So Sg 109^SoS9 2507 ^B^B^ 

1742400 44 ^ 44 ^ 63360 ^ 

69973 Si Sg 973^/^SiS9 17197 S2S9 197 ^B2Bg 

1188 1188 132 ^ 132 

120643 S3 Sg 1643y^S3S9 10939 S4S9 139 ^B^Bg 

396 396 891 ~'~ 891 

94717 B^Bg 1717 ^SsSg 331907 SeSg 2507 ^^SgSg 

2376 2376 ^ 63360 63360 ^ 

265651 ^/u^SgSg 3451^/^iJru^S6Sg 331907 S7S9 

63360 63360 ^ 63360 
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2507 57^9 265651 Juh B7 Bg 4062433 5,2 



+ 

6 3360 63360 580800 

3451711^^5759 695 58 59 169699 7^^592 

63360 ^ 1188 1742400 ' ^ ^ 

where Wi = 97 and W2 = 4801. 

Particularly useful expressions are obtained by recasting these formulae in terms of strict 
moments ^^^^"^"^^ = [riiO n^) via Eq. (|90D . The results up to order M = 2 (in correspondence 
to the nucleon distribution amplitudes on the orbit) are 

Ip = 196 4899 - 9851 (001) + 47231 (001)^ + 24954 (002) - 163098 (001) (002) 

+ 138696 (002)2 _ 39425 (iqo) + 95825 (001)(100) - 175500 (002)(100) + 67625 (100)^ 
+20790 (101) + 30720 (001)(101) - 30780 (002)(101) - 22350 (100)(101) 
+37800 (101)2 + 29310 (200) - 132120 (001)(200) + 231780 (002)(200) 

-161550 (100)(200) + 101400(200)2 , (B3) 



490 -485 - 982 (001) + 3799 (001)2 _ ^^^S (002) + 30300 (001) (002) - 26544 (002) 

+4792 (100) - 16184 (001)(100) + 34176 (002)(100) - 12680 (100)2 - 1728 (101) 
-6744 (001)(101) + 5976 (002)(101) + 1200 (100)(101) - 3420 (101)2 - 5400 (200) 
+26880 (001)(200) - 44520 (002)(200) + 30720 (100)(200) + 2520 (101)(200) 

-19020 (200)2 



(B4) 



Similar expressions can be obtained also for the other nucleon form factors displayed below. 

The axial form factor of the nucleon in terms of nucleon eigenfunctions, including next- 
to-leading (M = 3) contributions, is 



V 3 
175 5i 52 
9 

9625 532 _ 

108 
875 5i 55 

108 
125 5o56 _ 

192 
1195 52 Be 

576 
35 54 56 

144 
31891502 

18432 



N\ 



75 5o2 175 Bo B 



25 5i2 25 5o 52 

+ 

108 2 



4 6 

^ 75 522 ^ 625 Bo B^ _^ 1375 B1B3 _^ 125 52 53 _^ 



4 



6 



125 Bo 54 275 52 54 



27 

575 53 54 



625 Bo 55 



625 52 55 
108 

275^9750 Be 
192 

635^9752 Be 



54 

2875 53 55 



54 

10255455 



36 

15625 5. 



108 

715 5i56 _ 

1152 
3265 53 56 



972 3888 
265^9751 56 _ 

384 

4145^9753 56 



+ 



65 



576 
9754 56 



+ 



432 
^99A/97 5fi2 



1728 
95 55 56 985 

384 ^ 
125 Bo Bj 



1728 
9755 56 ^ 



18432 



192 



+ 



3456 
275^97 50 57 

192 
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715S1S7 ^ 265 V975i57 1195 ^2^7 , 635 V97 52^7 



1152 384 576 576 

3265 53^7 4145 79753^7 35 BiBj 65 VOTB^Bj 



+ + 



1728 1728 144 432 

95 55^7 985 ^9755^7 11356^7 31891 ^7^ 



384 3456 36 18432 

899^97 ^7^ , 8309 Bo Bs ^ 109 V^SOl Bo Bs ^ 17 Bi Bs 



18432 2112 2 112 792 

y/miBiBs 1719752^8 197^480152 58 

2376 ^ 6336 ^ 6336 ^ 

120643 53 58 1643 V480T53 58 9471755 58 



19 008 19008 114048 

1717V480T55 58 33190756 58 265651 56 58 



114048 3041280 3041280 

2507V480T56 58 3451 Bq Bg 33190757 58 



3041280 3041280 3041280 

265651^9757 58 2507^480157 58 3451 ^/iFIlI^57 58 



3041280 3041280 3041280 

3701 58^ ^ 17 V480T58^ ^ 8309 Bq Bg 109 V480T5o Bg 



2613600 871200 2112 2112 

UBiBg VA801 BiBg 17197 BiBg 197^4801 52 59 

792 2376 ^ 6336 6336 ^ 

120643 53 59 1643 V48OI5359 94717 B^Bg 

19 008 19008 114048 ^ 

1717V480T55 59 331907 56 59 265651 56 59^ 



114048 3041280 3041280 

2507^480156 59 3451,/w^BeBg 33190757 59 



+ 



3041280 3041280 3041280 

265651 ^9757 59 2507^480157 59 3451^/^^111^57 59 



3041280 3041280 3041280 

2 

9 



58 59 3701 59^ 17^4801 5^ 
1188 ^ 2613600 871200 



This expression verifies the leading-order results obtained by Carlson and Poor [ 126 
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